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American Society for X-Ray and Electron Diffraction 


N December 31, 1941, the American Society 

for X-Ray and Electron Diffraction voted ac- 
ceptance of an invitation extended by the American 
Institute of Physics to become an Associated Society 
of the Institute. This date not only marks the thresh- 
old of a new year but the beginning of a new type 
of relationship between the Institute and a new 
society operating in a relatively new field of physical 
science. 

The American Institute of Physics was founded 
in 1931 as a cooperative effort on the part of four 
societies active in various fields of physics and was 
shortly thereafter joined by a fifth. These five so- 
cieties—The American Physical Society, The Optical 
Society of America, The Acoustical Society of Amer- 
ica, The Society of Rheology, and The American 
Association of Physics Teachers—comprise its pres- 
ent membership. Among its functions is that of 
recognizing and encouraging the growth of new 
fields related to the parent science, Physics. Its 





structure permits extension into various specific 
helds of physics and into other fields which lie on 
the borderline between physics and other sciences, 
still preserving at the same time a proper sense of 
inter-relationship among these divergent fields. This 
policy of divergence without disunity is of funda- 
mental importance in the maintenance of a recog- 
nizable science of physics. 

The relationship of the American Society for 
X-Ray and Electron Diffraction to the Institute is 
that of ‘Associated Society.’’ This type of relation- 
ship was provided for in the by-laws of the Institute 
in anticipation of the organization of new societies 
in physics and closely allied fields. It is particularly 
well suited for such cases since, while expressing the 


fundamental unity of physics and related sciences, 
the relationship carries with it no onerous financial 
burden which might prove disastrous to a society in 
the early stages of its development. 

The American Society for X-Ray and Electron 
Diffraction came into formal existence on July 30, 
1941, through the adoption of a constitution and 
by-laws. The charter membership of 135 has since 
increased to over 200. In addition to this first official 
meeting, held at Gibson Island, Maryland, a meeting 
was held in December, 1941, at Boston. A feature of 
this meeting was a joint session with the Mineralogi- 
cal Society of America. Some twenty papers dealing 
with various aspects of x-ray and electron diffrac- 
tion, particularly in their relationship to the struc- 
ture of minerals, were presented. The society plans 
to hold two meetings each year one of which will be 
a joint meeting with some other society or group of 
societies. The desirability of such a joint meeting 
plan becomes apparent from a consideration of the 
close relationship between diffraction, methods and 
other fields of work such as chemistry, mineralogy 
and various branches of physics. 

The high caliber of the men responsible for leader- 
ship in the organization is indicative of what may be 
expected of it. Dr. Maurice Huggins was its first 
president and was succeeded by Dr. B. E. Warren. 
Dr. M. J. Buerger and Dr. George Tunell are the 
present vice president and secretary-treasurer, re- 
spectively. 

The American Institute of Physics wishes to take 
this opportunity to express its satisfaction in this 
new affiliation and to welcome the American So- 
ciety for X-Ray and Electron Diffraction to associa- 
tion with the present Institute membership. 


H. A. B. anno G. H. B. 
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X-Ray Absorption Structure as a Tool in 
Physics and Chemistry 


BY S. T. STEPHENSON 


Washington State College, Pullman, Washington 


—RAY absorption edges are explained in the 

following manner: As the wave-length of 
the x-rays passing through the absorber becomes 
shorter the energy of the quanta becomes great 
enough to eject one of the inner electrons of an 
atom in the absorber. The absorption coefficient, 
therefore, has a sharp discontinuity for this 
particular energy (wave-length) and an absorp- 
tion edge results. The K x-ray absorption edge 
for an isolated atom is due to absorption of the 
energy of an x-ray quantum by a K electron 
which jumps to the lowest unoccupied outer 
level of the atom. If the atom has alternate 
forbidden and allowed optical levels extending 
above this lowest allowed level, there will be 
variations in the probabilities of the ejection of 
electrons with energies larger than that necessary 
to jump to the lowest unoccupied level. Thus a 
fine structure will be found on the short wave- 
length side of the main edge! which can extend 
out no further than the energy corresponding to 


complete removal from the influence of the atom. 


(This will be about equal to the ionization 
potential of an atom whose atomic number is 
one greater.) Argon,? for example, shows a fine 
structure of this type extending out some five 
volts beyond the main edge. 

If the atom from which the K electron is 
ejected is part of a molecule, it can be shown’ 
that scattering of the ejected electron wave by 
the other portions of the molecule and subse- 
quent interference lead to the origin of energy 
regions of varying probability which are char- 
acteristic of the molecule. These energy regions 
extend sometimes a hundred volts or more above 
the lowest unoccupied level. Thus one expects 
for molecules an extended structure (any struc- 
ture due to energy zones outside of the atomic 
optical levels will be called extended) on the 
short wave-length side of the edge which may 
run out for a hundred volts or more. GeCl,! 
and Br>-* show such prominent structures. 
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If the K electron is removed from an atom in 
a crystal lattice it is scattered by the surrounding 
atoms of the lattice. The resulting interference 
leads to allowed and forbidden energy zones for 
the electron extending out,® in some cases, 
several hundred volts above the lowest unoccu- 
pied atomic level. It turns out that the forbidden 
energies of the electron are those for which the 
electron, upon removal from the atom in the 
lattice, falls on the crystal planes in such a way 
as to satisfy the simple Bragg condition (if the 
energies are high enough so that the electron 
may be considered to be free). Metal crystals, 
in particular, show a marked extended structure, 
e.g. Fe,!° Cu," Zn," and brass." The electrons 
removed by the absorption process in ionic 
crystals are subject to the electrostatic forces of 
the parent ion and to scattering by the sur- 
rounding atoms in the lattice, both of which 
affect the extended fine structure to a marked 
degree.'4—!6 

All the investigations of extended fine struc- 
tures to date have given results in agreement 
with the following general picture of the origin 
of the structure: 

(1) Electrons are ejected from an atom or ion 
with a definite energy and wave-length. 

(2) These electrons are scattered by the sur- 
rounding nearest neighbors of the parent atom 
or ion (if any). 

(3) The resulting interference pattern for the 
electron waves upon averaging over all directions 
of ejection has a larger amplitude for some 
energies than for others. 

(4) Thus, there is a greater probability for the 
ejection of electrons of some energies than for 
others. 

(5) Therefore, there is a variation in the 
probability of the absorption of x-ray quanta of 
different energies, resulting in extended fine 
structure. Essentially, then, extended x-ray 
absorption fine structure is an electron diffraction 
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phenomenon. The function of the x-rays is to 
remove the electrons from an inner atomic shell. 
One infers the electron diffraction pattern from 
the fluctuations in the intensity of the x-ray 
quanta having energies within a few hundred 
volts of the main absorption edge. 

Most experiments on extended structure up 
to the present time have been performed for the 
purpose of testing theories regarding the origin 
of the extended structure. It is the purpose of 
this paper to emphasize that if this process is 
reversed, i.e. if one adopts the present picture of 
the origin of the extended structure and uses it 
in designing new experiments to investigate 
certain problems in physics and chemistry, one 
has a powerful research tool. Electron diffraction 
alone offers a very sensitive method for studying 
certain problems but it is restricted in use 
because of the lack of penetration of the elec- 
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Fic. 1. Absorption structure of bromine (reference 5) 
in various compounds. The lower curve illustrates the 
very pronounced structure observed for bromine vapor in 
molecular form. The two curves in the middle are for the 
vapors of ethyl and methyl bromide. Note the absence 
of any structure. Apparently the H atoms and the C atoms 
scatter the electron wave to a negligible degree. The 
curves then resemble those for isolated atoms. The top 
curve is for solid silver bromide. 


trons. X-ray extended structure measurements 
combine the sensitivity of electron diffraction 
methods with the penetration of x-rays. Many 
problems in which the positions of the nearest 
neighbors of an atom or ion are required may 
be solved by this method. 

In the following sections are outlined some 
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applications of extended fine structure measure- 
ments which have been made or could be made 
to problems in physics and chemistry. 


GASES 


The determination of the arrangement of the 
component parts of a molecule in a gas or vapor 
is a problem on whose solution observations of 
extended fine structure may throw considerable 
light. Such observations offer an independent 
means of checking and extending the knowledge 
already gained from x-ray scattering, spectro- 
scopic data, and chemical observations. Very 
little work has been done in this field, however. 
Br is the diatomic molecule whose structure has 
been most extensively studied.5~* Snyder and 
Shaw,* in a theoretical discussion of the Bre 
structure, conclude that the fundamental idea 
of electron scattering is predominant in pro- 
ducing the observed structure. However, the 
electron scattering is affected by such things as 
molecular vibration, exchange and multiple 
scattering, and polarization of the scattering 
atom, all of which are hard to include in any 
theoretical calculations. 

Of the polyatomic molecules GeCl, has the 
best known structure.'”* A theoretical discussion 
from the standpoint of electron scattering by 
Hartree, Kronig, and Petersen!’ yields good 
agreement with experiment and these authors 
are convinced that interatomic distances, the 
character of the molecular bond, and the energy 
states of the atoms in the molecule may be 
determined by this method. By observing 
extended structures for GeCl, and AsCl; Coster 
and Klamer were able to determine the distances 
of Cl atoms from the Ge and As atoms. 

Bogdanovich"® obtained expressions for the 
extended structure expected to be observed for 
triatomic molecules and predicted that the 
isomeric configuration of such molecules could 
be determined. Stelling?® has actually studied by 
experiment the dependence of the main edge 
and the extended structure on the distribution 
of the atoms in certain complex molecules. He 
has discovered some rather broad empirical 
rules enabling him to determine the structure 
of other complex molecules. 

To date the situation for gases seems to be 
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this: If the observed extended structure of a 
molecule is known, it is still not theoretically 
possible to determine with precision the actual 
shape and configuration of the molecule from 
these data alone. However, by observing the 
extended structures of molecules of different 
types new information regarding those molecules 
is made available to correlate with information 
obtained from other sources. Furthermore, by 
comparing extended structures of molecules of 
unknown type with those of known type definite 
predictions as to the actual shape and configura- 
tion of the unknown molecules may be obtained. 


SOLUTIONS 


The determination of the size and shape of 
molecules or ions in a solution is of interest and 
importance. Extended structure measurements 
play the same role here as in the case of a gas. 
The x-rays penetrate the solvent molecules and 
the electrons are ejected from the atoms under 
study in particular molecules or ionic configura- 
tions as for a gas. The electrons are scattered by 
the immediate neighbors in the molecules or 
ions and a characteristic extended structure 
results. 

One problem of major importance in solutions 
of salts, acids, and bases is the determination of 
the degree of ionic dissociation which takes place. 
The physical chemist concludes, on the one 





hand, from data on osmotic pressures, melting 
points, and conductivities that a salt solution 
contains both molecules and ions and is only 
partially ionized. On the other hand, the physical 
chemist concludes from the fact that the salt 
crystal lattice is entirely composed of ions that 
a salt in solution is likewise completely ionized. 
Debye and Huckel have shown that the data on 
osmotic pressures, melting points, and conduc- 
tivities for dilute solutions of strong electrolytes 
are consistent with the latter view if account is 
taken of the electrostatic forces between ions. 
It is important to see from independent evidence 
such as x-ray absorption extended structure 
measurements whether molecules do exist in 
solutions of strong electrolytes. The observations 
and interpretations are clear cut. If there are 
no molecules present, no extended structure 
would be observed because the free ions could 
have only a very fine close-in structure similar 
to that of atoms in a monatomic gas or vapor. 
(Any effect of the surrounding polarized water 
molecules would be small.) If an extended 
structure is observed, then there are molecules 
present and an estimate of the degree of dissocia- 
tion can be made from the observed intensity of 
the structure, especially if this structure is 
compared with that furnished by molecules of 
the salt in the vaporized state. 

The only experiments with solutions carried 
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out from this point of view are those of Ray, 
Das, and Bagchi! on certain cobalt salts which 
indicate a very small degree of ionization for 
solutions as dilute as N/20. However, some of 
their conclusions are open to question. Work 
is progressing in this laboratory on solutions of 
several salts. 


LIQUIDS 


Most of our knowledge of the structure of the 
liquid state comes from x-ray scattering experi- 
ments. Much added information could be 
furnished by extended x-ray absorption structure 
measurements, which combine the penetration 
of x-rays with the higher sensitivity of electron 
scattering. The ejected electrons will be scattered 
by the surrounding nearest neighbors in the 
liquid. Because of the rather definite distance of 
nearest approach of such neighbors, some inter- 
ference effects are to be expected giving rise to 
a variation in probability of ejection of electrons 
of different energies and to an extended structure. 
More specific information regarding the con- 
figuration of nearest neighbors in the liquid can 
thus be obtained. 

This field is wide open at the present time for 
no measurements have been made to date on the 
extended structure” of the elements in liquid 
form. In fact the only measurements on any 
salt or base in the fused state, of which the 
author is aware, are those of Drynski and 
Smoluchowski! on liquid GeCl,. They found the 
extended structure of liquid GeCl, to resemble 
that for gaseous GeCl, and to be quite different 
from the extended structure presented by solid 
GeCl,. However, they did not attempt to 
interpret their results in terms of a definite 
liquid configuration for GeCly. Some work is in 
progress in this laboratory on fused salts and 
liquid elements and has proceeded far enough to 
demonstrate that these fused materials have an 
extended structure characteristic of the liquid 
state. 


SOLIDS 


Practically all of the work with extended 
structure measurements” has been on_ the 
absorption by solids; chiefly metals, alloys, and 
ionic compounds, all in polycrystalline form. 
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Almost always the emphasis has been to use 
known properties of the solid to explain the 
extended structure. However, extended structure 
measurements can now be made to yield infor- 
mation regarding the solid state. The most 


























Fic. 3. Extended fine structure of solid copper. The 
maximum e is 270 volts from the main edge illustrating 
the large range of extension of energy zones in metals. 


obvious use is in the determination of the energy 
zones of metals and alloys as well as in a study 
of the conduction bands. High resolution is 
required for this work. The experiments of 
Bearden* and his co-workers and of Parratt®® 
should be cited in this connection. 

Since the orientation and configuration of 
solids can be determined so well by means of 
x-ray scattering and diffraction methods, little 
information can be added by using the more 
sensitive electron scattering method except in 
solving such a problem as that of placing 
molecules of water in water of crystallization. 


SUMMARY 
This paper is not intended to be a compre- 
hensive review of the x-ray absorption fine 
structure field. It is pointed out that almost all 
extended structure measurements of the past 
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have been for the purpose of testing some theory 


concerning the origin of such structure. 
Furthermore, it is shown that if this emphasis 
is changed and we assume that the structure is 





caused by an electron scattering process, and. if 
we design and perform the proper experiments, 
a great deal of helpful insight can be thrown on 
several problems in physics and chemistry. 
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~ Relation Between Viscosity of Solutions and Physical 
Properties of High Polymers* 


BY ROBERT SIMHA 


Department of Chemistry, Polytechnic Institute, Brooklyn, New Yorkt 


INTRODUCTION 


HE investigation of high molecular weight 

substances by physical methods is based 
fundamentally, apart from such direct attacks 
on the structural problem as x-ray examination, 
on the ability of the dissolved molecules to 
generate an effect in the solution dependent upon 
their number or weight and size. Since cryoscopic 
effects are too small, as soon as we exceed mo- 
lecular weights of about 10,000, the means of 
investigation applicable are sedimentation and 
diffusion analysis by equilibrium and non-equi- 
librium methods, osmotic pressure determina- 
tions, and viscosity measurements. In this paper 
we shall deal mostly with the last method, dis- 
cussing results found in other ways only for 
comparative purposes. 

We may therefore ask ourselves: What can be 
said about the size and shape of the molecules of 
a high polymer on the basis of the change in 
viscosity of the solvent, brought about by dis- 
solving the high polymer in it? Several factors 
contribute to making a complete answer to this 
question difficult. At this point we need re- 
member only one factor very fundamental in the 
whole work on giant molecules: the necessity to 
interpret our results statistically, as in a gas or 
liquid, for two reasons: first, because of the 
large number of degrees of freedom possessed by 
a macromolecule with its various possibilities of 
intramolecular arrangement, which can be formu- 
lated only statistically; second, because of the 
fact that synthetic and natural materials, even 
when fractionated, show an inhomogeneity in re- 
spect to molecular size. 

1. THEORY 

For the quantitative treatment of the viscosity 
of solutions, the hydrodynamic method has 

* Presented at the Annual Meeting of the Society of 
Rheology, October 24 and 25, New York City, 1941. 


+ This work was started while the author held a Lalor 
Foundation Award at Columbia University, New York, 
® 2 
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proved so far to be the most effective, at least in 
the range of sufficient dilution. Its validity rests 
on the condition that the dimensions of the 
solute molecules are large compared with those 
of the solvent. This assumption is needed in 
order to treat the solvent as a continuum obeying 
the Stokes-Navier equations of motion. A second 
condition is that the solvent should be a New- 
tonian-fluid. In this theory then, the influence of 
the solute on the viscosity of the pure solvent is 
measured by the disturbance of the original flow 
of the solvent in the viscosimeter, due to the 
presence of foreign molecules. 

The relative viscosity increase will depend 
primarily upon the concentration. As first shown 
by Einstein' the volume of the dispersed com- 
ponent is the determining factor. We may expect 
for small values of the volume concentration ¢ an 
expression of the type: 


7 ‘no— 1 = Ns p= aie +420? +a 3c*+ --- (1) 


(nm viscosity of the solution, 7 viscosity of the 
solvent). The determination of these coefficients 
a; is the task of the theory. The various empirical 
equations formulated can all be expanded into 
such a power series.” The physical significance of 
the coefficients a2, a3 is that of interaction terms. 
In our hydrodynamic picture they measure the 
influence of the movement of one dissolved 
particle in the streaming fluid on the motion of 
its neighbors. In the case of spherical suspensions 
the a; are pure numbers. Up to the second term 
they have been calculated by means of a pertur- 
bation method, a,;=2.5, a2=14.1, by Ejinstein,' 
Guth, Simha, and Gold, respectively.* A factor 
which frequently has not been taken into account 
by experimenters is the type of viscosimeter 
used. At first all pertinent results were derived 
for the behavior of a solution subjected to a 
constant shear. An additional investigation by 
the author‘ showed that, in the linear concentra- 
tion range, solutions of sphere-like molecules 
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behave alike in a rotating cylinder and a capillary 
tube apparatus. Experiments by Eirich and his 
collaborators’ with models like glass spheres 
seem to confirm this fact. In other cases, how- 
ever, there is no justification for extrapolating 
results found in one apparatus to another one. 

Colloidal particles of spherical shape are very 
seldom encountered, thus prohibiting an applica- 
tion of the above relatively simple theory. In 
order to arrive at quantitative relations between 
viscosity increase and dimensions of chain mole- 
cules, we may approximate a long chain by an 
ellipsoid of revolution, a cylindrical rod, or by an 
array of spheres. As long as these idealized 
molecules do not change their shape under the 
influence of the shear gradient, the exact shape 
will not be important and will show up only in 
slight changes of numerical factors, if the mole- 
cules are long enough, or more exactly, if their 
ratio between longest and shortest axes is large 
enough. Due to this asymmetrical shape, we may 
expect a considerable dependence of the results 
on the degree of orientation in the solution as, 
for instance, measured by the birefringence of 
flow. The dependence is the more pronounced the 
greater the asymmetry. This is just what the 
theory predicts. Figure 1 shows the values of the 
coefficient a, in the above equation as a function 
of the axis ratio f for two extreme cases. The first 
one is that of complete orientation in the plane 
of flow, according to G. B. Jeffery.* The second 
curve corresponds to the other extreme of com- 
pletely random orientation.’ In this case the 
increase is much steeper. For f-values larger than 
about 50, the ratio of the coefficients for the two 
cases is approximately proportional to f. The 
case of complete orientation is, of course, never 
realized, except in experiments on models such 
as those carried out by Taylor* and more recently 
by Eirich and his collaborators.® In general, 
there will be partial orientation corresponding to 
a finite value of the ratio of velocity gradient to 
rotatory diffusion coefficient® 

«/ Dror. 

The viscosity increases with decreasing values of 
this factor. The smaller the gradient, and the 
greater the intensity of Brownian movement, the 
higher the viscosity of the solution. Evidently 
because of this ratio the results will largely 
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Fic. 1. Viscosity coefficient, a:, as function of axial 
ratio, f. 1: Without Brownian movement (viscosity scale 
multiplied by ten). Ia: Represents portion of curve | 
below f=10, with viscosity scale multiplied by 100 and 
axial ratio scale multiplied by 10. Il: With Brownian 
movement. Ila: Represents portion of curve II below 


f=10, both scales multiplied by 10. 


depend upon the type of apparatus used, as x is 
variable in a capillary tube apparatus. Extrapo- 
lation to zero velocity gradient, leading to the 
extreme of so-called ‘“‘complete Brownian move- 
ment,’’ should, however, lead to a result inde- 
pendent of the viscosimeter used as the disorien- 
tation tendency renders the peculiarities of the 
special flow type insignificant.” '° It appears 
very desirable to have more experimental ma- 
terial in connection with this point. On the other 
hand at a given flow rate, a larger depression of 
the viscosity due to orientation for longer mole- 
cules with their smaller D,., would result. Data 
of Bauer, Fuller, and Heiss, Jr.'' on polyesters in 
chloroform in capillary tubes are in agreement 
with this conclusion. 

Because of the very strong dependence of the 
coefficients a; upon the asymmetry of the chains, 
the longer the molecules, the smaller is the range 
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of validity of a linear viscosity equation. To 
speak of an absolute range of validity of ideal 
laws therefore has no meaning unless the value 
of the intrinsic viscosity or the molecular weight 
of the material involved is specified. In the case 
of molecular weights of several hundred thou- 
sand, a concentration of 0.1 g/l may already be 
much too high. This is all the more to be con- 
sidered since no theoretical work has been done 
on the viscosity of more concentrated solutions 
of asymmetrical particles. Therefore, conclusions 
on size must be based on values found by ex- 
trapolation to zero concentration. 

The fact that all effects in solution are so much 
more pronounced for asymmetrical particles than 
for spherical ones is, in terms of hydrodynamics, 
simply due to the ability of spherical particles to 
follow the rotational components of the solvent 
motion completely, while the long molecules 
rotate more or less independently in the solution. 
The shear gradient exerts deformational forces 
on an immersed particle. If the particle is able to 
respond to these forces, then the viscosity will 
be reduced as compared with that produced by 
rigid molecules of the same size. On the other 
hand, the more or less free rotation around each 
single bond, present in the long chain, gives rise 
to the familiar curling effect, due to intramolecu- 
lar Brownian motion of parts of the chain.” 
Such molecules are deformable under shear, that 
is, they follow more readily the original flow of 
the solvent, thus reducing the intrinsic viscosity 
of the solution at a given concentration. Further- 
more, the concentration effect, as compared with 
that of rigid rods, is reduced because these soft 
coils have smaller axis ratios. 

Another factor to be taken into account is the 
solvation effect. In. our Eq. (1) the quantity c 
contains the hydrodynamically effective volume 
of a unit moving independently in the solution, 
or influenced merely by forces of hydrodynamic 
origin. The kinetic unit, therefore, comprises the 
dry molecule plus bound solvent. Deviations 
from the laws found for spherical suspensions 
may be due to two causes. First, solvation and 
second, asymmetry. Unfortunately, at present it 
is impossible to separate these two effects on the 
basis of viscosity measurements alone. In the 
case of solutions of approximately spherical 
particles, like gamboge or cane sugar, it is 
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possible to estimate a degree of solvation with 
the aid of Eq. (1), for spheres including quadratic 
terms in c. However, if the solvation degree so 
computed is larger than about 50 percent, then 
the large viscosity effect must be due at least 
partially to the asymmetric shape. It would be 
possible in principle, to determine the axis ratio 
as well as the amount of combination with the 
solvent from a knowledge of the theoretical 
viscosity-concentration curve in the quadratic 
range. As mentioned before, this knowledge is not 
available. We must, therefore, for our purpose 
use a combination of different independent 
methods such as sedimentation-diffusion meas- 
urements or, for proteins, measurements of di- 
electric dispersion. An attempt to separate in 
this way asymmetry and influence of the medium 
has been made by Mehl, Oncley, and the author.’ 
Figure 2 shows the results for two proteins, egg 
albumin and thyroglobulin. The one pair of 
curves gives the various axis ratios (vertical axis) 
which, combined with the corresponding degree 
of solvation, account for the measured specific 
viscosity. There are two curves allowing for an 
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Fic. 2. The asymmetry of protein molecules calculated 
from viscosity and diffusion assuming various degrees of 
hydration. 
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error in experimental values of +10 percent. 
An analogous pair of curves assuming an error 
of +4 percent is derived from sedimentation 
data. The area common to both regions permits 
an estimation of axis ratio and solvation degree, 
or at least of a range of possible values. The set 
of curves in the lower half of the figure was 





TABLE I. 
1/f, elongated ff, flattened 
dif- Vis- Dif- Vis- 
Protein f/fo* a1 fusion cosity fusion cosity 

Egg albumin ie’ ie > Se.” eS ee 
Serum albumin 1.25 65 50 56 54 7.7 
Hemoglobin 1.16 53 3.7 46 39 60 
Amandin 1.28 70 54 60 60 8.5 
Octopus 

hemocyanin 138 90 7.2 73 £82 -11.4 
Gliadin 1.60 14.6 10.9 10.5 13.6 21 
Homarus 

hemocyanin 127 44 $2 55 $38 75 
Helix pom. 

hemocyanin 124 64 48 55 5.2 7.5 
Serum globulin 141 90 76 73 89 11.4 
Thyroglobulin 143 99 138 $I9 92 12.7 
Lactoglobulin > 2 ee” oe? ae oe 
Pepsin 1.08 5.2 25 45 26 5.8 
Helix hemocyanin 

6.6 12.0 23.9 26 


pH 8.6 1.89 18.0 1 


* Mobility ratio as determined from sedimentation-diffusion. 


drawn under the assumption that the molecule 
can be represented by a flattened ellipsoid. As can 
be seen, there is no agreement between the two 
types of results. It may be mentioned that 
dielectric results evaluated by Oncley™ in the 
same way give agreement with this set of data. 
Table I shows numerical results for some ten 
proteins.” The second column contains the asym- 
metry factor as determined from diffusion meas- 
urements. It is the ratio of the translational 
mobilities of the molecule and of the sphere of 
equal volume. The next column contains the 
values of the quantity a;. The fourth and fifth 
columns give the calculated axis ratios. If we 
discard the two last substances, pepsin and 
helix hemocyanin, the agreement is as good as 
can be expected. The last two columns give 
corresponding data for oblate ellipsoids. As be- 
fore, this possibility seems excluded. 


2. RELATION BETWEEN MOLECULAR WEIGHT, 
SHAPE, AND INTERNAL FRICTION 


The molecular weights of the proteins con- 
sidered are relatively well known. If we try to 
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plot these molecular weights versus the corre- 
sponding dimensions, no regularities can be 
detected. If we now turn our attention to chain- 
like polymers, containing a large number of 
repeating units, it is evident that there must be 
a distinct correlation between molecular weight 
and e.g. axis ratio or length. The hydrodynamic 
theory, on the other hand, gives a viscosity 
increment as function of the dimensions of the 
solute molecule only. A non-phenomenological 
approach would, of course, contain as one of the 
parameters the molecular weight. No complete 
theory of this kind exists as yet. 

The very important question of the determina- 
tion of the molecular weight of a high polymer 
together with the ease with which viscosity 
measurements can be carried out for this purpose 
make imperative an attempt to find a connection 
between internal friction of the solution and 
molecular weight of the solute. The problem is 
complicated by the previously-discussed fact 
that macromolecular systems, containing very 
asymmetric molecules, are inhomogeneous in re- 
spect to size, usually being more inhomogeneous, 
the bigger the molecules on the average. Several 
approaches to this problem have been made. 
They are, as to be expected, of a theoretical, 
empirical and semi-empirical nature. 

An attempt of the first kind is that of Kuhn, 
Guth, and Mark.” In their theory, the average 
length and diagonal dimensions of an isolated 
chain with given bond distance and bond angle, 
kinked at random because of free rotation, are 
calculated on the basis of statistical considera- 
tions. The length, for instance, is found to be 
proportional to the square root of the degree of 
polymerization. This leads, as shown in particular 
by Huggins,® to direct proportionality between 
specific viscosity and molecular weight, if the 
concentration c, is expressed as weight per 
volume, and if the molecules are sufficiently long. 
Unfortunately, the difficulties are such that the 
quantitative results are still open to the intro- 
duction of correction factors, like restricted 
rotation, influence of solvent and others, the 
exact nature and extent of which are not known. 
Purely empirical is the well-known approach of 
Staudinger.'* He found that in certain polymer- 
homologous series, the specific or intrinsic vis- 
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cosity increases in proportion to the molecular 
weight : 
Nsp=Ku Mew. (2) 


The quantity Ay, characteristic for each series, 
is, according to Staudinger, independent of the 
molecular weight. This is in agreement with the 
theoretical result mentioned before. Kraemer and 
Lansing'® were the first to point out that a 
molecular weight determined by means of this 
rule represents a weight average, which is 
identical in type with the one found in the 
ultracentrifuge, while the osmotic method gives 
a number average value quite different from the 
other one in polydispersed systems. It is therefore 
possible, by combining osmotic and viscosity (or 
ultracentrifuge) data, to gain information about 
the size distribution curve, large discrepancies 
being indicative for a broad distribution. 
Fordyce and Hibbert’s'* homogeneous polyoxy- 
ethylene glycols of relatively low molecular 
weights, however, obey in the range of 800-8200 





TABLE IIT. 
M P Km104 Literature 
Cellulose triacetate in m-cresol 
1542 5 10.3 ) 
Ber. Dtsch. chem. Ges. 67, 
22500 79 «= «6.7 |_«S 479 (1934). 
50000 175 57 
100000 370 6.65 
H. Lachs und A. I. Gros- 
- man, Bull. Acad. Polonaise 
| 3A, 11, 220 (1939). 
272000 900 3.6 


Polyvinyl Acetate in Acetone 
75000 870 0.52 


——— oa on | J. prakt. Chem. 158, 278 


315000 3650 0.26 | (1940). 
640000 7400 0.25 | 


Methyl Polymethacrylate in Acetone 
25000 250 0.56 | 


54000 540 0.43 | ai . 
113000 1130 0.35 | J. prakt. Chem. 155, 292 


158000 1580 0.25| (1940). 
260000 2600.19. 


a linear relationship with an additional constant. 
As a consequence, the Staudinger constant de- 
pends upon the molecular weight. Flory and 
Stickney’s work" on inhomogeneous polyesters in 
the range of 1500-30,000 gives a similar result." 
Recently K. H. Meyer’ collected very extensive 
material on various high polymers of widely 
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differing molecular weights. Table II shows some 
of his data. The decrease in the Ky constant 
is, in some of the series, rather marked. The 
molecular weights of the various fractions were 
determined by independent methods. In spite of 
these discrepancies, the rule is being used ex- 
tensively for a rapid relative estimation of mo- 
lecular weights. It is very useful in this respect 
when calibrated by means of other methods, as 
done by Kraemer and Lansing," and Signer and 
his collaborators.'® 

A third attempt of a semi-empirical character 
could give a relative measure for the internal 
flexibility of the molecule and in the case of 
branched polymerization products, like poly- 
styrene at higher reaction temperatures, a rela- 
tive measure of the extent of branching, if the 
distribution curve of molecular weights is known. 
The average relation between axis ratio f and 
molecular weight M will be for sufficiently long 
molecules of the nature: 


f=C-M?, (3) 


where the constants C and p are characteristic 
for a given series. In the case of a rigid, linear 
rodlike molecule we have p=1, while in the 
actual molecule we may expect values smaller 
than one, varying with the degree of curling and 
with the extent of branching. This equation 
holds for each species in the mixture. The in- 
homogeneity of a system is usually described in 
terms of molecular weights rather than mo- 
lecular dimensions. Kraemer® has extended the 
previously described viscosity results of the 
author’ to non-homogeneous systems by repre- 
senting the axial ratio distribution by a function 
identical in type with that used by him and 
others to represent the molecular weight distri- 
bution of a fractionated product. Of course, any 
distribution can be used for this purpose. Now 
with the aid of the above relation and the weight 
distribution parameters, those of the f-distribu- 
tion can be found. These data, combined with 
viscosity measurements of several members of 
the series, give the characteristic constants C 
and p. It is hoped to apply this procedure to 
various high polymers at a later stage. Schulz 
and Jirgensons* derived semiempirically a rela- 
tion between the precipitability of a high polymer 
and its molecular weight, in which the negative 
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exponent of M is equal to or smaller than one, 
varying in the same manner as our coefficient p. 
It would be interesting to compare the exponents 
found from solubility and viscosity measure- 
ments, respectively. 

Finally, applications of the theory to some 
measurements of polystyrene and methyl cellu- 
lose by Staudinger, Signer and collaborators! '® 
may be given. These results were used because 
average molecular weights of several of the 
higher fractions were determined in this work by 
means of careful rate and equilibrium measure- 
ments in the ultracentrifuge. Furthermore, 
Burgers® made use of the same results in his 
report, although arriving at differing quantita- 
tive conclusions. The correct procedure would, of 
course, be to use distributions instead of average 
values as outlined above. However, for this com- 
parison, it may suffice to use average values. 
Table III shows the molecular weight, the vis- 
cosity ratio, and the axis ratio /]/d following 
therefrom. Assuming a partial specific volume of 
0.92 for polystyrene independently of the molecu- 
lar weight, the absolute values of length and 
thickness can be found. The last column shows 


TABLE III. 
Size and Shape of Some Polystyrene Molecules 
ul nep/€ l/d lA lmax4 
519 4.00 3.3 25 7.7 12.5 
2460 7.64 6.4 67 10.4 60 
5225 11.84 9.0 108 12.0 125 
7600 14.78 10.6 136 12.8 182 
13000 23.70 14.6 200 13.7 312 
30000 28.15 16.5 289 17.5 720 
80000 58.70 26.3 544 20.7 1922 
270000 250.00 61.0 1427 23.4 6490 


I/d~O12M2, 1~0.34M?3, d~2.9M"%, partial specif. 
volume v=0.92. 


Size and Shape of Some Methylcellulose Molecules 


VU 1 l/d lA dA lmax4 
14100 0.683 32 313 9.9 410 
24300 .720 45 486 10.8 670 
38100 731 58 665 11.3 1040 


L/d =0.09 M*'8, 1~0.26M*"*, d~3.0M"8. 


the values of the maximum length, calculated 


from Staudinger’s value of 2.5A for the unit. If 


we discard the first rows, we see an increasing 
degree of curling with increasing molecular 
weight, as one might expect. Approximate rela- 
tionships between molecular weight and dimen- 
sions according to (3) can be derived. These val- 
ues do not lead to Staudinger’s law. Below are 
data for methylcellulose in water. However, they 
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extend over a small weight range only. Again the 
molecules have a smaller length than could be ex- 
pected on a purely geometrical basis. The effect, 
however, is not as pronounced as in the former 
case. This can also be seen from the values of the 


TABLE IV. 


Sedimentation Constants 








(a) of some Polystyrene Molecules: 


] — §-10'3 (obs.) — S-10% (calc.) 
30000 3.8 3.7 
80000 6.5 6.0 
270000 16.5 9.6 
(b) of some Methylcellulose Molecules: 
J — S-10" (obs.) — S-10"% (cale.) 
14100 0.83 1.03 
24300 0.79 1.10 
38100 0.89 1.28 


exponents (3 and ?), indicating a somewhat stiffer 
molecule. In Table IV are some sedimentation 
constants calculated from the data of the previ- 
ous table. The center column shows the measured 
values extrapolated to zero concentration. As far 
as one can judge from a few data, the agreement 
becomes worse for polystyrene, the higher the 
molecular weight. One way of accounting for this 
discrepancy is the following :** In the viscosimeter 
and especially in a capillary tube, the shear 
gradient may extend a flexible coil somewhat, as 
pointed out by Burgers. The observed sedi- 
mentation constant would become greater than 
the calculated one, which depends inversely upon 
the length, because the stretching effect is 
negligible in the ultracentrifuge. In the case of 
methyl cellulose, we have the opposite behavior. 
Calculated values are always large. This may be 
indicative for a concentration effect which lowers 
the rate. Because of the difference in the type of 
hydrodynamic flow such interactions will start 
in sedimentation experiments at concentrations 
lower than the corresponding ones in a viscosim- 
eter. Also the approximation of using average 
values may be responsible. However, these con- 
siderations are at present hardly more than a 
guess.”* 

The author wishes to express his indebtedness 
to the Lalor Foundation for a Fellowship. He 
had the privilege to discuss this work before a 
meeting of members of the Research Labora- 
tories of the Dow Chemical Company, Midland, 
Michigan. 
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3. CONCLUSION 


As a summary, we may say that it is possible 
today to formulate quantitative relations be- 
tween shape factors of the solute molecule and 
viscosity increment, in a linear concentration 
range, as long as appreciable orientation effects 
are absent. This can be done in ascribing to the 
molecule an effective shape corresponding to an 
ellipsoid, or as done by Kuhn and Huggins® to 
an array of spheres, which cannot change their 
shape to a considerable extent while moving in 


the viscosimeter. In these cases, it is possible to 
undertake a comparison with sedimentation data. 
Orientation due to shear as well as flexibility 
tends to decrease the viscosity, the latter factor 
presumably to a lesser degree. It is possible on 
this basis to interpret the Staudinger rule and 
deviations therefrom, and to derive correlations 
between molecular weight and shape from experi- 
ment. Because of imperfections of the theory and 
the scarcity of sufficiently accurate data, in the 
results found independently there are discrep- 
ancies which await further elucidation. 
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Deferment of Academic Physicists and Students 


S a result of the declaration of war, local draft boards throughout the nation have lost sight 
of the paramount importance of occupational deferment of physics students and physics in- 
structors. To correct this situation a directive has recently been issued from the National Head- 
quarters of the Selective Service System (21st Street and C Street N. W., Washington, D. C.). It 
is quoted below so that department heads and other employers of physicists may refer to it when 


writing to local boards. 


January 12, 1942 


MEMORANDUM TO ALL STATE DIRECTORS 
LOCAL BOARD RELEASE 


The attention of local boards is again invited to the 
necessity of seriously considering for deferment students in 
certain specialized professional fields in which dangerously 
low levels of manpower are found to exist. This memo- 
randum is in addition to and does not rescind those previ- 
ously issued which apply to students in other critical fields. 

Subsequent to the declaration of war, local Selective 
Service agencies have in many instances proceeded to 
classify registrants without regard to the fact that they are 
in training or preparation for activities, the maintenance 
of which is essential to the national health, safety, or 
interest, and war production. This is particularly true in 
cases of engineering, chemical, physics, medical and dental 
students. 

Admittedly there is an overlapping of the military and 
civilian requirements of a nation at war; however, it must 
be borne in mind the one is dependent upon the other. It 
is estimated that the expanding Army will eventually re- 
quire doctors and dentists in numbers heretofore unknown. 
They will not be available if those students who show 
reasonable promise of becoming doctors and dentists are 
inducted prior to becoming eligible for commissions. 

War industries are undergoing a hitherto unknown ex- 
pansion, Aeronautical, Civil, Electrical, Chemical, Mining, 


Metallurgical, Mechanical, and Radio Engineers; together 
with Physicists and Chemists are essential to insure a suffi- 
cient flow of material for the armed forces, and industry 
must look to the engineering, chemical, and physics stu- 
dents now in training to meet their present and future 
requirements. 

It is equally important that instructors in these fields 
be seriously consider.d for occupational deferment. Short- 
ages of qualified instructors are known to exist. The edu- 
cational institution employing the instructor should be re- 
quested to file DSS Form 42A in all cases where defe-ment 
is sought. 

In considering student deferment cases, certain local 
boards are requiring the execution of DSS Form 42A in 
addition to the affidavit of the college or university con- 
tained in Bulletin No. 10 issued by the American Council 
on Education. DSS Form 42A should not be required 
where the American Council on Education affidavit has 
been submitted. 

Local Boards will be informed when the manpower re- 
quirements necessary to the national health, safety, or 
interest, and war production become static. Until such 
time, the policy set forth in the Memoranda to All State 
Directors I-62, I-91, I-99, and I-150 remains in force. 

Sincerely yours, 
(Signed) Lewis B. HERSHEY, 
Director 
LBH/gtg 





Calendar of Meetings 


March 
12-13 Society of Automotive Engineers, New York, New York 
17-19 American Railway Engineering Association, Chicago, I}linois 
23-25 American Society of Mechanical Engineers, Houston, Texas 
31—April4 American Physiological Society, Boston, Massachusetts 
April 
3-4 American Mathematical Society, New York, New York 
10-11 American Physical Society, Southeastern Section, Oxford, 
Mississippi 
11 American Mathematical Society, Berkeley, California 
16-17 National Petroleum Association, Cleveland, Ohio 
17-18 American Mathematical Society, Chicago, Illinois 
19 American Ceramic Society, Cincinnati, Ohio 
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April (continued) 


20-24 American Chemical Society, Memphis, Tennessee 
20-24 American Foundrymen’s Association, Cleveland, Ohio 
23-25 American Philosophical Society, Philadelphia, Pennsylvania 


27-29 National Academy of Science, Washington, D. C. 


30-May 2 American Physical Society, Baltimore, Maryland 
May 
4-8 Society of Motion Picture Engineers, Hollywood, California 


11-13 American Institute of Chemical Engineers, Boston, Massa- 
chusetts 
15-16 Acoustical Society of America, Ann Arbor, Michigan 


31-June 5 Society of Automotive 
Springs, West Virginia 


Engineers, White Sulphur 
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New Books 








High Polymeric Reactions, Their Theory and 
Practice 


By H. MARK AND R. Rarr. Pp. 476+-xiii, Figs. 49, 
15233 cm. Interscience Publishers, Inc., New York, 
1941. Price $6.50. 

This book is Volume III of the series of monographs on 
the chemistry, physics, and technology of high polymeric 
substances, the first volume of which consists of a collection 
of the papers of W. H. Carothers on the subject. The 
second volume is entitled Physical Chemistry of High 
Polymeric Systems and is written by H. Mark. The pub- 
lishers indicate that other volumes are to follow. 

The first and general part of this volume consists of: 
(a) a short resume of the structure of high polymers; (b) 
some data on experimental methods, 55 pages; (c) the 
fundamentals of reaction kinetics; and (d) the general 
theory of the mechanism of polyreactions. The last subtitle 
occupies 100 pages. The special section, Part II, of 200 
pages, consists of a discussion of the polymerization of 
various individual substances. Over 100 different organic 
molecules are treated separately as starting materials. The 
polymerization of inorganic compounds occupies less than 
ten pages with ten examples, and polycondensations, which, 
distinguished from polymerization in that H2O, NHz; or 
HCl or other simple molecules are split off in the formation 
of the polymer, are treated in twenty-five pages in the back. 

The field is one in which the classical methods of physical 
chemistry are not completely applicable without some 
modification in point of view. This is primarily due to the 
fact that a “polymer’’ is not an aggregation of identical 
molecules, but a mixture of molecules having many common 
characteristics, with a considerable range of molecular 
weights. Thus no polymer is a pure compound. For in- 
stance, various methods of determination of molecular 
weights lead to different average values such as average 
molecular weight or root mean square; and apparently 
conflicting values obtained by differing methods are not 
necessarily inconsistent but lead to information concerning 
the homogeneity of the sample. The development and 
emphasis of this point of view and its consequences have 
occupied a considerable portion of Volume II of this series. 
The same type of consideration, namely, the extension of 
well-known methods of approach of the classical kineticist 
to a field in which unfamiliar considerations are of primary 
importance, occupies the larger portion of the general part 
of this book. 

Although the general methods of theoretical approach 
to simple reaction kinetics are given in sufficient detail for 
the general reader, the greater portion of the discussion is 
devoted to the complexities introduced by chain or step 
reactions in which the product formed is distributed over 
a range of molecular weights. 

A fundamental distinction is drawn between step reac- 
tions and chain reactions. In the former the initial and 
subsequent steps by which the polymer is built up, are of 
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essentially identical nature and equal rates. This is usually 
the case in the polycondensations for which the essential 
step consists of the combination of two radicals by the 
elimination of some simple molecule such as water. Direct 
polymerization, in which molecules unite by pure addition, 
are generally characterized by chain reactions, in which the 
first step is different in character, and inherently slower, 
than the subsequent growth of the polymer, which is rapid, 
until some special reaction stops the sequence. More em- 
phasis is placed on the discussion of the chain reactions 
than on that of the step reactions. 

The treatment is clear and concise: It is quite readable 
for the lay scientist with no special knowledge of the field, 
and complete enough to be valuable for the specialist. 
There are occasional slips in the mathematical formula- 
tions, but few of these are serious enough to confuse the 
reader. An occasional ambiguity of language is probably a 
failing preferable to a pedantic and difficult style. 

Especially in the discussion of the experimental methods 
and of particular compounds the authors were in the diffi- 
cult danger, inherent in their task, of being either incom- 
plete in their treatment of special cases or boring in their 
repetition of details. The danger has been skilfully avoided, 
in the main. A summary of a large number of methods has 
been given, and the literature references cited should suffice 
to aid the reader seeking precise details to locate the infor- 
mation he needs. The same statement applies to the treat- 
ment of individual compounds. 

The book should be valuable to the specialist as a refer- 
ence book, or to the learner seeking new information. The 
field is one which is becoming of increasing importance, 
not only technologically, but also to the pure scientist. 
This book should do much to aid in uniting the two 
interested groups. 

JoserpH E. MAYER 
Columbia University 


How Things Work 


By GEoRGE RusseL_L Harrison. Pp. 301+ ix, illus- 
trated, 15323 cm. William Morrow & Company, 
New York, 1941. Price $2.75. 

“There are two kinds of people to whom the working 
of such a device as the telephone is not somehow mysteri- 
ous. The first is the person who uses a telephone every day 
without ever stopping to wonder how his voice travels over 
a wire. The second is the person who has wondered, and has 
found out. How Things Work has been written to help 
change the first kind of person into the second kind.’’ Thus 
does Dr. Harrison describe his aim in writing this book. 
The subject matter is not limited to any particular branch 
of natural science, although most of the emphasis is placed 
upon the physical sciences. 

The principal problem confronting the author of such a 
book is the attainment of a text that is comprehensible to 
the youthful reader and at the same time remains within 
reasonable limits of scientific accuracy. The book was ex- 
amined with particular attention to this question, and in 
practically every case where a passage seemed to be either 
not clear or incorrect, the difficulty disappeared upon re- 
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reading. Numerous line sketches are used, for the clarifica- 
tion of the text both by simple illustration and by analogous 
situations within the experience of the intended reader. 


The descriptions of scientific instruments—for example, 
the spectroscope—are excellent. Some of the drawings are a 
little far-fetched, and apparently in an attempt both to 
amuse the reader and to make a single point clear, in- 
accuracies occur in a few of them. Throughout, there is a 
certain animism in the drawings: the representation of the 
sun producing kinetic energy of translation by means of 
hands, of electrons as a crowd of tiny human beings, and 
of birds using propellers to outfly their fellows are examples 
of this questionable technique. Aside from this relatively 
minor point, the book is highly commendable for its ac- 
curacy of text, clear typography, and excellent index. It 
makes pleasant reading and is to be recommended for 
inquiring children. 

D. H. D. RoLLer, EnsiGn, U.S.N.R. 


Diffusion In and Through Solids 


By RicHarp M. Barrer. Pp. 464+xiii, Figs. 158, 
14522 cm. Cambridge University Press and the 
Macmillan Company, New York, 1941. Price $6.50. 

Diffusion in solids is commonly discussed by analogy 
with the classical theory of the flow of heat. When, instead 
of heat, the diffusing solute is material particles which take 
up space and associate themselves chemically with the 
solvent, the phenomena are vastly more diverse and their 
physical interpretation is correspondingly more difficult. 
No previous text, as the author says in the Preface, has 
attempted to correlate and summarize diffusion data in 
condensed phases, save briefly and in relation to either 
chemical kinetics or sorption equilibria. Dr. Barrer is well 
qualified, by reason of his continued interest and sub- 
stantial studies in the field, to undertake a survey. He 
has succeeded in organizing the widely-scattered material 
and presenting it in a book that will be interesting to the 
general reader as well as useful to the specialist. 

The volume begins, properly, with a chapter on the 
formal solutions of the diffusion equation. These are listed, 
for a variety of coordinate systems and boundary condi- 
tions, in explicit form so they may readily be applied to 
experimental data. 

Chapter 2, after a discussion of the several types of flow 
that may occur in single capillaries, describes transport 
through composite capillary systems such as sand masses 
and refractory bodies. Proceeding in the direction of more 
condensed systems, we find in Chapter 3 an account of 
diffusion in zeolites and in glasses. The flow of gases through 
metals, the migration of ions in ionic crystals, and the 
interdiffusion of metals are treated in the next three 
chapters. Factors influencing diffusion constants in these 
systems are summarized, and several proposed models for 
diffusion processes in crystals are criticized in turn—‘*The 
very variety of the attempts which have been made to 
obtain satisfactory models for diffusion, and the limited 
application of so many of them, demonstrate how little prog- 
ress has been made towards a comprehensive treatment.”’ 

The non-equilibrium faults that are inevitable in real 
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crystals and the role of these in transport phenomena are 
the subject of Chapter 7. Chapter 8 discusses migration in 
the surface layer of solids. The two final chapters are de- 
voted to the passage of gases and vapors through mem- 
branes of organic materials. 

More than 100 tables, scattered through the book, 
summarize diffusion and permeating data for various sys- 
tems. The information collected in these tables and the 
formulae of Chapter 7 are alone sufficient to make the book 
valuable as a reference. Citations of original papers are 
copious, averaging 75 per chapter. This fact, together with 
the reflection that the migration of matter through solids 
is a very wide topic to be reviewed in a single volume, 
should console the reader who finds that his own particular 
field of interest has been touched on only lightly. 

R. P. JOHNSON 
General Electric Company 








* 
Here and There 
New Officers 


Dr. Arthur Holly Compton, Charles H. Swift Dis- 
tinguished Service Professor of Physics at the University 
of Chicago and Dean of the Division of Physical Science, 
was elected president of the American Association for the 
Advancement of Science for 1942 at its Dallas meeting. He 
succeeds Dr. Irving Langmuir of the General Electric 
Company. Dr. A. J. Dempster, also of the University of 
Chicago, was elected Chairman of the Physics Section. 


* 


Recent Awards 


Dr. Frederick Ernst Giesecke, Professor Emeritus at the 
Agricultural and Mechanical College of Texas, College 
Station, Texas, received the F. Paul Anderson Gold Medal 
awarded by the American Society of Heating and Ventilat- 
ing Engineers at the 48th Annual Meeting of the Society 
held in Philadelphia in January. The Committee on Award 
honored Dr. Giesecke for his notable contributions to the 
advancement of heating based on his research work in the 
fields of heat transfer and hot water heating, and for his 
eminent services to the Society as member and officer. 


Of particular interest are the following awards of the 
Royal Society: a Royal Medal to Professor Edward Arthur 
Milne, Professor of Mathematics at the University of 
Oxford, for his researches on planetary and stellar atmos- 
pheres, on the internal constitution of the stars, and on the 
theory of relativity; the Hughes Medal to Professor Nevill 
Francis Mott, Professor of Theoretical Physics at the 
University of Bristol, for his application of the principles 
of quantum theory, especially in the fields of nuclear and 
collision theory, theory of metals, and theory of photo- 
graphic emulsions. 
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A Theoretical Approach to the Problem of Yarn Strength 


R. R. SULLIVAN 
The University of Tennessee Agricultural Experiment Station, Knoxville, Tennessee 
(Received May 19, 1941) 


An idealized yarn composed of fibers with specified properties is treated analytically with 
the aim of determining the yarn strength at any degree of twist. The results are presented in 
the form of equations and curves which relate the yarn strength to the fiber properties and 
the degree of twist. Two cases are studied: (1) All fibers alike; (2) fiber properties variable 
from fiber to fiber. In the latter case the mathematical expectation of the yarn strength at 
any yarn cross section is obtained. It is found that the optimum twist multiplier is largely 
determined by the fiber length, fiber fineness, and coefficient of friction, whereas the maximum 
yarn strength (corresponding to the optimum twist multiplier) is more strongly dependent 
upon the intrinsic fiber strength than upon the other fiber properties studied. 


INTRODUCTION 


HE manner in which the physical properties 

of a group of textile fibers contribute to 
the properties of a yarn fabricated from the 
fibers has long been a source of much speculation 
and frequent bafflement to the textile industry. 
Much work has been done, largely of a statistical 
nature, in the attempt to correlate fiber prop- 
erties with yarn properties.'~® It has been found, 
in general, that finer, longer, stronger fibers 
produce stronger varns. The findings of some 
workers are at variance with the findings of 
others. 

The present paper makes no attempt to 
discuss these problems in all their details, but 
simply studies, analytically, the strength of an 
idealized yarn which involves concepts that are 

‘Jerry H. Moore, Agri. Exp. Station of The North 
Carolina State College of Agriculture and Engineering, 
Tech. Bull. No. 58 (1938). 

? Malcolm E. Campbell and Roland L: Lee, Jr., Mimeo- 
graphed Report of U. S. D. A., Agricultural Marketing 
Service and Bureau of Plant Industry (May, 1940). 

’ Malcolm E. Campbell, Address to Am. Soc. of Agron., 
New Orleans, Louisiana (March 23, 1939). 

‘OQ. Johannsen and H. K6éb, Klepzig’s Textl. Zeits. 42, 
333-42 (1939). 

5 C. Underwood, J. Text. Inst. 26, T309 (1935). 

6 Ahmad N. Venkataraman, Proc. 1st Conference Res. 
Workers on Cotton in India, I. C. C. C. (March, 1937). 

7A. J. Turner et al., J. Text. Inst. 19, T286 (1928); 
21, T511 (1930); 21, T561 (1930); 25, T1 (1934). 

8’ Robt. W. Webb, A. S. T. M. Standards on Textile 
Materials, Prepared by Committee D-13, pp. 321-339 
(October, 1940). 

® Sigurd Kohler, J. Text. Inst. 25, T141 (1934). 
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analogous to those occurring in a consideration 
of an actual yarn. 


THEORY 
Part I. Fibers All Alike 
Case A. Fibers Parallel 


Let the yarn under consideration be an 
extended bundle of fibers which are aligned 
parallel to the lengthwise axis of the bundle and 
in which the right end of any particular fiber is 
just as likely to lie at one place as another along 
the bundle axis. Let BB’ in Fig. 1 be the cross 





Fic. 1. Sketch of yarn composed of parallel fibers. The 
x coordinate is taken along the axis of the yarn. BB’ 
indicates the cross section at which the yarn strength is 
desired. 
section at which the yarn strength is desired. 
Let the strength of the yarn at this point be 
contributed to by each fiber of the cross section 
in an amount equal to that necessary to cause 
the particular fiber to break or slip, as may be 
the case. In a very long bundle of length X 
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containing N fibers of length Z there will be 


Ndx/X fibers with their right ends in the 
element of yarn length dx. Also n=NL/X will 
be the number of fibers in a cross section. 
Therefore the number of fibers with their right 
ends in dx will be ndx/L. 

Define J as the breaking strength of a fiber. 
Then if f, the frictional force per unit length of 
fiber, is such that f<2J/L the cumulated 
frictional force over half the length of a fiber 
will be less than the fiber strength and no fibers 
will break when tension is applied to the yarn. 
In this case the strength contributed to the yarn 
by those fibers which extend to the right of BB’ 
by not more than L/2 will be given by 


Li? fxn fnL 
f —de= 
0 L 8 


Similarly, the strength contributed by fibers of 
the cross section BB’ which extend more than 
L/2 to the right will be given by 


& n fnL 
f f(L—x)—dx =—. 
L/2 L 8 


The strength of the yarn when no fibers may 
break then becomes 


Fay=fnL/4, f<2J/L. (1) 


Now let fxo=J, where xo L/2, so that a 
fiber extending a distance of at least xo on each 
side of BB’ will break before slipping. The 
strength contributed by those fibers which do 
slip before breaking will then be 

—s bs” n 

F=f f—«xdx+ f—(L—x)dx = f—x,’, 

o £ t-x~ L L 
whereas the strength contributed by those which 
break before slipping will be 


L—zo n n 
r= f J—dx = J—(L—2x»). 


Therefore the total strength in this case where 
fibers may break will be given by 
Fao=F,+F,=nJ(i-—J/fL), f22J/L. (2) 


Now the frictional force, f, along a fiber will 
depend upon the forces which press the fibers 
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together and upon the coefficient of friction. 
Let p be the average pressure exerted upon a 
fiber by the adjacent fibers. Then f=ypy, where 
y is the perimeter of the fiber cross section taken 
normal to the fiber axis and uy is the coefficient 
of friction. But y=aéS where a@ is the fiber 
cross-sectional area, 6 the fiber density and S 
the surface per unit mass of the fiber.'® (In the 
latter part of this paper where the fibers are not 











a b Cc 


Fic. 2. Geometrical relations for a yarn containing twist. 


considered to be all alike, S=total surface of 
all the fibers/total mass of all the fibers.) The 
expression for f then becomes 


f=pabSp. (3) 


Upon combining Eq. (3) with Eq. (1) and with 
Eq. (2) the yarn strengths may be written, 
Fa,=npabSpL/4, Faz=nJ(1—J/pLadbSp). These 
expressions may be presented in a more useful 
form by introducing the quantities J = mass per 
unit length of yarn, A\=strength of fiber in 
strength units per mass per unit length of fiber. 
Thus for a yarn in which fibers which may not 
break, 


Fa,;= MpSpL/4= MBp/4, p<2d/B; (4) 

and for a yarn in which fibers which may break, 

Faz= MX(1—2X/Bp), p22Xd/B8. (5) 

The symbol B=LyzS is used for convenience in 
later portions of the theory. 


Case B. The Bundle of Fibers Contains Twist 


Let the bundle of fibers now contain twist so 
that the fibers lie as in Fig. 2a. Q=rg=tan 86, 
since a unit length of yarn is taken. g=2zr7, 





1 R. R. Sullivan and K. L. Hertel, Text. Research 11, 
30-38 (1940). 
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where 7 is the number of turns of twist per unit 
length. Let T in Fig. 2b be the average force or 
tension per unit area where the force is taken 
along the fiber axis and the area is taken normal 
to the fiber axis which makes an angle @ with the 
x direction. Let z be the area width while dr is 
the area height normal to the yarn axis. The 
tension over the area 2dr is Tzdr. It will have a 
component, 72dr cos 6, in the x direction and a 
component, 7 2dr sin 6, normal to the x direction. 
As one proceeds along a fiber the change of 
direction of the latter component will tend to 
compress the yarn. If df in Fig. 2c be the 
small force toward the yarn center due to the 
change of direction of this latter component 
over a small angle dg, then df= Tadr sin @rdg/r, 
and the small pressure increase produced at this 
position will be 
df T sin? 6dr Tgy?*rdr 
—dp=—___—___= ——_ = ——. (6) 
rdgz/sin 6 r 


The negative sign is used with dp since the 
increase is toward rather than away from the 
yarn center. 

If t is the tension per unit cross-sectional area 
where ¢ is in the x direction, and the area is 
normal to the x direction, then 


Tzdr cos 6 T 
t{=———_—_——_-=T cos? @=— 


—., (7) 
drz/cos 0 g*r?+1 


The total strength of the twisted yarn is given by 


ro rm=ro 1 dp 
F=f ianrdr= J —— —2rrdr 
- gr dr 


0 r=0 


P=P0 De 2r 
=,  =dp==[p*- pul. 
2 y? 


p=p* 


1 
F=—[p*—po], (8) 


2nr? 


where p* is the pressure at the yarn center 
(r=0) and fp is the pressure at r=ro=radius 
of yarn. 

In order to obtain the value of p*, Eq. (6) 
must be integrated. Since different fibers at the 
same distance r from the yarn axis will, in a 
given yarn cross section, have different tensions 
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in accordance with their extension from the 
cross section, they will have different values of 
dp associated with them. The average value 
which dp would have for the elementary ring of 
radius 7, if the ring were a representative sample 
of fibers, will now be used for all parts of the 
ring. This involves using the average value of T 
for a given ring. Average values of T will be 
obtained from Eq. (4) and Eq. (5) for cases 
where fibers may not break and where fibers 
may break, respectively. Thus, 


Fa, MBp 8 Bp 
T= = —_ = ——* 
A A4 4 


Fao Md d r 
T.= = -(1-—)=0n(1-—), (10) 
A A Bp Bp 


where A is cross-sectional area normal to fiber 
axis for the non-twisted yarn and 6’ is the density 
of the yarn. 

Combining Eq. (9) with Eq. (6), 








(9) 








dp 8B g*rdr 


p 4 g*41 





Upon integrating and imposing the condition 
that p= pp» at r=ro, the following expression for 
p as a function of r when p=2)/8 is obtained: 


gro? +1 5’B/8 
p=p(———) . (11) 
g’*r?+1 

Combining Eq. (10) with Eq. (6), 
pdp rg? ~—s rdr 


(8p—r) B (y%+1) 





Upon integrating and imposing the condition 
that p=2d/8 at r=r,=radius at which fibers 
just begin to break, the following relation 
between p and r when p> 24/8 is obtained: 


erty"? Bp 
( —-) = (~-1 ons, (12) 
grt d 


where e is the base of natural logarithms. 7; is 
obtainable from Eq. (11) by setting p=2d/8. 


1 BPo 8/6’B \' 
‘= ( *) (g*ro?+1)—1}. (13) 
gl 2d 
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In case no fibers may break, Eq. (11) gives 


dar? M 6’A/8 
p*= pol +1) » p*t=2xr/B. (14) 
6’ 


In obtaining Eq. use was made of the 


relation 


(14) 


gro? =497r°M r. 


(15) 


In case fibers may break, Eqs. (12), (13), and 
(15) combined give 


BPo ‘7497? M 6’A/2 
GO) (+) 
2d 5’ 
Bp* 
-( -1)ern-s, (16) 
nN 


Equation (16) may not be solved explicitly for 
p* but may be solved for r. 


(= ) 5's —) 6's 
= —1 
LIN (— 


, nN 
«K e2p*/b’n 4/6’8 1 


] , p*>2d/B. (17) 
4nM 





Let 7; be the number of turns per unit length at 
which fibers may just begin to break at the yarn 
center. 7; may be obtained from Eq. (17) by 
setting p* = 2d/p. 


| Qr \ 8/88 s’ 7 

Ti=141- —1} ? 

| (—) ks] 

In order to obtain values of yarn strength as 

a function of +, p* is obtained directly from 

Eq. (14) for r<7r; and used with Eq. (8). For 

values of r>7r, values of p* are postulated and 

corresponding values of + obtained from Eq. 

(17). These are then used with Eq. (8) to obtain 
yarn strength. 


(18) 


Part II. Fibers Not All Alike 


Where the fibers in the yarn are not all alike 
the analysis becomes less simple. A rigorous 
treatment demands a detailed knowledge of the 
joint distributions of the variable quantities. 
A further source of difficulty is the fact that 
where fibers are not all alike some fibers will be 
subject to breakage at a given pressure while 
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other fibers may not be subject to breakage at 
this pressure. Thus the transition from a rela- 
tionship of the form of Eq. (4) to one of the 
form of Eq. (5) is not abrupt but takes place 
over a relatively wide and important range of 
pressures. Another source of difficulty arises 
from the fact that in a yarn with a finite number 
of fibers in a cross section the mean value of any 
particular fiber will be somewhat 
different at different points along the yarn. In 
spite of these difficulties it is felt that the 
following treatment by approximate methods is 
sufficiently complete to indicate the manner in 
which the variability of the fiber properties 
studied affects the strength of the yarn at a 
representative cross section. 

In this study the fibers considered will all 
have the same density, 6, and the same coefficient 
of friction, uw. The properties of fiber length, 
fiber perimeter, fiber shape, and strength per 
mass per unit length will vary from fiber to fiber 
but will be uncorrelated. These properties will 
be represented by L, y, y and AX, respectively. 
y, the shape factor, is by definition the ratio of 
the fiber cross-sectional area to the square of 
the fiber perimeter. The yarn strength at a 
representative cross section of the yarn will be 
obtained by finding the mathematical expecta- 
tion of the strength at any cross section. 


property 


Case A. Fibers Parallel 
When a yarn with n fibers in a cross section is 
composed of parallel fibers and the pressure is 
such that no fibers may break, Eq. (4) indicates 
the mathematical expectation of the strength 
to be 


nu . 
E(Fa;) =—pE(wLs), (19) 
4 


where the symbol E is used to indicate mathe- 
matical expectation, w is the mass per unit 
length of a fiber, and s is the value of S for a 
particular fiber. Since the mathematical expecta- 
tion of a product of stochastic variables is the 
product of the expectations when the variables 
are independent, Eq. (19) may be written 


M Sup 


app 
E( Fa,) = oe ‘. ‘E(L). (20) 
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It is noted that in Eq. (20), nE(ws) = nE(surface 
of fiber per cm length) = (total fiber surface in 
varn per cm length)=(mass of yarn per cm 
length) X (fiber surface per gram of yarn) = MS. 
The symbol, E(L), refers not to the value for 
the whole yarn but rather to the biased sample 
caught in the cross section. Since measurements 
of mean fiber length are usually made with 
respect to the whole mass of fiber from which 
the yarn is spun," it is desirable to convert E(L) 
to E*(L) where E* refers to the expectation of the 
length for the population of fibers composing the 
varn. 

Let the subscript 7 refer to a particular length 
of fiber. Then E(L)=>on,Li/>on;, where > 
represents the summation over all 7. By referring 
to case A of Part I it may be observed that 
E*(L)=>-NiLi/>-N; and that n;=N;L;/X. 
Thus E(L)=>0N;,L2/>0NiL; and E(L)- E*(L) 
=(VN.L2/UNiLi)(TN:iL:i/CNi)=LUNL?2, 
> N;=E*(L*). From this, E(L) = E*(L*) E*(L). 
But the variance of L for the whole population 
is, by definition the square of the standard 
deviation, o, and is given by o®?= E*(L—E*(L))? 
= k*(L*) — | E*(L)}*. Thus 

E(L) = E*(L)-{1+[o/E*(L) F} 

= F*(L)-{1+C,7}, 


where C, is written for the coefficient of variation 
of L, which is by definition the ratio of o to 
E*(L). 

In view of these findings, Eq. (20) for the yarn 
composed of fibers which may not break becomes 


MypSE*(L)-(1+Cxz? MB’ p 
E( Fa,) = - - p= 9 (21 
4 4 
where 
B’ = pSE*(L)-(14+C_,?) 
It should be remarked at this point that 


Hertel, by a geometrical treatment of the curve 
plotted by the fibrograph, has already indicated 
in different terminology that E*(L)-(1+C,?) 
rather than E£*(L) should be used for Ptah 
varn strength for those fibers which may not 
break." 

Returning now to the case where fibers of all 
types may break, Eq. (5) leads to the following 


"kK. L. Hertel, Text. Research 10, 510-525 (1940). 
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value for the expectation of the strength: 


9 


1 wr* 
up sL 


Use will now be made of the relations, s = (yéy)~'; 
w= (v5)y*. Thus, 


1 
E( Fas) = n| E(w) -E(A) _ —E(\?*) 


up , 
noevy-2()| 
L 
| 6° E(\?) 
= nE(w)-E(A)\; 1 -——-— 
| upE(aA) 


It follows from the definition of S that 


S=nE(ws) /nE(w) = E(y)/{ E(y)-E(y*) - 6}. 
Also E(w)=6E(y)-E(y*). The expression for 
E( Faz) then may be written 

7 E(x’) 
E( Faz) = M- E()) 
r ~ upS-E(n 
EM): -E(y*®)- E(y) 1 
= 5 ? 4) +. (23) 
{E(y)}2-{E(y)}2 ALI 


It will now be necessary to obtain E(1/L) in 
terms of E*(L). From case A of Part I, n; 
=N;L;/X and n;/L;=N;:X where subscript i 
refers to a particular length group. Then if >> 
represents the summation over all i, E(1/L) 
=> (n:/Li)/Doni=>Ni/TNLi=1/E*(L). 

At this point, for simplicity, the fiber perime- 
ters will be considered normally distributed. 
This type of distribution of perimeters, in fact, 
is very nearly that actually found for cotton 
When the normal distribution of y is 
assumed, it is not difficult to show that 


fibers. 


E(y*)- E(y) 
| E(y?) }? 


1+3C,? 
~ (1+4C,2)? 


where C, is the coefficient of variation of y. 
There also exist the identities, E(d*)/E(A\) 
= E(A)-(14+G2); E(y’)/(E(y) }?=(1+C,’) where 


C, and C, are the coefficients of variation of \ 
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and y, respectively. In view of these findings 
Eq. (23) for the expectation of yarn strength 
where all types of fiber may break becomes 














| E(x) -(1+C)*)-(1+C,*) (1 +3Cy’) | 
x}1- bettie Ativan 
| E*(L) -pSp(i+C,?)? 
= ME(\)(1—E(A)/6’’p), (24) 
where 
” E*(L)pS(i+C,*)? 
(1+ C,2)(1+C,2)(14+3C,2)_ 
J soo} tee 
Ss 
= 
i] 
o 
F 200; 
H 
> Ke 
100 L i i 1 
65 70 75 80 


LOG, PRESSURE (LOG, ¢/cu*) 


Fic. 3. Linear relation between yarn strength and log, p 
for the intermediate range of pressures, pi=p=p2. The 
values used are for the yarn of Fig. 4 which is composed 
of fibers whose properties vary from fiber to fiber. 


Equations (21) and (24) provide means for 
studying the low strength and high strength 
regions of a yarn strength vs. pressure curve. 
The intermediate region of the curve may be 
studied readily enough, when only one fiber 
property is variable, by establishing the distribu- 
tion function of that property, and then from 
such distribution function determining the 
number of fibers which may break, and the 
number of fibers which may not break at any 
specified pressure. The proper mean values are 
then calculated for these groups and applied to 
Eqs. (21) and (24), the total strength being the 
sum of the contributions from the two sources. 
This type of treatment has been applied to 
each of the properties separately where, for 
definiteness, a normal distribution of the prop- 
erty was postulated in each case. As might be 
supposed from Eqs. (21) and (24), when the 
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variabilities are not too great the results are not 
particularly sensitive to the type of distribution 
used. This type of treatment becomes very 
unwieldy, however, when all the fiber properties 
are allowed to vary simultaneously. 

The study of the intermediate portion of the 
yarn strength vs. pressure curve is very much 
simplified by the fact, discovered by trial, that 
this intermediate range is well represented, for 
the range of variabilities ordinarily encountered, 
by an expression of the form 


F;= Ma log (bp),  pi=p=pr. (25) 
The values of p: and p2 and of the parameters a 
and b of Eq. (25) were determined by imposing 
two sets of conditions: (1), Eqs. (25) and (21) 
must, for p=p:, give the same values of yarn 
strength and also the same values of rate of 
change of yarn strength with pressure ; (2), Eqs. 
(25) and (24) must, for p=po, give the same 
values of yarn strength and also the same values 
of rate of change of yarn strength with pressure. 
These conditions require that 


4a [E(A) }? 8° 
pi=—; pPo=——; b=—e; (26) 
3’ B’'a 4a 
and that 
E(A) /8'\3 
2 B” 


where e is the base of natural logarithms. That 
the value of a must be obtained by trial from 
Eq. (27) is unfortunate but does not preclude 
its use. Equation (25), where a and } were 
determined as indicated, was found by trial to 
accurately the intermediate range 
when a single fiber property was variable, or 


represent 


when all the fiber properties studied were varied 
simultaneously. Normal distributions of the fiber 
properties were postulated. Figure 3 shows a 
linear relationship between F; and log p for the 
intermediate range when all the properties were 
varied simultaneously, the values of the fiber 
properties being the same as for the yarn C of 
Fig. 4. 

Summing up then, the expectation of the yarn 
strength at any cross section of a yarn composed 
of parallel fibers with properties varying from 
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fiber to fiber is given by 


F,=E(Fa:) = MB’'p/4, p<4a/s’; (28) 
4a {E(A) }? 
F;= Ma log (6p), —=p=———; (29) 
, ap”’ 
E(a) 
‘, = E( Faz) => MEQ)(1 -—), 
B''p 
| E(A)}? 
$2——,, (%) 
ap’ 


where the symbols used are as previously defined. 

Figure 4 shows the form of a curve given by 
Eqs. (28), (29), (30), and compares it with a 
curve for the same yarn obtained by combining 
the results of separately calculating the effect of 
variability of each fiber property when each 
property was normally distributed. A third curve 
is given for a yarn whose fibers are all alike and 
whose properties have values equal to the 
mathematical expectations of those for the 
previous yarn. 

The hypothetical 20’s yarns of Fig. 4 were 
composed of fibers with E*(L)=3.0 cm, E(A) 
=1.52X10° grams/grams/cem, M=2.95X10~ 
gram/cm, S=2.84X10* cm?/gram, p=0.299. 
For the yarn where the fibers were not all alike 
the fiber properties had coefficients of variability 
as follows: C,=33.3 percent, C,=25 percent, 
C,=20 percent, C,=40 percent. All these values 
except the strength, strength variability, and 
coefficient of friction were obtained by actual 
measurement of cotton fibers. The values of A 
and yw were calculated from data published by 
C. Nanjundayya and N. Ahmad,” and K. R. Sen 
and N. Ahmad." The value of C, was postulated. 


Case B. The Bundle of Fiber Contains Twist 


When the twisted yarn is composed of fibers 
with different properties there are three regions 
corresponding to Eqs. (28), (29), (30) to be 
considered rather than the two regions corre- 
sponding to Eqs. (4) and (5), when all fibers 
were alike. Aside from this difference, the 
treatment is similar. Thus, when the integrations 


T75-109 (1938). 


3K. R. Sen and N. Ahmad, J. Text. Inst. 29, T258-279 
(1938). 


VOLUME 13, MARCH, 1942 


are performed it is found that 


F= 





Lp*— po] 


2nr? 


as before, and that 


T°42rM 6°6'/8 
p= ro( +1) reas 





(31) 
Aer 2 |* 5’ 
[© bec! lal thew 
pi=p*=p2; (32) 
SL ey 8” (Bp t— EQ) eee" 
“5 | 8"p:—E(a)] 
POP ey 
Xexp (p caw ~ | —, 
p2<p*; (33) 


where the symbols have the same meaning as 
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Fic. 4. Values of yarn strength for 20’s yarns are given 
for different pressures. The yarns are composed of parallel 
fibers. Points A are for all fibers alike. Points B are for 
fiber properties variable as indicated in text, and are 
calculated from Eqs. (28), (29), (30). Points C are for 
fiber properties variable and are obtained by combining 
the results of separately calculating the effect of variability 
of each fiber property where each property was normally 
distributed. 
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before, and 


|p p g(bp) 
I -{ 
PI p log (bp) 


=log (log (bp)) +log (bp) 


(log (bp))? (log (op))’ » 
+ ss le | (34) 


2-2! 3-3! 


Pl 


Figures 5 and 6 show curves for yarn strength 
vs. twist as calculated from Eqs. (8), (31), (32), 
(33) for variable fiber properties and as calcu- 
lated from Eqs. (8), (14), (17) for fiber properties 
not variable. Except for the twist, the yarns of 
Figs. 5 and 6 are the same as those of Fig. 4. 
The external pressure postulated to hold the 
fibers together was po=0.6 g/cm?. This value 
was chosen in order that the initial strength of 
the non-twisted yarn might be equivalent to 
that found experimentally on a piece of un- 
twisted sliver in the writer’s laboratory.* It is 
observed that over the upper range of twists 
the effect of the variable fiber properties on the 
strength-twist relationship is very small for the 
representative yarn cross section here considered. 
In an actual yarn, however, variable fiber 
properties would contribute to the ‘‘unevenness”’ 
of the yarn, since, statistically, not all yarn 
cross sections could be representative. It is 
encouraging that the curve of Fig. 5 has the 
characteristic form observed in actual practice." 


Part III. Maximum Strength and Optimum Twist 


In view of the smallness of the effect produced 
by fiber-property variability for a representative 
cross section when high twists are used, the 
discussion of maximum strength and optimum 
twist will be confined to the simpler case where 
all fibers are alike. A precise method is to plot 
a strength vs. twist curve for each case and take 
values of maximum strength and optimum twist 
_directly from the curves. A more practical 


* This was done by replacing Fa, of Eq. (4) by the 
strength of the untwisted sliver and then solving for p, 
the value of 8=LyuS and of M, the sliver mass per unit 
length being known. The resulting value of p was used 
for po. It might be pointed out that the upper portion 
of the curve in Fig. 5 is not very sensitive to the value 
used for po. 

4 A. A. Mercier and Chas. W. Schoffstall, J. Research 
Nat. Bur. Stand. 1, 733-750 (1928). 
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method, however, is to replace Eq. (17) by the 
approximate relationship, 
| 2X ‘ee 5’ 
ta] -— e2P*/5's _ i. — (35) 
| spo! 4nM 


The omission of (8p*/A—1)*/°’8-e-*/°'8 from Eq. 
(17) is the neighborhood of 
400 


not serious in 
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Fic. 5. The parallel-fiber yarns of Fig. 4 have been 
twisted and the resulting yarn strengths plotted as func- 
tions of the twist. The constant external pressure, po, is 
0.6 g/cm?. 6’=unity. Points A are for fiber properties 
variable. Points B are for all fibers alike. 
optimum twist, since the expression is nearly 
unity for values of the parameters ordinarily 
used. Equation (35) gives values of 7 in the 
neighborhood of optimum + which differ by 
about 1 percent from those given by Eq. (17). 
Equation (35) has the attractive feature that it 
may be solved explicitly for p*. Thus, 


5X 4x Mr? 
pra? og eM 41) 
2 5’ 


Dr 8/6’B 
~loe(—) ). (36) 
Bpo 


An approximate value for the yarn strength is 
then obtained : 


5’ | 4x Mr* 
F ts —— loe( +1) 
| 5" 


4rr® 
2d \ 8/58 
—loe(—) ; (37) 
Bpo 
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The logarithms are to the base e. Equation (37) 
gives values of F for a given r in the neighbor- 
hood of optimum 7 which agree within about 
2 percent with those obtained by use of Eq. (17) 
and Eq. (8). This degree of approximation, while 
perhaps not so high as might be desired, is 
sufficiently good to allow a study of the manner 
in which optimum twist and maximum strength 
vary with fiber properties. 

Consequently, upon differentiating Eq. (37) 
and setting dF/dr=0 for a maximum value of 
F, the following relationship is obtained: 


Bpo\ 8/98 4n Mr? 
(2a 
2r 4nMr?+8’ 
42x Mr? tah 
x(-—+1) . (38) 


At this point it appears desirable to introduce 
the concept of twist multiplier. By definition, 
the. twist multiplier, K, is the number of turns 
per inch of yarn divided by the square root of 
the yarn number. For cotton yarns the yarn 
number, m, is by definition equal to the length 
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Fic. 6. The values of Fig. 5 have been plotted on semi- 
logarithmic paper in order that the lower portion of the 
curve may be given in more detail. Points A are for 
variable fiber properties. Points B are for all fibers alike. 
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of yarn in yards divided by 840 times the mass 
of the yarn in pounds. After suitable conversion 
of units, 


K =33.07rM!}, (39) 


where r=turns/cm, and M=grams of yarn per 
cm length. It is now worth noting that Bpyp>=4Fo, 
where Fy is the yarn strength per mass per unit 
length of the original non-twisted bundle of 
fibers. In the case here treated this strength 
was imparted to the bundle by means of an 
external pressure fo. In an actual yarn the small 
initial strength is due to other causes. Equation 
(38) may then be written, 


2F»o 8/B b'vK2, vK,? —_ 
(Yano) 
r vKo?+ 5’ 5’ 


where Ko is the optimum value of K; v is written 
for the constant, 42/(33)*; and where grams are 
used for mass and for strength, and cm for length. 

By applying Eqs. (38) and (39) to Eq. (37) 
the expression for maximum yarn strength 
becomes 


vK" —1 or vKo2 a 
Pus ® AM (1+ ) = (14 ' +) i 
5’ m 8’ 
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where m is the yarn number and ¢=0.00590 is a 


constant to take care of the units involved. 
Fuax Will be in grams. 

Equation (40) cannot be solved explicitly for 
Ko. However, the right-hand member may be 
plotted as a function of Ky for any particular 


value of 6’ and when the left-hand member is 
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Fic. 8. Calculated values of optimum twist multiplier 
are given for yarns composed of fibers with various 
properties. 


known the corresponding value of Ky may be 
read directly from the curve. Curves correspond- 
ing to three values of 6’ are given in Fig. 7. 
It is seen that to some extent Ky depends upon 
the value of 6’. The density of a hard-twisted 
cotton yarn is in the neighborhood of unity 
although it varies somewhat with the degree of 
twist. In the remainder of the discussion 6’ will, 
for definiteness, be taken as unity for the given 
yarns. 

When this theory is used, the optimum twist 
multiplier for the desired value of (2Fo/A)**, is 
taken from the curve of Fig. 7. The maximum 
strength is then obtained by applying the value 
of Ky so found to Eq. (41). The manner in which 
optimum twist varies with fiber properties is 
shown by Fig. 8. The manner in which maximum 
strength varies with fiber properties is shown by 
Fig. 9. It is observed that over the range of 
fiber properties studied, the value of Ko is largely 
determined by the value of 8, whereas the value 
of Fimax is much more sensitive to \ than to 8. 


CONCLUSION 


It is hoped that the foregoing analysis will aid 
in clarifying certain aspects of the complicated 
relationship existing between the properties of a 
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yarn and the properties of the fibers composing 
the yarn. For the sake of simplicity, and since 
yarn strength was the yarn property under 
consideration, several fiber properties were 
omitted from the analysis. How serious these 
omissions may be, remains to be seen when 
sufficient data are available for testing the 
theory. It appears to the writer that a likely 
cause for disagreement of theory and experiment, 
particularly for cotton fibers, is the postulate in 
the theory that any particular fiber has uniform 
properties throughout its length. As already 
pointed out, an actual yarn will break at its 
weakest point, whereas the theory gives the 
strength at a representative point. Furthermore in 
the yarn here studied all parts of a given fiber 
lie at the same distance r from the yarn axis, 
whereas the slight tangling of fibers necessary to 
give an actual untwisted yarn its initial strength 
will have an undetermined effect on the fiber 
orientation, and hence on the yarn strength, at 
any degree of twist. Also, certain of the fiber 
properties, such as length and fineness in a cotton 
yarn, are most likely correlated. The theory, 
however, yields reasonable values of Ko and Fyax. 














10 
VALUE 
OF @ 
8 4400 
2200 
ra 1100 
o 
° 
x 
6E 
3 
< 
= 
34 
3 
u 
2 = 
re) l —_— | 
te) 1 2 3 4 


A (10 ®x Grams /Grams Cm) 


Fic. 9. Calculated values of maximum strength for 20's 
yarns are given for yarns composed of fibers with various 
properties. 
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(See Figs. 8 and 9.) It appears, therefore, that for 
high twists the fibers behave in an actual yarn 
much like those in the idealized picture. 

If one wishes to use this theory to predict the 
optimum twist multiplier, Ko, and the maximum 
yarn strength, Fax, which the present idealized 
picture would predict for a 20’s yarn spun from 
a given group of fibers the simplest procedure 
would be as follows: 

(1) Determine LZ in cm, the mean fiber length 
by number; determine S, the surface per unit 
mass in cm* per gram; determine \, the fiber 
strength in grams breaking load per gram per 
cm of fiber; determine yp, the coefficient of 
friction for fiber against fiber (use »=0.3 if a 
more precise value is not available) ; determine 
Fo, the strength of the original non-twisted 
sliver in grams breaking load per gram per cm 
length of sliver (if no more precise value is 
available use Fo=8po/4 where B=yuLS and 
po=0.6 g/cm?*). 

(II) Calculate B=yLS; calculate Fo/X. 

(III) From Fig. 8 find the value of Ko corre- 
sponding to 8 and Fo/X. 

(IV) From Fig. 9 find the value of Fryax 
corresponding to 8 and X. 


It appears desirable at this point to restate 
the simplifying assumptions upon which the 
above mathematical treatment is based. Any 
given fiber was considered to have uniform 
properties along its length. Fibers were supposed 
to be flexible but not extensible. The properties 
of length, fineness, and strength were allowed to 
vary from fiber to fiber, but these properties 
were uncorrelated. The yarn was assumed con- 
structed so that any given fiber had all its 
elements at the same distance from the yarn 
center. An elementary ring normal to and 
concentric with the yarn axis was supposed to 
contain a representative sample of the fibers in 
the yarn cross section. Only the mathematical 
expectation of the yarn strength at a cross 
section was calculated. 

It is hoped that further work will remove some 
of the above restrictions from the analysis. For 
example, the incorporation of elastic properties 
would be desirable." 

It is a pleasure to acknowledge the very 
helpful discussions with Professor K. L. Hertel 
during the course of the study. 


16H. P. Gurney, J. Text. Inst. 16, T269-T289 (1925). 





An Extension of Saint-Venant’s Principle, with Applications 


J. N. GoopreR 
Cornell University, Ithaca, New York 


(Received November 18, 1941) 


Simple energy considerations, which have previously been employed to provide a rational 
basis for the principle of Saint-Venant, are shown to lead to the conclusion that forces applied 
in the neighborhood of a rigidly fixed portion of an elastic solid can cause only local stress and 
strain. This principle indicates useful features of certain problems encountered in engineering, 
and is used here in the calculation of the stiffness of a thin layer of elastic material confined 
between rigid plates, under any kind of load applied to the plates. 


I 


CCORDING to Saint-Venant’s principle, 

the stress and strain produced in a solid by 
applying to a small part of it a system of forces 
statically equivalent to zero force and couple, 
are appreciable only in the immediate neighbor- 
hood of the loaded part. It follows that any 
redistribution of a load applied within such a 
small part, leaving the resultant load unchanged, 
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changes the stress only in the immediate 
neighborhood. 

It has been pointed out! that this principle is 
a simple consequence of conservation of energy. 
For, if the loaded part has linear dimensions of 
order a, and the force per unit area is of order p, 
there will be stress components (excluding sharp 

1]. N. Goodier, Phil. Mag. [7] 8, 607 (1937), and 24, 


325 (1937). See also O. Zanaboni, Rend. Accad. Lincei [6a] 
25, 117, 596 (1937). 
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concentrations due to discontinuities) all of order 
p in the loaded part and in its immediate 
neighborhood. The associated strain components 
of all types will then be of order p/E. E is 
Young’s modulus if the material obeys Hooke’s 
law, and otherwise it represents merely the order 
of magnitude of the slopes of the stress strain 
curve. 

In order to eliminate rigid-body movement, 
we may regard one infinitesimal surface element 
of the solid as remaining undisturbed in position 
and orientation. Let this element be taken within 
or close to the loaded part. Then the strain 
components of order p/E will entail displace- 
ments, within the loaded part, of order pa/E. 
The load, forming zero resultant, is not all in 
the same direction, so we do not underestimate 
the order of the work done during loading if we 
multiply pa*, the order of the force of a given 
type to be found in any subdivision of the 
loaded part, by pa/E. Thus we arrive at p*a*/E 
for the order of the work done, and so of the 
strain energy. 

The strain energy corresponding to stress of 
order p is of order p?/2E per unit volume. An 
amount of order p*a*/E is thus adequate only 
for a volume of order a*. This is evidently 
equivalent to saying that the stress and strain 
are appreciable only within distances of order a, 
which is Saint-Venant’s principle. 


II 


So far we have considered a distribution of 
force with zero resultant. Consider now a solid 
with either a finite or infinitesimal element 
rigidly fixed, and forces forming non-zero 
resultant and couple, applied to a small region 
adjacent to the fixed part or element—adjacent 
meaning within a distance of the same order as 
the dimensions (of order a) of the loaded part. 
Equilibrium is then maintained by forces induced 
at the attachment of the fixed part. It remains 
true that, with force per unit area of order p 
on the loaded part, there will be displacements 
of order pa/E and the work done will be of order 
pa*’Xpa/E or p*a*/E. Again this can provide 
strain energy, of appreciable value, only within 
a volume of order a*. We may therefore make the 
following addition to Saint-Venant’s principle. 

“If forces are applied to a body over a small 
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part adjacent to a rigidly fixed part, the stress 
and strain induced are appreciable only within 
distances comparable with the dimensions of the 
loaded part, whether the applied forces form a 
non-zero resultant or not.”’ 





Fic. 1. 


An exception in the case of incompressible 
material is made later. 

If the rigidly fixed part is small, of order a, 
Saint-Venant’s original principle applies, since 





Fic. 2. 


the loads would be balanced by reactions 
induced at the fixed part. But where the latter 
is extensive it is not self-evident that the 
reaction is confined to the neighborhood of the 
load. 


Ill 


Such a principle of course gives formal 
expression to what in many instances might be 
regarded as intuitively evident. It might be 
expected, for instance, that if the sides of a 
strip of elastic material are rigidly fixed, and 
forces are applied at the two ends (Fig. 1), only 
localized end stress will be induced. However 
when only one side of the strip (Fig. 2) is fixed, 
the answer is not so evident. The fixed side 
suggests the case of Fig. 1 but the free side 
suggests the strip under simple compression. 
The principle of the preceding paragraph of 
course indicates localized end stress only. This 
question arises in the idealized problem of a 
wall attached to a massive foundation, with 
stress due to thermal contraction of the wall, 
the foundation being treated as rigid. If the 
wall is free, a (uniform) fall of temperature T 
causes contraction a7 per unit length, for a 
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coefficient of expansion a. Tensile stress EaT 
(E is Young’s modulus) applied to the ends, 
annuls this contraction, and the rigid foundation 
may then be supposed attached. Removal of 
the tensile forces on the ends is equivalent to 
superposition of the system of Fig. 2. Thus, 
since the stress of the latter system is only 
appreciable at the ends, to an extent governed 
by the height of the wall, the wall will be under 
thermal tension EaT except near the ends, 
where there will be shear and normal reaction 
from the foundation. If the wall is sufficiently 
long, so that the end zones do not interfere, 
these reaction stresses will be independent of 
the length. 


IV. THE ELASTIC LAYER CONFINED BETWEEN 
RIGID PLATES 


The layer is supposed attached to the plates 
(Fig. 3). There are four fundamental types of 
loading: 

(i) The plates are pressed together, remaining 
parallel. 

(ii) The upper plate is sheared over, remaining 
parallel to the base plate. 

(iii) The upper plate is twisted about some 
vertical axis by a pure couple. 

(iv) The upper plate is rotated about some 
axis in its plane by a “‘bending”’ couple. 

In (i) the stress is usually taken as 


o:=—p, o:=o,=—vp/(1—»), 


rog Vas™ Tyz=0 


(1) 


(where v is Poisson’s ratio) corresponding to a 
force pA applied at the centroid of the plate, 
A being the area. This system satisfies the 
conditions of equilibrium and compatibility. It 
yields zero for the lateral strain components 
€z, €y, and the uniform value 
(1+v)(1—2v) p 
~~, (2) 
1—yp E 


for the strain e., thus satisfying the condition of 
attachment to the plates. But the condition 
that the edge is free from stress is violated, for 
there is the normal pressure vp/(1—yv) round the 
edge. 

To complete the solution we have to find the 
effects of applying a tensile normal edge stress 
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of this amount, when the top and base plates 
are completely fixed, or alternatively when the 
base plate is fixed and the other free to move. 
In either case the stress and strain induced are 
local to the edge. Since no appreciable strain 
will be thus induced in the interior, no appreci- 
able displacement of the upper plate occurs in 
the case where it is free to move. The approach 
of the plates is therefore determined by (2) very 
closely if the layer is thin. 

The tensile stress yp/(1—yv) applied round the 
edge will set up stress in the edge zone, involving 
shearing stress between the plates and the layer, 
but not in the interior. The layer is thus re- 
strained from lateral expansion by a circum- 
ferential ring of shearing force. It is immaterial 
whether there is attachment to the plate or not 
within this edge zone. J 

An exception must be made when the material 
is incompressible. In such material a hydrostatic 
state of stress can exist with zero strain energy, 
and the argument employed in paragraphs II and 
III, and the principle of Saint-Venant and its 
extension, do not apply. It will be evident, 
however, that any non-hydrostatic part of the 
stress must be localized. The general problem of 
the elastic layer, compressible or incompressible, 
attached to rigid surfaces, has been extensively 
discussed by Synge* and Hay.* 






LLL CLL LLL LL 


f CELL LL ALLL 


Fic. 3. 


In (ii) we may first regard every element of the 
plate as in a state of pure shear (Fig. 4, upper 
plate removed). Then the condition of rigid-body 
displacement of the upper surface is satisfied, 
and the amount of displacement corresponding 
to a given resultant force through the centroid 
of the area can be found. The edge is, however, 
subject to complementary shears. The removal 
of these, the upper plate being unconstrained, 
then induces only local disturbance at the edge. 


3]. L. Synge, Phil. Trans. Roy. Soc. A231, 435 (1933); 


Phil. Mag. [7] 15, 969 (1933); Trans. Roy. Soc. Canada 
[3] 31, 57 (1937). 
4G. E. Hay, Can. J. Research 17, 106, 123 (1939). 
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Since the interior remains practically undis- 
turbed, the force-displacement relation for the 
upper plate as obtained for the pure shear 
stands as a close approximation when the layer 
is thin. 

In (iii), consider first a circular disk form. It 
is clear that all conditions are satisfied if we 
regard the layer as a slice of a round bar in 
simple torsion. If the layer is not circular the 
torsional solution of a bar is not applicable 
since it involves warping of cross sections, or, in 
the layer, of the top and base surfaces. However 
a possible state of the layer is found by regarding 
it as a part of a larger circular layer. There will 
be shearing stress proportional to the radius, 
say kr, directed along circles concentric with the 
axis of rotation, from which r is measured. 
These shearing stresses intersect the non-circular 
boundary, and so on this boundary there are 
complementary shearing stresses. Again removal 
of the Jatter merely disturbs an edge zone. In 
the interior therefore we have the shearing 
stress kr. Because the layer is thin, this prevails 
except for a small edge zone. It forms a resultant 
force with components {kxdA, /{kydA, and 
since our loading is to be a twisting couple only, 
the origin must be taken at the centroid of the 
area, and therefore the axis of rotation is the 
axis normal to the layer through its centroid. 
The shear stresses form a moment about this 
axisk fr'dA or kI, where I, is the polar moment 
of inertia of the area about the centroidal axis. 
The value of k is thus M,/I, where M, is the 


Fic. 4 


applied couple. The shearing stress is M,r/J>. 
The rotation of the upper plate corresponds to 
that of a circular slice in which the torsional 
stress is kr, that is, M,h/GI,, h being the thick- 
ness and G the shear modulus. 
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In (iv) it is natural to consider displacement 
components 
u=0, 


v=0, w=6xze/h 


corresponding to rotation @ of the upper plate 
about a line in it parallel to and directly above 
the y axis (Fig. 5), and to fixture of the lower 


Pe 8 


a a 
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plate. On calculating the stress implied, however, 
it is found that it does not satisfy the differential 
equations of equilibrium, and it appears that 
the simplest way to satisfy them is to assume 


u=azt+bz*, v=0, 


with b= —(@/2h)1 (1—2v). The constant a is 
determined as 2v@/1—2y by the condition that 
the stress r,. should produce zero resultant force 
on the upper plate. The initial system assumed 
for this problem is thus 


6 12 
u= —a(»-—=), v=0, w=0Ox2/h 
1—2p 2h 


and the corresponding stress components are 


w= 6xz2/h, 


Ev Ps 
~ YL) (1 — 20) ‘ 
E(1—v) Pg 





O:>= re 


(+9) (12s) 7 


o( --), Tey =T,:=0. 


The action between the layer and the upper 
plate is that of a, only. If it is to be a ‘‘bending 
moment,” the z axis must pass through the 
centroids of the faces of the layer, in order that 
the resultant force should vanish. 
done, the moment formed is 





oS 


This being 


1-—yp 0 
M,=———————-EI, 
(1+v)(1—2v) h 
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about the y axis, and 


1-—yvp 6 
= —_______—_FE],,- 
(1+y7)(1—2p) h 
about the x axis, J, and I,, being §x*dA and 
JSxydA evaluated for the face of the layer. 
These agree with the corresponding values for 
a slice of a beam except for the factor 
(1—v)/(1+yv)(1—2p). 
In this solution the edge of the layer is subject 
to certain forces corresponding to the edge values 
of o:, dy, Tzz. The removal of these, with the base 


M, 


completely fixed, causes only changes localized 
round the edge. In a thin layer the solution 
prevails except in this edge zone, and the values 
obtained for the moments will still be a close 
approximation. 

It may be observed that the movement of the 
upper plate under the action of a ‘“‘bending”’ 
moment is a rotation of amount @ about an axis, 
agreeing with that of the couples applied to the 
plates if the axes are principal, and lying half- 
way between the plates; that is, the line through 
x=0, z=h/2 parallel to the y axis. 





The Phase of Arcback 


ALBERT W. HULL AND FRANK R. ELDER 
Research Laboratory, General Electric Company, Schenectady, New York 


(Received December 31, 1941) 


Thyratrons of FG-41 type were tested at over-voltage for 
arcback in a circuit which simulates operating conditions 
very closely. Memory oscillograms showed the time in the 
cycle at which the arcbacks occurred. At a mercury-control 
temperature of 45°C, the distribution was nearly random, 
with a maximum at the time of highest negative anode 
voltage, and depending on voltage in the same manner as 
the average arcback rate, taken over a long period. When 
the mercury-control temperature was raised to 64°C, all the 
arcbacks occurred at the beginning of negative voltage 
application. These results are interpreted on the theory 
recently advanced by Kingdon and Lawton, which attri- 


1. INTRODUCTION 


HE time in the cycle at which arcback 

occurs is of importance because of the light 
which it sheds on the cause of arcback. Early 
experiments, in which it was found that arcback 
occurred near the middle of the cycle, at the 
time of highest negative anode voltage, pointed 
to glow discharge as a probable cause.' Later, 
an observed correlation between arcback and 
residual ionization led to a thermionic theory,” * 
according to which loose particles on the anode 
were heated to incandescence by bombardment 
with residual ions. In agreement with this theory, 


1 Prince and Vogdes, Mercury Arc Rectifiers and Circuits 
(McGraw-Hill, 1927), p. 71. 

* Marti and Winograd, Mercury Arc Power Rectifier 
(McGraw-Hill, 1930), p. 41. 

3 J. V. Issendorff, E.T.Z. 50, 1085 (1929). 
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butes arcback to charging up of small non-conducting 
particles on the surface of the anode. At low vapor pressure, 
where the ions remaining from the discharge have time to 
diffuse to the electrodes before negative voltage is applied, 
the charging is produced by processes such as glow dis- 
charge and field emission, which depend on voltage; hence 
arcback is correlated with voltage. At higher pressure, 
diffusion is interfered with by molecular collision,.so that 
ions are present when negative voltage is applied. In this 
condition arcback occurs at the beginning of the negative 
cycle. 


arcback was observed to occur at the beginning 
of the cycle. More recently, Pakala and Batten‘ 
have reported experiments with commercial mer- 
cury arc rectifiers, in which about half of the 
arcbacks were found to occur at the beginning of 
the period of negative anode voltage, while the 
remainder were distributed at random through- 
out this period. They conclude that “there are 
several types of arcback; and that one type, 
which constituted 54 percent of the total, is 
occasioned by the momentary positive ion cur- 
rent.”’ 

A theory of arcback has been proposed recently 
by K. H. Kingdon® of this laboratory, which 


4W. E. Pakala and W. B. Batten, Trans. A.I.E.E. 59, 
345 (1940). 

5K. H. Kingdon and E. J. Lawton, Gen. Elec. Rev. 42, 
474-8 (1939). 
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Fic. 1. Circuit for testing the phase of arcback. 


appears capable of accounting for all arcbacks. 
The experiments reported in this paper were 
made under controlled experimental conditions, 
and will be interpreted in the light of this theory. 


2. METHOD OF MEASUREMENT 


Thyratrons of FG-41 type* were tested at 
over-voltage in an artificial circuit which simu- 
lated practical operating conditions as closely as 
possible. The circuit is shown in Fig. 1. A current 
of 90-amp. peak value from a sixty-cycle source 
was passed through the tube to be tested and 
two or more “hold-off” thyratrons, in series. It 
was limited to 10 amp. av. by first delaying the 
starting in each cycle, and then cutting off the 
current just after passing its peak value. The cut- 
ting off was accomplished as quickly as possible 
by starting a “‘commutating”’ thyratron, the re- 
duction from 90 amp. to zero being completed in 
4 electrical degrees. A definite time was then 
allowed, during which no current flowed, for ions 
to diffuse to the walls. This time could be varied 
from 6 to 30 electrical degrees. At the end of 
this time a negative voltage of definite value, and 
of a wave form similar to that encountered in 
rectifier circuits, was applied to the anode. This 


* The FG-41 tube has a rating of 75 amp. max., 12.5 
amp. av., 10,000 volts max. 
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—— “Memory” oscillograph deflection plates 









voltage was obtained from 
a charged condenser, CiC2, 
with the help of two thyra- 
trons. The first, connected 
to an intermediate tap on 
the condenser, delivered its 
charge to the anode through 
a low resistance (7000 
ohms), thus producing a 
steep voltage rise to about 
+ maximum. The second 
thyratron had a high series 
resistance (1.1 megohms) 
and produced a slow further 
rise of voltage to maximum 
value in about 45 electrical 
degrees. From this point the 
voltage decreased slowly, 
as the charge leaked off 
through a thyrite resistor. 
The wave forms of current 
and voltage are shown in Fig. 4. The voltage wave 
form can be seen also.in the oscillograms of Fig. 
6. It will be noted that the first part of it, which 
is important for arcback, is a fairly good simu- 
lation of the inverse voltage in a rectifier under 
operating conditions. 

Arcback was induced by over-voltage. In some 
of the tests the voltage was raised slowly until 
arcback occurred; in others it was set at a value 
at which arcback occurred about once in five 
minutes. The voltage at which arcback occurred 
varied between 38 and 43 kilovolts in the differ- 
ent tests. 

The phase of the arcbacks was determined 
from records of the anode-cathode voltage during 
the cycle in which the arcback occurred. These 
records were made by a memory oscillograph. 





Thyrite 
resistors 








3. THE MEMORY OSCILLOGRAPH 


The memory oscillograph is a high voltage 
cathode-ray oscillograph with willemite screen, 
operated in such a way that the photographic 
record covers the period from just before to just 
after the event which is of interest. To accom- 
plish this, the deflection plates of the oscillograph 
are connected permanently, through capacity 
dividers, to the anode and cathode of the tube 
under test, and the oscillograph allowed to write 
continuously. When arcback occurs, the shutter 
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of a camera is tripped by a thyratron-operated 
magnet, causing a photograph of the screen to 
be taken. The time lag of the magnet in opening 
the shutter is about } cycle; hence the ‘‘event”’ is 
over before the shutter opens. But since the per- 
sistence of the phosphorescent image is about 
1/25 sec., the record which the camera sees covers 
the period from approximately 2 cycles before the 
event to as long after as the camera remains open. 
This was usually 1/60 second. In the present 
tests a special relay (S, Fig. 2) stopped the 
writing at the instant the shutter opened, since 
events subsequent to arcback were of no interest. 
The photograph obtained, therefore, showed the 
arcback cycle and about two cycles prior thereto. 

The essential features of the memory oscillo- 
graph are shown in Fig. 2. It contained three 
standard RCA 912 cathode-ray tubes, with 5 
inch screens, rated 15,000 volts. They were 
operated at about 11,000 volts. In the present 
tests only one of the three tubes was used. 

Figures 3 and 5 show front and rear views of 
the oscillograph. The camera, with the relay 
which operates it, can be seen at the top in each 
photograph. The three cathode ray tubes can be 
seen in the center of Fig. 5. 


4. EXPERIMENTAL RESULTS 


Seventeen films of memory oscillograms were 
taken, with 16 records on a film, making a total 
of 272 arcback records. Typical oscillograms are 




















Test circuit 


Jt 



















































































E RCA 
P 912 Anode #2 <2 
i 20 sss 
= o’oO 
Anode *! 2: 8 
—_ med °ouy 
es a rd * — 
' LA. ea =z | | 
Cathode of Sweep 
ifier . | 
er fest circ t sss 


Additional bias to prevent 
cathode ray tubes from 
writing while camera 
shutter is open 


Fic. 2. Circuit for operating the memory oscillograph. 
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shown in Fig. 6. The ordinates are anode-to- 
cathode voltage, the abscissas linear time sweeps 
of 1-cycle length. The first record on each film is 
an undistorted wave, without arcback, showing 
the form of the impressed negative voltage. In 
the records which follow, the beginning of this 














Fic. 3. Front view of memory oscillograph. The camera 
and tripping relay may be seen at top. 


wave occurs at random times in the cycle, since 
the sweep was not synchronized. It can be recog- 
nized easily, beginning with a steep rise of volt- 
age, too fast to write, followed by a slower, well 
recorded rise. From this peak value the voltage 
decreases slowly, until arcback occurs; this is 
evidenced by an abrupt fall, followed by a partial 
rise. The rise is due to recharging of the line, from 
the still partially charged condensers C,C, (Fig. 
1), after the arcback is extinguished. The amount 
of recharging depends on the phase of the arc- 
back. In the records of film 8 (lower film, Fig. 6) 
in which arcback is seen to have occurred late and 
random in phase, the recharging is only +5 to } 
of the full voltage, and appears like a rise in the 
base line. In film 13 (upper film, Fig. 6) on the 
other hand, where the arcbacks occurred at 
the beginning of the cycle, the recharging is 
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Fic. 5. Rear view of upper half of memory oscillograph, 
open, showing cathode-ray tubes. 


almost complete; it appears as a doubling or 
broadening of the slow discharge trace ; the upper 
portion of the double line being the memory of 
the previous cycle. 

The results of the tests are shown in Fig. 4. 
The lower portion of this figure gives the form 
and sequence of current and voltage used in the 
tests. The sequence is: 

1. Ninety electrical degrees of anode current, 
being a portion of a sine wave of 90 amperes 
maximum value. 
commutation time, during 
which the anode current dies to zero. 

3. Six degrees (10° from beginning of commu- 
tation) of time delay, before applying negative 


2. Four-degree 


voltage, to allow left-over ions to diffuse to walls. 
This time was increased to 18 degrees and 35 
degrees, respectively, after beginning of com- 
mutation, in the case of films 12 to 17. 

4. Negative voltage applied to anode, rising 
steeply at start, then more slowly to a maximum 
in 45 degrees, decreasing slowly thereafter for the 
remainder of the cycle. 

The times in the cycle at which the 272 arc- 
backs occurred are shown in the upper part of 
Fig. 4. Each dot represents an arcback, plotted 
against time in the cycle, in electrical degrees. 
It will be noted that the arcbacks recorded on 
films 1 to 8 are distributed nearly at random. 
On the other hand, in each of the films 9 to 12, 
all sixteen arcbacks occurred in the first 10 
degrees, which is indicated by the number 16 in 
the square bracket. Film 16 had 6 arcbacks at 
the beginning, the rest distributed ; in film 17, all 
were distributed. The interpretation of these 
results is discussed in Section 6. 


5. KINGDON’S THEORY OF ARCBACK 


K. H. Kingdon has proposed that the cause of 
arcback is the charging up by positive ions of 
small particles of insulating material on the sur- 
face of the anode, thereby causing local fields 
strong enough for runaway field emission. For 
example, a particle of 4X10~-*-cm diameter on a 
negative anode to which a total positive ion 
current of 8X10-" amp. per cm? flows, will 
receive 12 ions, charging it to 2.8 volts and pro- 
ducing a field of 0.710 volts per cm at the 
metal surface just below it. Kingdon and Lawton® 
have given data which show that this field can 
cause arc-formation. 





Film 8. 





Film 13. 


Fic. 6. Typical memory oscillograms of arcback, from films 8 and 13. A normal cycle, without arcback, 
is shown at the beginning of each film. 
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That fields of this magnitude can cause field 
emission is well known. The transition from field 
emission to arc under these conditions has been 
demonstrated by experiments of Giintherschulze 
and Fricke® and of Koller and Johnson.’? Both 
pairs of experimenters, the first with a glow dis- 
charge and the second with secondary electron 
emission charging in vacuum, observed that the 
emitting particles are quickly destroyed unless a 
high series resistance is employed, the discharge 
taking the form of a shower of sparks. The field 
emission in this case is a runaway effect, burning 
up the particle and thus providing the ‘‘cathode 
spot”’ of arc-formation. 

While Kingdon and Lawton considered only 
the charging up of the particles by residual 
ionization, it is evident that any other mecha- 
nism of charging, such as glow discharge, photo- 
electric effect, or field emission, would have the 
same effect. For example, in multi-anode recti- 
fiers there is evidence that arcback is facilitated 
by allowing the anode to ‘‘see’’ a discharge in 
another part of the tube; this points to photo- 
electric effect as the charging means. On the 
other hand, many tubes are observed to arc back 
when glow discharge appears in the anode stem, 
although the glow itself does not become an_are. 
































GRID " 
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ne CARBON PLATE 
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Fic. 7. Outline structure of experimental thyratron. The 
distance between grid and anode is 6 mm. 


6 Giintherschulze and Fricke, Zeits. f. Physik 86, 463, 
821 (1933). 

7L. R. Koller and R. P. Johnson, Phys. Rev. 52, 521 
(1937). 
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In this case the charging is caused by either posi- 

tive ions or ultraviolet light from the glow. 
While the charging up of small insulating 

particles as a cause of arcback was proved by 
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sop tiie 


Kilovolts 


Fic. 8. Average rate of arcback as a function of voltage. 


Kingdon and Lawton only for arcbacks which 
occur immediately after commutation, there is 
good reason to attribute all arcbacks to this 
mechanism, with variations only in the means of 
charging. We shall assume this point of view in 
interpreting our results. 


6. INTERPRETATION OF RESULTS 


In mercury vapor at 45°C, the mean free path 
of a mercury atom is 2.92 mm. This is approxi- 
mately equal to the maximum diffusion distance, 
3 mm, to be traversed by a mercury ion in the 
grid-anode space of the FG-41 tube, in order to 
reach a discharging surface (Fig. 7). Diffusion 
under these conditions is very rapid and is 
essentially complete in 463 microseconds. Films 
1 to 8 (Fig. 4) were taken under these conditions. 
It is seen that no arcbacks were observed in the 
first 15 degrees after application of voltage. We 
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interpret this to mean that the means for charg- 
ing insulating particles on the anode were absent ; 
the residual ions had diffused to the walls, and 
any charges on the particles resulting therefrom 
had leaked off. (Very little net charging of the 
particles is to be expected from diffusion under 
these conditions, i.e. in the absence of voltage 
between the electrodes, since the plasma con- 
sists of equal numbers of ions and electrons.) 

Arcbacks were observed on these films (1-8) 
only after about 25 degrees, when the negative 
voltage had risen nearly to its maximum. At this 
point some ionizing process occurs, purely volt- 
age-dependent, which is capable of charging the 
insulating particles. We believe it to be long-path 
glow discharge in the anode stem, since such 
glow was observed frequently in these tubes at 
the limit of rectifying voltage. On this view, the 
limiting voltage of a rectifier, operated under 
these conditions, is a function of the geometry 
of the anode stem. 

The correlation of arcback with voltage for 
films 1-8 is shown in Figs. 8 and 9. Figure 8 gives 
the average rate of arcback as a function of 
voltage, determined by a scale counter over a 
long period. The product of the rates in Fig. 8 by 
the voltage in Fig. 4, gives the probability of arc- 
back due to voltage, as a function of time in the 
cycle, for the data of Fig. 4. This product is 
plotted as the solid curve in Fig. 9. The circles 
and bars in Fig. 9 give the observed arcbacks 
from films 1-8, by ten-degree intervals. The 
correlation is seen to be as good as should be 
expected from the number of observations. 

Films 9 to 11 were taken with the same short 
time allowance (463 microseconds) for deioniza- 
tion, but a higher vapor pressure (33 microns). 
In this case every arcback occurred within the 
first 10 degrees, at a time when the voltage was 
only half its maximum value. These arcbacks are 
attributed to residual ionization. 

Even when the time allowed for deionization 
was increased to 815 microseconds (film 12) all 
the arcbacks occurred at the beginning. With a 
further increase of deionizing time to 1620 micro- 
seconds, a few voltage-induced arcbacks occurred 
(films 13-15). This represents the transition 
region, in which residual ionization is reduced so 
that its effect is comparable to that of voltage. 
Finally, when the speed of deionization was 


VOLUME 13, MARCH, 1942 


10° Interval 
© 


Percent Arcbacks in each 


slightly increased by lowering the pressure to 27 
microris (film 16), with the long period of 1620 
microseconds still allowed, the arcbacks were 
nearly all of the late, voltage-correlated type; 
residual ionization was insufficient to charge the 
particles to the critical voltage, and this charging 
did not occur until high voltage caused glow 
discharge or field emission, or some form of 
voltage-dependent charging. 


7. CONCLUSION 


According to Kingdon’s theory, the condition 
for arcback exists whenever the ion current to a 
negative anode is large enough to charge insu- 
lating particles on its surface to a critical voltage, 
sufficient to cause runaway field-emission. In our 
experiments, when the mercury vapor pressure 


100 





Hg Pressure 33 Microns (64°C) 


13 High Voltage Applied 10 
Electrical Degrees After 
12 Beginning of Commutation 


3° 


Hg Pressure 9 


20 40 60 80 100 140 0 220 260 
Electrical Degrees after beginning of negative voltage —» 


Fic. 9. Distribution of arcbacks in the cycle. Solid line, 
expected distribution at 9 microns calculated from average 
rate. Circles, observed values. The vertical black band 
indicates that all the arcbacks at 33 microns and short 
deionization time occurred in the first 10 degrees. 
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was as high as 33 microns, and the time allowed 
for deionization not more than 815 microseconds, 
this condition was fulfilled at the beginning of 
the period of negative anode voltage; a large 
ion-current flowed momentarily to the anode 
when it became negative, consisting of ions left 
over from the conducting period. In these circum- 
stances, all the arcbacks were observed to occur 
at the beginning. This condition is present in any 
practical rectifier when it is operated at high 
vapor pressure, or at moderate vapor pressure in 
a circuit of low reactance (short time for de- 
ionization), or at any practical vapor pressure if 
the distance between the anode and wall or grid 
is too great (long diffusion distance). 

With a pressure of only 9 microns, the left- 
over ions were gone in 463 microseconds. Hence 
the condition for arcback did not exist at the 
beginning of the negative cycle, and no arcbacks 
were observed at this time. This situation exists 
in practical rectifiers which are properly designed, 
with closely-spaced deionizing grids, and operated 
at moderate vapor pressure (e.g. 45°C Hg tem- 
perature), in circuits having average reactance. 
Arcback is not to be expected under these con- 
ditions until voltage or some other operating 
parameter reaches a critical value, at which a 
new charging mechanism appears. In our experi- 
ments we believe, from visual observations, that 
this new mechanism was glow discharge; al- 
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though other voltage-dependent mechanisms are 
not excluded. Whatever the mechanism, voltage 
appears to have been the parameter responsible, 
since the phase distribution of arcbacks was 
correlated with voltage. Other parameters, how- 
ever, are possible sources ; for example, resonance 
radiation from the discharge to other anodes, in 
the case of multi-anode rectifiers. 

These interpretations are consistent with the 
conclusions of Kingdon and Lawton’ and Pakala 
and Batten,‘ and we believe that they are of 
general application. They are also in agreement 
with earlier investigations pointing to residual 
ionization*® and glow discharge! as the chief 
causes of arcback. We conclude, tentatively, 
that: 

1. The majority of arcbacks, and perhaps all, 
in practical rectifiers, are caused by the charging 
up of small non-conducting particles on the 
anode, as suggested by Kingdon, ‘‘pure”’ arcback, 
from clean metals, occurring only at much higher 
voltages. 

2. The charging up is mainly due to one of two 
mechanisms, viz., left-over ions or glow discharge. 
If operating conditions favor the former, arcback 
will occur at the beginning of the cycle; if the 
latter, it will occur most often near the middle, 
where the voltage is highest. If conditions are 
controlled so as to favor neither, arcback-free 
operation may be expected. 
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The Depth Dependence of Earth Conductivity Upon Surface Potential Data 


W. G. Keck* Ano W. F. CoLsy 
University of Michigan, Ann Arbor, Michigan 
(Received December 16, 1941) 


The problem of determining electrical earth conductivity from the surface, in the case 
where it is a function of the depth only, is solved by a perturbation method, which formally 
at least, allows the unperturbed functions to be perfectly arbitrary. Numerical work, however, 
is restricted to those functions which are available as solutions of the so-called “inverse” 
problem, a group of which is given. Of this group the case in which the unperturbed conduc- 
tivity varies exponentially is treated in detail and two examples showing the success of the 
method are presented. A numerical method of solving integral equations of the Laplace type, 


which occur in the above, is also submitted. 


INTRODUCTION 


N the prospecting for shallow geologic struc- 
tures by electrical resistivity methods, the 
problem of determining earth conductivity as a 
function of depth is of the utmost importance. 
In its simplest form the problem consists of 
finding the conductivity function from a knowl- 
edge of the surface potential around a single 
point electrode. Since more elaborate electrode 
arrangements may be dealt with once the method 
for the point electrode has been worked out, the 
majority of papers on the subject deal with this 
simple arrangement. These papers may be di- 
vided into two groups, depending upon the 
manner in which they approach the problem. 
In the “‘direct’’ method the purpose is to deter- 
mine the form of the conductivity function from 
a knowledge of the potential produced at the 
surface. This is inherently the more difficult and 
has received attention only recently. The so- 
called ‘inverse’ problem assumes that the con- 
ductivity as a function of the depth is known, 
and has as its objective the computation of the 
surface potential. Practical application of the 
inverse method would consist of having available 
a large number of solutions giving the surface 
potentials for chosen conductivity distributions 
and comparing the experimental data to them. 
The assumed conductivity function whose sur- 
face potentials best agree with those obtained 
experimentally is then taken as the correct 
solution. The limitations of such a method are 
obvious. 
The “‘inverse”’ problem is of the type ordinarily 


*Now at Michigan State College, East Lansing, 
Michigan. 
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encountered in potential theory. Thus, means for 
attacking the problem were already available and 
it is understandable that most of the early 
literature on the subject should fall into this 
classification. Hummel! supposed the earth to 
consist of horizontal layers of different thickness 
and conductivity and obtained a solution by the 
method of images so familiar in classical electrical 
theory. Peters and Bardeen,? with the same 
assumption, gave a direct solution of Laplace’s 
equation in terms of Bessel functions. Ehrenburg 
and Watson*® modified Hummel’s method by 
assuming equal thickness of strata, thus obtain- 
ing superposition of images, in an attempt at 
solving the many-layered earth. Roman‘ pub- 
lished numerical tables for computing the poten- 
tial distribution in the case of a two-layered 
earth. Finally, Muskat® made a detailed analysis 
of the formal solutions for a many-layered earth, 
obtaining special formulas appropriate for com- 
putation at large and at small distances from the 
electrode. 

It is to be noticed that all the above references 
deal with a layered earth. The problem of a 
continuously varying conductivity seems to have 
been neglected until Slichter,* although primarily 
interested in the direct problem, indicated the 
solutions for a limited number of cases in the 
form of infinite integrals. 

An early attempt at a “direct’’ solution was 


1}. N. Hummel, Zeits. Geophysik 5, 89, 228 (1929). 
2L. J. Peters and J. Bardeen, Physics 2, 103 (1932). 
3. O. Ehrenburg and R. J. Watson, Trans. A. I. M. 

M. E. 97, 423 (1932). 

‘I. Roman, U.S. Bur. Mines, Tech. Paper 502 (1931). 
5 M. Muskat, Physics 4, 129 (1933). 
6 L. B. Slichter, Physics 4, 307 (1933). 
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made by Tagg’ who prepared a set of charts for 
the two-layered earth, from which he showed 
that one could determine the depth to the 
discontinuity, as well as the conductivity of the 
second layer, from a knowledge of the potential 
measured at two different distances from the 
electrode. Pirson* extended Tagg’s method with 
slight variation to the case of three layers. 
Slichter® and Langer® collaborated on a formal 
solution of the “direct’”’ problem. They conclude 
that a knowledge of the surface potentials alone, 
suffices to determine uniquely the variation of 
the conductivity with depth. Unfortunately their 
method demands the representation of the experi- 
mental data in the form of a reciprocal power 
series. Thus if is likely to be limited in its prac- 
tical application. In fact Stevenson'® has shown 
that in certain cases the method must fail since 
an expansion of the form postulated is not 
capable of representing the experimental data. 
Pekeris," on the assumption of a layered earth, 
modifies Slichter’s method to a determination of 
the depth to the discontinuities involved and the 
conductivities of the various layers. It is a 
graphical method designed to eliminate the 
power series representation mentioned above. 

Stevenson” has developed a method of suc- 
cessive approximations in which the zero approxi- 
mation is that of a uniform earth. He then 
carries out computations for the first approxima- 
tion, showing that although the results obtained 
are not in very good agreement with the true 
solution, he is able to apply the method to 
certain problems which the more exact method of 
Slichter and Langer cannot handle at all. 

The present paper attempts to supplement the 
already existent methods. The basis of the solu- 
tion rests upon a method of perturbations which 
is distinct from that which Stevenson has used. 
In- Part I we shall develop the formal solution 
when the unperturbed functions for the potential 
and conductivity are perfectly arbitrary. In 

‘Part II we shall restrict ourselves to a certain 
class of the unperturbed functions, showing the 
type of integral equations which they yield. Of 

7G. F. Tagg, Trans. A. I. M. M. E. 110, 135 (1934). 

&S. J. Pirson, Trans. A. I. M. M. E. 110, 148 (1934). 

*R. E. Langer, Bull. Am. Math. Soc. 39, 814 (1933). 


10 A. F. Stevenson, Physics 5, 114 (1934). 
"C, L. Pekeris, Geophys. 5, 31 (1940). 
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this class the type in which the conductivity is 
assumed to be exponential wili be analyzed with 
the objective in mind of developing a simplified 
procedure for numerical treatment. In Part III 
a method is proposed for numerically inverting 
integral equations of the Laplace type which 
occur in Part II. Finally, in Part IV, two ex- 
amples are treated numerically. These show what 
success to expect from the perturbation method 
when the actual conductivity is either continuous 
or discontinuous. The numerical inversion out- 
lined in Part III is applied in these problems. 


PART I. SOLUTION OF THE GENERAL PROBLEM 


We will consider the case of an isotropic, 
infinitely extendant half-space containing a point 
source of current at its surface. In this half- 
space the conductivity is assumed to be a func- 
tion of the depth beneath the surface only, with 
the further restriction that it shall not vanish at 
any finite depth. The symmetry in the potential 
function resulting from such an assumption 
makes it convenient to use cylindrical coordi- 
nates with origin at the electrode. The condition 
that the divergence of current be zero every- 
where, except at the electrode, yields us the 
partial differential equation : 

ht 


or? 


ldg o dg Be 


r or 





(1.1) 


where it is understood that the potential and 
conductivity functions are represented, respec- 
tively, by: 


g=¢(r,z) and o=a(z) with o’ =da/dz. 


The above equation must be solved, subject to 
the boundary conditions: 

(a) That ¢ approach zero as either r or z 
approach infinity. (1.2) 

(b) That a disk of small, but finite, radius a be 
assumed as the electrode and is maintained at 
constant potential. Then d¢/dz(r, 0) =0 if r>a; 
and shall have the same value as would be ob- 
tained if the earth were a uniform half-space for 
r<a. This is equivalent to saying that the 
surface layer is uniform to a depth large com- 
pared to a. 

(c) That the surface potential g(r, 0) is known 
everywhere. 
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One may observe here the distinction between 
the so-called “‘inverse’”’ and ‘‘direct”’ problems. 
The former follows the usual pattern of problems 
in potential theory. In it the conductivity func- 
tion is known, which with the differential equa- 
tion governing and a knowledge of the potential 
or its normal derivative at the boundaries, 
enables us to determine the potential every- 
where. In our “direct” problem we have both 
the normal derivative and the potential given for 
all points on the surface. This appears to be an 
overdetermination until we remember that the 
conductivity function o(z) is now also unknown. 
Thus we might suspect that the added boundary 
condition contains the information necessary to 
determine o(z). Indeed Slichter has shown that 
the surface potentials do uniquely determine the 
conductivity as a function of the depth. A com- 
plete solution will then yield us both the poten- 
tial and the conductivity everywhere. However, 
the potential at depth has little interest for us 
and it will be shown that we may solve for the 
conductivity directly from our knowledge of the 
surface potential. 

Equation (1.1) is separable and yields the two 
ordinary differential equations: 


R”+R’'/r+?R=0, 
| Z"+0'Z/o—NZ=0. 


(1 3a) 
(1.3b) 


The proper solving function for (1.3a) is the 
zero-order Bessel function, since the requirement 
(1.2a) must be met. The solution for (1.3b) is not 
available until a choice of o’/¢ is made. 

We will now make use of a perturbation 
method in which it is assumed that 


(1.4a) 
(1.4b) 


¢=¢got ¢1, 
oe /o=t=T0+71, 


where g is the solution of (1.1) corresponding 
to a choice of the conductivity function such 
that oo’/oo=79. Substitution of (1.4) into (1.1) 
leaves us with: 


ag, 1dgy O¢1 
—4— —4 (rot 171) — 
Or? or Or Oz 
O go 0791 
+7r,—+——=0. (1.5) 
dz =z” 


This equation is also separable, yielding the same 
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R part as in (1.3a) and the equation in z: 


Z1"+(ro+71)Z1'—NZ1=—71Z0'. (1.6) 


We shall now express the complete solutions for 
¢o and g; in the form of Fourier-Bessel integrals, 


e) 


ol’, 2) = Ao(A 
oo(r, 2) J () 
sin (Aa) 
x—— Jo(Ar)Zo(a, z)dx, (1.7a) 
Aa 
a(r,s)= f Ax(n) 
* da Qs) 
he z)dX, (1.7b) 
a 


where the normalizing factors Ao(A) and A;(A) 
have been introduced so that Zo(A, 0) and Z,(A, 0) 
become unity. If we multiply the Z equations 
(1.3b) and (1.6) by Z,; and Zo, respectively, and 
subtract, we obtain: 


{Zo'Z1—Z1"'Zo} +10{Z0'Z1—Z1'Zo} 


—1Z)'Zo=71Z0' ZoAo(ar)/A (A), (1.8) 


in which we have used the normalizing functions 
introduced in (1.7). The term 7:Z,/Zo is the 
product of two perturbing functions which we 
shall neglect in comparison to the other terms in 
the equation. Equation (1.8) may then be re- 
written : 


d 
Flic +10{Zo'Z1—Z1'Zo} 
Zz 


=11Z0ZoAo(d)/Ai(A), (1.9) 


a differential equation of the first order, which 
when integrated over all positive z gives us 


« 


g0{Z9'Z1—Z1'Zo} 
0 
AA) ¢” 
= f o071Z0 Zodz. (1.10) 
Ax(A)~ 





From conditions at the infinite boundaries we 
know that Zo(A, ©) and Z,(A, ©) are both zero. 
We may further invoke condition (1.2b) to show 
that 

Ao(A) = —/Cao(0)Z0'(A, 0) 


C=2x/I. 


where (1.11) 
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Finally, if o(0) be set equal to oo(0), which is a 
convenient choice, we may show that 0¢;/d2(r, 0) 
must be zero everywhere. This leads immediately 
to the conclusion that Z;’(A, 0)=0. The above 
results are put into (1.10) to obtain: 


w 
f oo71Z0' Zodz 
0 


= Cao?(0)[Z0'(A, 0) 7A (A) /A. (1.12) 


If we look back to (1.7b) we see that the 
function A,(A) comprises our known information 
since that integral may be inverted to yield 


= f gilt’, 0)Jo(Ar)rdr. (1.13) 
sin Aad y 


At this point it is convenient to allow a to be- 
come vanishingly small, under which condition 
\a/sin (Aa) approaches unity. By substitution 
we then have 


f oo71Z9 Zodz= Cao"(0)[Zo'(A, 0) }* 


0 


xf gi(r, O)Jo(Ar)rdr. (1.14) 
0 


Now remembering that 
/ 4 / 
Ti=d70 /F—do /GO 


and making use of a partial integration of (1.14) 


: a | d 
f log |= —{aZo'Zo\dz 
0 oe) dz 
= Cay?(0)Zo' (A, of gilr, 0)Jo(Ar)rdr, (1.15) 


which gives us our solution in the form of an 
integral equation of the first kind. In this equa- 
tion ¢gi(r,0) represents the experimental data. 
The indicated integration of the right-hand 
member may, as is usually necessary in practical 
applications, be carried out by numerical quadra- 
ture and will yield a function of \ only. The 
unknown conductivity function, which is the 
objective of our analysis, rests under the left- 
hand integral in the form log ¢/o. The kernel of 
the equation, it will be noticed, changes with 
each new choice of the unperturbed conductivity 
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function. In the event that the inversion of 
(1.15) is accomplished the required solution for 
the conductivity is: 


(1.16) 


a(z) =a9(z) exp {log (¢/ao)}. 


PART II 


A. Some Special Unperturbed Conductivity 
Functions and the Integral Equations 
Which They Yield 


In determining the kernel of (1.15) it becomes 
necessary to integrate (1.3b). Unfortunately very 
few solutions of this differential equation which 
involve physically occurring conductivities have 
been reported. Slichter gives a few which occur 
as special cases under a more general group 
which we shall now develop. Along with each 
solution the corresponding integral equation will 
be given. 

If one 


relation 


transforms (1.3b) by means of the 


Y(A, 2) =a0'(2)Zo(A, 2), (2.1) 


we obtain 


1 oo . 
alts.) 
4 00 


which under the condition that 


le AA 
Mek -»|¥=0, (2.2) 
2 Ln) 


1 (a0 /o0)?— 300 /oo= —a” (2.3) 
has the solution in terms of Zo 
oo'Zo=exp [ — (a?+°*)!s J. (2.4) 


In the above oo must satisfy (2.3) which may be 
rewritten to read . 


d?(ao') /dz*+a7o,'=0, (2.5) 
the solution of which is 
oo! = Ae**+ Be-:. (2.6) 


For this case then, (1.15) becomes 


a o 
f log (=)}+| (a*+)?)! 
0 ni) 


Ae**— Be-** ' es 
+a( . ——-) | I natn hedg 
Ae**+ Be-= 
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= —C(A+B)*{(a?-+2?)'(4 +B) 
+a(A -B))*f gi(r,0)Jo(Ar)rdr. (2.7) 


One may specialize (2.6) and obtain consider- 
ably simpler expressions for the above. Consider : 


Conductivity Uniform 


Taking a=0 and oa(0) as unity, 
f log (¢/ao)e~**dz 
0 
1 ° 
= --cf gil(r, O)Jo(Ar)rdr. (2.8) 
2 Yo 


It is interesting to note that Stevenson’s method, 
which employs this approximation, leads him to 
an integral equation which cannot be solved by 
ordinary means. He then transforms this equa- 
tion into the somewhat more tractable equation 
of the Laplace type. If one then performs a 
partial integration of Stevenson’s final equation 
the result is our (2.8). Thus it appears that the 
present method, in addition to providing more 
latitude in the choice of unperturbed functions, 
leads more directly to equations which may be 
treated numerically. 


Conductivity Exponential 
With oo(z) =e? 
i) log (a /ao)e-200* edz = 


| (a?+ ?*)i +a 


1 ” 
soni | f gi(r,0)Jo(Ar)rdr. (2.9) 
0 


2 | (a?-+d2)) 


Conductivity Hyperbolic 
With oo'(z) =2 cosh (az) 


J log (o/o0)[A?+ {(a?+A*)! 
0 


+a tanh (az)} 2 Je-2attr begs 
= —Clat+n) f gilr, O)Jo(Ar)rdr. (2.10) 
0 


VOLUME 13, MARCH, 1942 


If we go back to (2.5) setting 


9 


—(oo')=0 weobtain op=az+b. (2.11) 
dz? 
Then our integral equation becomes 
ra a 2 
f log (o/a9| x4 +——| Jena 
! \"" (@e+0)* 
= — C(b*\?+2ab\+<a?) 
xf gi(r,0)Jo(Ar)rdr. (2.12) 
0 


All of the above integral equations are closely 
related to the Laplace type. The inversion of the 
Laplace integral has been accomplished when the 
integral is known as a function of a complex 
variable. This is of no help in the present 
situation however, since we know the (A) function 
only along the real axis. In Part III we shall 
propose a numerical inversion, not only for the 
equation which is strictly of the Laplace type, 
but also for integrals like the above which are 
closely related. 


B. Application of the Method when the Con- 
ductivity is Assumed to be Exponential 


It will facilitate numerical computations con- 
siderably if we can find a function oo(z) such 
that by simply varying one of its parameters it 
can be made to approximate any conductivities 
encountered in practice with sufficient closeness 
to make our perturbation method effective. 
Since Stevenson has had some success with an 
approximation in which the earth is considered 
of uniform conductivity we have reason to 
expect that 


oo(2) = 09(0)e#" (2.13) 


will answer our needs. It has the desirable 
feature of being able to approximate conductivi- 
ties either increasing or decreasing with depth as 
well as the situation of a uniform earth. 

In what follows it will be convenient to choose 
o (0) as our unit for measurement of conduc- 
tivity. Then according to (2.9) 


f g(z)e~2(a*+a) egg = Gy), (2.14) 
0 
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where 


g(z) =log a/a 
and 
1 (a°+X*)! +a 


G(A) = —- GC . ao f gir, 0) Jo(Ar)rdr. 
2 | (a?+ ?)! 0 


The function G(A) may be determined from our 
knowledge of g;(r,0) by quadrature. It now 
becomes necessary to establish some criterion for 
the selection of a. It must be remembered that 
(2.14) is an approximation which is best when 
oo(z) is most nearly like o(z). To make these 
two functions coincide for all values of z with a 
single choice of @ will, in general, not be possible. 
However we may so choose a that the assumed 
fits the actual conductivity well in that range on 
z in which we are most interested. Correspond- 
ingly then our solution will be most reliable in 
the same range. However, we have no previous 
knowledge of the actual conductivity so we shall 
have to study G(A) for this information. Inspec- 
tion of (2.14) shows that \ controls the interval 
over which the integration is effective. Thus for 
large \ information concerning g(z) is restricted 
to shallow values of depth, while as \ becomes 
smaller the range of integration becomes deeper. 
Thus the reciprocal correspondence between the 
metrics \ and 2z is evident and in the numerical 
inversion process described in Part III it is 
possible to obtain the approximate depth to 
which the integration is effective for any value 
of \. Suppose then that we have decided upon 
the depth to which we wish to determine the 
conductivity. From it we obtain the correspond- 
ing value of \,, and so select a that G(A,,) =0. 
This means then that contributions to the 
integral where g(z) is too small have been offset 
by contributions where it is too large. Due 
allowance must be made of course for the fact 
that the weighting factor in the integral has 
exponential form. Similar reasoning may be 
applied to show that the point or points at which 
the curve for the assumed crosses the true con- 
ductivity will be revealed by the values of \ at 
which 0G/d\=0, since if extension of the effective 
range of integration has not changed the value 
of the integral, we must conclude that g(z) is 
zero at the extremity of this range. 

Applications of the above criteria in the selec- 
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tion of a, since they are somewhat a matter of 
judgment, are facilitated considerably by ex- 
perience. A set of curves showing G(A) against \ 
for various choices of @ is obtained with little 
additional labor and is very helpful. Generally 
it is readily apparent which values of a@ are 
entirely outside the range of consideration. Ex- 
perience has shown that a considerable variation 
of a near the correct choice does not seriously 
affect the result. 


PART III. THE NUMERICAL INVERSION OF INTE- 
GRAL EQUATIONS OF THE LAPLACE 
AND RELATED TYPES 


Consider the integral equation 


D 


{ g(z)e**dz=G(£). 


0 


(3.1) 


This equation may be inverted under suitable 
restrictions to give 


1 ct+in 
g(z) =- -f G(B)e*-dB. (3.2) 
2m c—iw 


A rather complete discussion of the inversion of 
this integral is given by Pipes” in a paper on the 
operational calculus. Unfortunately the method 
requires a knowledge of G(8) as a function of a 
complex variable whereas in our problem we 
know it for real values of 8 only. An inversion 
which can be employed however, assumes that 


x n! 


g(z)=> b,2", G(8)=X b,— 


n=0 n=1 prt! 


(3.3) 


and if one can expand G(@) in reciprocal powers 
of 8 it is merely necessary to identify coefficients 
in the two series. This method has the limitation 
that regions of rapidly varying g(z) cannot be 
well represented by a power series and that a 
large number of coefficients will have to be 
evaluated if any detail in the function sought is 
to be obtained. As an alternative the following 
numerical method of inversion is submitted. 
Assume that an approximate solution of (3.1) 
is available. If we call this solution go(z) we may 


2 Louis A. Pipes, J. App. Phys. 10, 3-5 (1939). 
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write 


{ | 2(%) — go(z) }e-**dz 
=G(8)—Gi(8)=G.(8) (3.4) 


or 


x 


f gi(zje**dz=G,(8), (3.5) 


0 


which is a new integral equation of the same 
form as (3.1). Obviously we may operate upon 
it as we did upon (3.1) to obtain its approximate 
solution gi(z) and a new residual function G2(8). 
This process is to be repeated until the residual 
function G,(8) has been reduced to zero every- 
where, within the experimental uncertainty of 
determining G(8) originally. The solution may 
then be written 


g(z) =go(z) +g1(2) +--+ gn(z). (3.6) 


The number of terms one will have to determine 
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Fic. 1. The two-layered earth. Showing a set of curves 
useful in determining the most favorable value of the 
parameter @ in the unperturbed conductivity function 
e*2?_ Note that certain of the curves show two points at 
which the slope is zero meaning that there are two places 
— the unperturbed conductivity curve crosses the 
actual, 
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in this series depends upon the closeness of the 
approximate solution. 


The Approximate Solution 
Consider Eq. (3.1) to be rewritten to read: 


m 


 ang(2x)e-P**-+-g (2m) f e-btdz=G(6,), (3.7) 
k=1 tm 


where the subscripts k designate the position of 
the ordinates at which the value of g(z) is de- 
sired. This represents a numerical integration in 
the interval 0<z<z,, from which point on, g(z) is 
considered constant and the integral evaluated 
directly. The coefficients a, will be determined by 
the particular scheme of numerical quadrature 
employed. The expression (3.7) is then a set of 
simultaneous equations, equal in number to the 
number of points at which we wish to determine 
g(z). A numerical inspection of these equations 
reveals that the choice of 8, because of the 
exponential decrement, influences each in such 
fashion that it may be terminated at any desired 
value of k. Thus an approximate solution sug- 
gests itself, for with some large value of 8, say 61, 
we might solve (3.7) for the constant value of 
g(z) which would yield G(6:) and ascribe .it to 
g(zi) since this ordinate surely had greatest 
weight in determining it. Proceeding to some 
smaller value of 8 we use the previously deter- 
mined g(z,;) and solve for the constant value, 
g(z2), which would yield G(2) if it were effective 
in the remaining interval. This process is re- 
peated until we have determined all the g(z)s to 
g(Zm). The set so determined we shall call the 
function go(z). Insertion of it back into Eq. (3.6) 
enables us to compute the function Go(8) of (3.4) 
which when subtracted from G(8) leaves the 
residual function G,(8) and the new integral 
Eq. (3.5). The procedure from this point on is 
evident, the process being terminated when 
G,(8) is everywhere less than the uncertainty in 
the original function. 

It should be pointed out that care in the 
selection of the 6’s and the position of the 
ordinates of g(z) will influence the success of the 
approximate solution, reducing the labor in the 
whole inversion process. In the examples pre- 
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Fic. 2. The two-layered earth. 
Curve A represents the true con- 
ductivity. Curve B is the solution 
4 obtained by the present pertur- 
bation method for which curve D 
is the unperturbed conductivity. 
Curve C is the solution by 
Stevenson’s method in which the 
unperturbed conductivity is as- 
sumed to be uniform. 
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sented in Part IV four terms in the series (3.6) 
were sufficient to represent g(z) adequately. 


PART IV. NUMERICAL EXAMPLES 


In the preceding sections it has been con- 
venient to suppose that G(A) is to be obtained 
from an evaluation of the integral over ¢(r, 0). 
However we may obtain it more directly by 
noting that, utilizing an inversion of (1.7) 

D 


i) gilr, 0) Jo(Ar)rdr 


0 
D 


-{ g(r, 0)Jo(Ar)rdr—Ao(A)/X. (4.1) 
0 

This step is often, as is the case for a) exponential, 
a necessity since we cannot obtain ¢go(r, 0) except 
by mechanical quadrature of (1.7). The numerical 
integration over g(r, 0) may of course be facili- 
tated by taking out the part due to a uniform 
earth and integrating separately, as is shown in 
our example. 


The Two-Layered Earth 


_The methods of Tagg or Pekeris could solve 
this problem with less labor than our own, since 
they have been designed especially for it. In our 
own method the requirement is that o(z) be 
everywhere continuous. In applying it to a 
discontinuous case we presume that the potential 
produced is not far different from that for some 
continuous ¢(z) which is everywhere equivalent 
except in the neighborhood of the discontinuity, 
at which place it has merely rapidly changing 
values of conductivity. We should like to check 
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how well this assumption is borne out. In the 
event it is substantiated we shall feel free to 
apply the method to more practical problems, 
complicated by having both continuous and dis- 
continuous combinations. 
By method of images we know that the poten- 
tial due to a two-layered earth is 
1 (1 x k” | 


g(r, 0) 7 ae —+2 
Co(0) |r 


th. 
n=1 (r?+4+4n2h?)! 
where the reflection coefficient k=(¢1—¢2)/ 


(o,+02) and the depth to the second layer is h 
It may then be verified that 





(4.2) 


1) (a?+X*)? +a 
onal etre 
2| (a?+?*)} 
1 2 « 1 
x|—+- Tie. ___..—}, 49 
XN Aeet (a?+22)!+a) 
In the example at hand we have taken k= —1/2 


which means that the conductivity of the lower 
layer is three times that of the upper, while / is 
the unit of measurement of depth. A series of 
curves for G(A) plotted against 1/A appear in 
Fig. 1 from which the selection a=0.1 was made. 
The inversion of (2.14) was then accomplished 
by the method outlined in Part III, a system of 
twelve ordinates being used. The final solution 
for o(z) appears in Fig. 2 as curve B. For pur- 
poses of comparison the same figure contains the 
actual conductivity (curve A), the unperturbed 
conductivity dictated by our selection of a 
(curve D), and the solution obtained by an 
application of Stevenson’s method (curve C). 
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Conductivity a Continuous Function of Depth 


We shall now consider the case where the 
actual conductivity obeys the equation 


o}(z) = Ae“ + Be-™. (4.4) 


By a proper choice of constants this function 
may be made to have a minimum at any desired 
depth. Now if this problem were given us as an 
unknown, we should be free to select a function 
A,(A) as close as possible to the actual A(A) 
observed and thus be led to a oo(z) of the same 
form as (4.4). The application of our method 
would then be quite pointless. As an alternative 
we will again take oo= e?** but apply it somewhat 
differently from what was done in our first 
example. The variation we are about to describe 
will greatly enlarge the field of application of this 
particular function. In fact it is possible that, 
for most practical purposes, it will yield results 
sufficiently accurate so that no other function 
need be considered. This will be an advantage 
since a single, rather simple, set of equations for 
the inversion process will suffice regardless of the 
nature of the true conductivity. 

Proceeding as in the first example it may be 
verified that 


1 (a?+X*)} +a 


G(A) =-; ————_—— 
21 (a2+r2)! | 
| : : | 
x a ———}, (45) 
| (a?+22)!+a(A —B) (o2+22)!+a} 


The constants A =0.2, B=0.8 and a=0.3 were 
taken so that the minimum in o(z) would occur 
at a convenient depth. Again the function G(A) 
was plotted for various values of a. These are 
shown in Fig. 3. We notice that there is no a 
which makes G(A) small over any considerable 
range in \. That this should be the case is not 
surprising, since we are attempting to approxi- 
mate a function which has a minimum by means 
of one which has none. However, we notice that 
successively increasing a@ moves the point at 
which 0G/0X is zero to smaller values of \, 
meaning that the point at which the assumed 
and true conductivities coincide is moving to 
greater depths. The possibility thus exists that 
one might obtain a number of solutions, each 
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utilizing a different a, then build g(z) as a com- 
posite of these individual solutions properly 
weighting each for the range in which it is most 
reliable. However this would require a great 
amount of labor. Approximately the same result 
may be achieved by preparing a composite G(A) 
before performing the inversion. This will consist 
of using a different a for each of the equations 
(3.7) used in the inversion process. As explained 
in Part III, these equations represent integra- 
tions to successively greater depths, depending 
upon the value of \; assigned to each. If, having 
predetermined these values of \; we force a; to 
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Fic. 3. Curves for determining the most favorable a 
for the continuous conductivity case. Note that no curve 
shows more than one place where the slope is zero. 


meet the requirement that 0G(A;)/dA=0 it im- 
plies coincidence of the assumed and actual con- 
ductivities at the lower extremity of the range of 
integration for each equation. These equations 
are then solved simultaneously. Thus, since in 
this process the successive ordinates of g(z) are 
determined from successively deeper integrals, 
we have achieved the situation that each ordinate 
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is determined predominently from an equation 
utilizing its own most favorable unperturbed 
function. 

In Fig. 4 we have plotted the solution for o(z) 
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Fic. 4. Curve A is the actual conductivity for the case 
o' = Ae*?+ Be“? while curve B shows the solution obtained 
by the present method. The dotted curves show the 
successive unperturbed conductivity functions employed 
in building the composite function G(A). 


so obtained (curve B), comparing it with the 
true conductivity (curve A). The dotted curves 
represent oo(z) for the various a’s employed, 
showing the characteristic predicted that the 
successive approximations are in better agree- 
ment to the true conductivity at greater and 
greater depths. 


CONCLUDING REMARKS 


The perturbation method here described, 
though in principle capable of being repeated as 
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often as desired until convergence to the actual 


‘conductivity is attained, is in actual numerical 


computation restricted to the first approxima- 
tion. This is so because one knows only a limited 
number of solutions of the inverse problem. This 
limitation has made it important that oo(z) be of 
such form that, by merely changing a parameter, 
a wide variety of actual earth conductivity dis- 
tributions may be approximated. The reasons for 
selecting an unperturbed conductivity of expo- 
nential form have been given in Part IIB. 

It is quite in order, then, that one ask con- 
cerning the magnitude of the rejected term in 
Eq. (1.9). From the steps immediately following 
Eq. (1.9) it may be shown that what one re- 
jects is 


x 


f 07121 Zodz 


0 
in comparison with 


x 


f ooT Zo Lodz 
0 


which makes it clear that the success of the 
approximation depends upon the interval on \, 
or what is equivalent, the interval on z in which 
we are interested. Now we cannot evaluate the 
above integral directly because of our inability to 
compute Z,(A, 2). However in Part IIB we have 
a criterion for selecting the parameter in o9(z) so 
that the integral above be small or zero for any 
particular value of A, an extension of which 
results in the method of the second example 
where continuously shifting the parameter makes 
the approximation good over a_ considerable 
range in depth. 
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Large Scale Extraction of Polonium from Active Lead and Formation of 
Radioactive Electrode Alloys 


J. H. DILLon anp J. N. STREET 
Firestone Tire & Rubber Company, Akron, Ohio 


(Received January 5, 1942) 


Various possible sources of Po for industrial use are discussed and reasons for selection of 
active lead salts for this purpose are given. Methods of measurement of the amounts of Po 
present on plated metal foils and in electrode alloy are described. The RaD/Pb ratio was 
found to vary from 1.21078 to 1.9 10-8 for various samples of active lead from the Great 
Bear Lake region. Methods for extraction of Po from active lead salts are discussed and a 
large scale plating method is described. Results of experiments on the rate of Po-plating 
from boiling solutions of active lead salts as related to the acid concentration, exposed foil 
area, foil composition, etc., are presented. Chemical methods for producing soluble active 


lead salts are described. 


A method of introducing Po into electrode alloy is described. In a typical experiment, it 
was found that 94 percent of the Po introduced into the melt was retained in the final alloy. 
Contact radiographs of electrode alloy samples showing the affect of various heat treatments 


and surface condition are presented. 


INTRODUCTION 


ESULTS obtained with spark plugs having 

electrodes containing Po have been de- 
scribed in previous communications.'! This paper 
will consist of a description of experiments on 
the extraction of Po from the residues of a 
radium refinery and methods of incorporating 
Po into spark plug electrode alloys. 


SELECTION OF A SOURCE OF POLONIUM 


Polonium is the last radioactive element in 
the uranium series. It has a half-life of about 
138 days and changes into a non-radioactive 
isotope of Pb with the emission of alpha- 
particles at the rate of about 1.810" sec.—! 
gram~!. Po emits only a slight amount of gamma- 
radiation (gamma-radiation energy/alpha-radi- 
ation energy =10~°).2 Hence, it is particularly 
adaptable to an industrial process where danger- 
ous penetrating radiation is undesirable. 

For practical purposes, RaD may be con- 
sidered the parent substance of Po since the 
former is the only long period element (half-life 
= 22.3 yr.) between Ra and Po in the uranium 
series. (See Fig. 1.) Although RaD may be ex- 
tracted from old sealed Ra samples by the use of 
Pb as a carrier, it is much more commonly ob- 
tained from old radon tubes. In the latter case, 


‘J. App. Phys. 11, 291, 561-562 (1940); U. S. Patent 
2,254,169 (1941). 

fk! Bothe and H. Becker, Zeits. f. Physik 66, 307 
(1930). 
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the RaD exists in almost a pure state and the 
Po generated therefrom can be removed with- 
out great difficulty. The methods used for 
extraction of Po from radon tubes include 
chemical separation, electroplating, electrochem- 
ical deposition from solution upon a stationary 
or rotating metal plate, and deposition on a 
hydrogen—charged platinum foil. An excellent 
discussion of these methods with a fairly com- 
plete bibliography is given by Haissinsky.* 

RaD also exists in slight amounts in the waters 
of oceans, lakes, streams, and springs. A few 
calculations based on the known radium and 
radon contents of these waters, however, indicate 
that they cannot be considered as _ practical 
sources of Po. 

Po has been produced by deuteron bombard- 
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Fic. 1. Polonium and RaE present with RaD vs. age of 
RaD. Tp=22.3 years, Tg=5.0 days, Tr =138 days. 


3M. Haissinsky, Trans. Electrochem. Soc. 70, 343 
(1936). 
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ment of Bit-® but the efficiencies of the processes 
thus far employed appear tooo low to warrant 
consideration at this time. 


of RaD is the active Pb 
residues of a radium refinery. RaD and Pb are 
isotopic. Hence, most of the RaD generated 
from the extracted Ra is present with the Pb 
residues since only slight amounts of Rn and 
the other short life products are lost from the 
ore by diffusion and solution in eroding waters. 
Haissinsky* describes several methods of re- 
moving Po from solutions of active Pb residues. 
In general, because of the large amounts of 
active Pb (radio-lead) with which it is necessary 
to deal,. this source has not been considered 
very suitable, especially where very concentrated 
Po preparations are desired for fundamental 
work. Hence, most of this fundamental research 
has been carried on with Po obtained from the 
more convenient radon tube source. The require- 
ments of the particular industrial process here 
described were such that it was not necessary to 
obtain Po in a concentrated form. However, 
large amounts were needed at regular intervals. 
Hence, the active Pb residues appeared to offer 
promise and methods were developed for ob- 
taining Po from them at intervals as it was 
developed from the RaD present with the Pb. 


Another source 
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Fic. 2. Rate of plating vs. plated foil area. Solutions— 
30 g active PbCl, in 1 liter 0.03N HCl (boiling); age 
PbCl.—765 days. 


‘J. J. Livingood, Phys. Rev. 50, 425 (1936). 

5D. G. Hurst, R. Latham, and W. B. Lewis, Proc. 
Roy. Soc. A174, 126 (1940). 

6 J. M. Cork, J. Halpern, and H. Tatel, Phys. Rev. 57, 
371 (1940). 
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METHODS OF MEASUREMENT 


Most of the measurements of amounts of Po 
plated upon foils were made with one of several 
calibrated gold-leaf electroscopes, equipped with 
horizontal parallel-plate ionization chambers at 
atmospheric pressure. Circular Po-plated foils 
were laid upon a larger circular metal disk 
which was placed centrally on the lower plate of 
the chamber. For routine control measurements, 
the plates of the chamber were 1.6 cm apart and 
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765 days; 


the upper plate was maintained at —800 volts. 
The rate of rise of the gold leaf over a range of 
10° to 30° from the vertical was observed 
through a 2-cm diameter window at a distance 
of 20 cm from the leaf. For most laboratory ex- 
periments (see Figs. 2, 3, and 4) another similar 
apparatus was used where the top plate potential 
was —1200 volts and the plate separation could 
be increased to 4.0 cm. This instrument was 
equipped with a microscope with an eye-piece 
scale to facilitate measurements with foils plated 
with small amounts of Po or with thin machined 
plates of Po electrode alloy. 

The gold leaf apparatus was calibrated 
originally by means of a counter which responded 
to each entrant alpha-particle from a Po-plated 
foil mounted in an evacuated tube sealed to the 
‘window end.’’ The counter was similar to that 
employed by Rutherford and Geiger? in their 
early experiments with Ra. 


"7, Rutherford and H. Geiger, Proc. Roy. Soc. A811, 141 
(1908). 
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Fic. 4. Rate of plating in 1200-liter factory tanks, 
60,000 cm? Ni foil surface in tank. Curve A—active 
PbCO; added to HCI of strength to give 0.03N HCl 
boiling sol. of 22.7 kg PbCl.; age of PbCO;—300 days. 
Curve B—22.7 kg active PbCl. added to boiling 0.03.N 
HCl; age of PbCl,—180 days. 


Confirming calibration measurements were 
made with foils plated with Po from 0.05N HC! 
solutions containing crushed radon tubes of 
known age and initial Rn content, kindly 
furnished by Dr. G. Failla.* A curve relating 
alpha-particle emission from these foils to 
electroscope charging rates is given in Fig. 5. 
For the sake of brevity, further discussion of 
electroscope calibrations will be omitted. 

Early control measurements of the Po content 
of electrode alloy were made with the gold leaf 
apparatus as follows: A straightened sample of 
electrode wire or ribbon was cut into a number 
of 6.4-cm lengths and laid in a shallow metal 
tray so as to form a 6.4 cmX6.4 cm-square 
surface covered with alloy wire or ribbon. This 
“pad” was placed on the lower plate of the 
chamber and the rate of leak compared with the 
“zero leak’’ obtained with a similar “pad”’ of 
non-Po wire or ribbon. The differential leak 
thus obtained was translated into absolute units 
from knowledge of counter results on wire or 
ribbon of the same size. The obvious corrections 
for internal absorption in alloy samples were 
made in interpreting counter results. Approxi- 
mate measurements were also made by Dr. F. 
W. Stallmann with electrode wire samples in a 
cloud chamber. 

Later control measurements of the Po content 
of the electrode alloy were made by a photo- 
graphic method. Wire or ribbon samples were 
taped to the surface of a photographic plate 


* Physicist, New York Memorial Hospital. 
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(Eastman contrast lantern slide plate) and 
exposed for a definite period of time. After a 
standard development the densities of the lines 
on the plate were compared against those 
obtained with wires of known Po contents. This 
method was much less tedious than the gold 
leaf method for factory control purposes where 
every coil of electrode wire is tested and was 
convenient for qualitative experimental studies 
(see Figs. 6-9). Contact photographs of this 
type will be termed ‘‘radiographs.” 

A slightly modified Westinghouse ‘“‘Radium 
Hound” is now being used experimentally for 
control measurements of the Po content of wire 
and ribbon. 


CONCENTRATION OF RaD IN ACTIVE LEAD 


The amounts of Po and RaE present with 
RaD of known age (since last separation of Po 
and RaE) can be calculated from the curve of 
Fig. 1. In plotting this curve, values of the half- 
lives of RaD and. RaE of 22.3 years and 5.0 
days, respectively, have been selected. The 
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fifth edition, third supplement. 
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value of 138 days for Po has been taken as the 
average of three recently determined values.* 
The method of measuring the RaD/Pb ratio 
was as follows: A sample of active PbCl, of 
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Fic. 6. 


Exposures of 12-96 hr. of filed and polished 
ribbon. 


Melt A— 47.5 pug Po/g Ni 
Melt B—190 pug Po/g Ni 


known age /, i.e., a sample from which all of the 
Po and RaE had been extracted at time ¢=0, 
was placed in nearly saturated solution in boiling 
0.03N HCI (30.0 g PbCl,/1 0.03N HCl). In most 
cases, the solution was contained in a Pyrex 
beaker. Several identical clean Ni or Cu foils 
were suspended in the boiling solutions for 
consecutive 2—4-hour periods. The Po and RaE 
plated out electrochemically on the foils as 
reported previously by Cook.’ The amounts of 
Po on the foils were measured with the calibrated 
electroscope and a curve of Po plated out as a 
function of time plotted. Typical curves of this 
type are shown in Fig. 2. Where the slope of 
these curves was very nearly zero after a plating 
period of 8 hours, it was assumed that all the 
Po had been removed within a 12-hour plating 
period. The approximate truth of this assumption 
was confirmed by the close agreement between 
measured values of the RaD/Pb ratio for a 
given PbCl, sample with different aging times. 
The amount of RaD present in a given sample 
of PbCl: was calculated from the curve of Fig. 1 
simply by noting the Po/RaD ratio for the 
known age of the sample and dividing the 
measured mass of Po by this Po/RaD ratio. 
Then, from knowledge of the dry weight and 


~ *D. L. Cook, J. Chem. Ed. 11, 313 (1934). 


192 





purity of the PbCl., the RaD/Pb ratio was 
computed. The values of the RaD/Pb ratio 
varied from 1.2*10-° to 1.910-8, for various 
samples from a given locality. The average from 
a large number of samples of active Pb obtained 
from the Great Bear Lake region was 1.6 10-5. 
A calculation based on the Ra and Pb contents 
of the ore, assuming RaD and Ra to be at 
radioactive equilibrium, gave a value of 2.0 
xX 10-5 for this ratio. 


SELECTION OF CONDITIONS FOR LARGE 
SCALE PO PLATING 


Electrochemical deposition on Ni or Cu foils 
from boiling dilute acid solutions of active Pb 
salts was selected as a production process from 
the standpoint of simplicity. It found 
necessary to use glass or glass-lined vessels: for 
the plating solutions. All types of stoneware 
and rubber-lined vessels employed were found 
to adsorb Po strongly during the plating opera- 
tion. Glass-lined tanks of 1200-liters 
capacity, equipped with welded steam jackets, 
were found to be well-adapted to large scale Po 
extraction. 


Was 


steel 


The choice of a Pb salt for Po-plating solutions 
depends primarily on its solubility in dilute acids. 
Only three common salts of Pb appear practi- 
cable, from the solubility standpoint, namely, 
the acetate, nitrate, and chloride. The acetate 
has a very high solubility, 2000 g-I-' at 100°C, 














Fic. 7. Radiographs of filed and polished samples of 
Po alloy ribbon. 
A—Age= 53 days 
E—Age = 258 days 
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and hence appears excellent for this purpose. 
Several Po-plating experiments with saturated 
HC2H;0O¢ solutions of active Pb(C2H302)2 showed 
that the rate of plating was much lower than 
for saturated HCl solutions of active PbCl». 
For example, the rate of plating, measured by 
the apparent RaD/Pb ratio, from a saturated 
0.25N HCsH;Q0s¢ solution of active Pb(C2H;02). 
was less than 20 percent of that of a saturated 
0.25N HCl solution of active PbCl.. Further- 
more, dilution of the Pb(C2H;O2)2 solutions in 
the ratio 10:1 failed to increase the rate of 
plating appreciably. Addition of small amounts 
of HCI also failed to increase the plating rate. 
Further efforts to remove Po more rapidly from 
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Fic. 8. Radiographs of filed and polished flat samples. 
C—Non-Po drawn ribbon 
B—Po alloy drawn ribbon (Melt B) 
D—Poured sample of Melt B 


Pb(C2H 302) solutions, which included the use 
of a Cu disk spinning at 3200 r.p.m. and applied 
voltages on Cu foils, were successful only to a 
small degree. | 

The solubility of Pb(NOs3). is also high 
(1390 g-I-! at 100°C). However, extremely 
dilute HNOs solutions must be used in order 
to avoid rapid dissolution of the Cu or Ni foils 
during the plating process. Since too close 
approach to neutrality of acid solutions of Pb 
salts results in precipitation of the Pb, the use 
of Pb(NOs3)2 has not as yet been found practical. 


RATE OF PLATING OF PO 


The selection of PbCl. as the most convenient 
of these active salts for spontaneous plating, 
from boiling solution imposed definite limits on 
the range of acidities which could be used, since 
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Fic. 9. Radiographs of ribbon of Melt B with heat 
treatments—24 hr. exposure. 
. No heat-treatment 
. Brown oxide (torch) 
Blue oxide (torch) 
. 15 min. torch heating (straw color) 
. 24 hr. at 550°C in furnace 
. 20 hr. at 616°C in furnace 


Due whe 


the solubility of PbCl. decreases rapidly in the 
range (0.0N-1.0N HCl). Hence, the rate of 
plating was studied only in a narrow range of 
acidities (0.005N to 0.11N HCl). These acidities 
were the averages of the initial and final meas- 
ured values at the beginning and end of the 
plating period, respectively. The results of this 
study, plotted in Fig. 3, indicate that the rate 
of plating reaches a maximum at a normality of 
0.025N in the range 0.005—0.11N. The rate for 
the 0.005N solution is extremely low but this 
result has very little meaning since the solution 
became slightly basic after a few minutes because 
of a slight amount of PbO» which was present 
initially with the PbCl. and a large fraction of 
the Pb precipitated as the hydroxide. 

Since it was desired to introduce Po-plated 
foil directly into the electrode alloy melt which 
is essentially Ni, it was desirable to use Ni foil. 
However, the position of Ni in the electromotive 
series suggested that Pb and hence RaD might 
plate out on it as well as Po. Spectrograms of 
scrapings from Po-plated Ni foil surfaces, 
however, failed to show Pb although some Bi 
and Ag were present. Since Pb and RaD are 
isotopic, the absence of Pb indicated that RaD 
was not plated on the Ni. Measurements of 
activity of Ni foils over a period of two years 
failed to show appreciable deviations from the 
known rate of decay of Po. Hence, it was 
concluded that Ni foils could be used without 
danger of abnormal depletion of the RaD 
present with the PbCle. 

The rate of plating as a function of exposed 
Ni foil area under laboratory conditions is shown 
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in the curves of Fig. 2. In obtaining these data, 
a foil was plated for four hours and then removed 
and replaced by a clean foil for the following 
four hours, etc. A curve for Cu is also given to 
demonstrate that, under these conditions, the 
rate is about the same for Cu and Ni. Since the 
rate of plating should be a function of amount 
of agitation and the mean ion-foil distances, it 
is not surprising that the rate of plating for 
corresponding Po surface densities is lower for 
large-scale factory conditions shown in Fig. 4, 
curve B. For factory plating, the foils remained 
in the solutions throughout the entire plating 
period except for short periods necessary for 
cutting samples for test. The rather long factory 
plating period required with this process might 
appear disadvantageous. However, in practice, 
this is not serious since platings can be run over 
night with very little attention except that 
necessary for occasional additions of distilled 
water to the plating tanks. 

The method of preparation of foil surfaces has 
not been found critical in regard to the rate of 
plating. Commercially pure Ni foils seldom 
require cleaning. Etching for a few seconds in 
concentrated HNO; at room temperature has 
been employed in all factory work and for most 
laboratory experiments but is a safety measure 
rather than a necessity. The use of organic 
solvents prior to etching has never been found 
necessary with pure Ni foils employed thus far. 

The rate of plating depends greatly on the 
complexity of the ionic state. For example, 
where active PbCO; was added to HCl of 
sufficient concentration to give 0.03N HCl 
solutions of PbCl2, the rate of plating was much 
greater than for the condition where active PbCl. 
was added directly to 0.03N HCI (see Fig. 4). 
It is natural to assume that the ionic state was 
different for the two conditions. The exact 
nature of this difference has not yet been 
established. Drops of the two solutions were 
_ placed upon a 0.004-cm chlorinated rubber film 
which was laid upon a contrast lantern slide 
plate and exposed for various periods. No 
evidence of radio-colloid aggregates could be 
seen on the exposed areas of the plates, the 
distribution of alpha-particle tracks being very 
uniform for all exposures, when viewed at a 
magnification of 400. 


194 





PREPARATION OF ACTIVE PbCl, 


The several chemical methods which have 
been employed in this work for transforming 
the Pb residues of a radium refinery into useful 
active PbCl, were developed with the coopera- 
tion of C. M. Doede and B. L. Johnson. These 
methods will be discussed only so far as necessary 
to describe the action of the variously prepared 
salts in factory Po plating. 

The active lead was originally obtained from 
the radium refinery of Eldorado Gold Mines, 
Ltd., Port Hope, Ontario in the sulfate form. 
This PbSO, was boiled in a solution of NasCO,; 
and thus converted to PbCO3;. After washing, it 
was converted to Pb(NQ3)2 by addition of HNO; 
and filtered. The hot Pb( NOs). solution was then 
poured into a hot 20-percent solution of NasCOs3. 
The precipitated PbCO; was washed by decan- 
tation until free of sulfate and carbonate ions 
and concentrated HCI was added to the PbCO; 
slurry. The resulting PbCl. was washed until 
just faintly acid and dried. The PbCl,. prepared 
by this process was reasonably pure. The 
process, however, was quite complicated for 
large amounts of PbCls. It was first simplified 
by omitting the last step and simply adding the 
PbCO; to HCI in the plating tanks of strength 
sufficient to convert the PbCO; to PbCl. and 
leave the plating liquor acidity 0.03N. The 
plating was then: made with Ni foils in the 
boiling PbCl. solution (see curve A, Fig. 4). 

It was later found that the active PbSO, 
could be placed in caustic solution and the Pb 
precipitated as PbO, by passing Cle gas into the 
solution. The PbO, was then added to concen- 
trated HCl to form PbCl.. The PbCl. was 
washed so that it was only slightly acid and 
added directly to the 0.03N HCl in the plating 
tanks. In addition to being simpler than the 
earlier methods, this process has another ad- 
vantage in that nearly all the Po present in the 
PbSO, when received from the supplier is carried 
with the Pb and precipitated with it in the oxide 
form. The Po can then be extracted by plating it 
out upon Ni foils suspended in the waste excess 
HCl liquor of the PbO2:—PbCl, transformation. 

A method for precipitating the Pb as PbO: 
directly from the waste caustic liquor of the 
radium refinery was developed by the Port 
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Hope organization of Eldorado Gold Mines, Ltd. 
The active Pb was then shipped in the form of 
PbO:s. The PbO: could then be introduced into 
HCl to form PbCl, and the associated Po 
extracted from the excess HCI liquor. 


LARGE SCALE PO PLATING 


The most practical and efficient large scale 
plating processt thus far developed by this 
laboratory may be described as follows: Aged 
active PbCl, is stored in earthenware crocks. 
At the beginning of each plating cycle, about 
11.5 kg of unaged active PbCl. is in solution in 
1200 1 0.03N HCl contained in each glass-lined 
tank. The HCl is raised to boiling by means of 
steam in the tank jackets and 23 kg of aged 
PbCle is added. Previously etched Ni foils are 
then suspended in the solutions and the boiling 
continued for 12-14 hours. Distilled water is 
added at intervals to compensate for evaporation 
losses. At the end of this period, the steam is 
turned off and 15°C water admitted to the 
jackets. When the liquid temperature has 
dropped to about 25°C, the Po-plated foils are 
removed and the 23 kg of ‘“‘cold-insoluble”’ PbCl. 
is raised from the bottom of the tank and 
replaced in a wet condition in the earthenware 
crock. The Po-free PbCl, is then set aside for 
6-12 months, depending on the aging cycle in 
use. At the end of this aging period, the Po 
developed is again extracted. In the meantime, 
other 23-kg batches of aged PbCl. are run 
through the process at bi-daily intervals. 

A glance at the curve of Fig. 1 shows that the 
aging period should not be too long. Fluctuations 
in Po requirements are easily handled by 
variations in the total number of platings per 
month with consequent variation in the aging 
time. The RaD present with the Pb decreases 
at the rate of about three percent per yr. Hence 
the active PbCl, is considered as machinery 
with a three percent annual depreciation rather 
than as a raw material. The obsolescence period 
for this “machinery” depends, of course, on 
future improvements in the extraction process. 

This large scale plating process appears 
somewhat laborious but is quite satisfactory in 
practice. Consecutive bi-daily platings have been 
carried on by a factory operator with no technical 


t See U. S. Patents 2,254,170; 2,254,171 (1941). 


VOLUME 13, MARCH, 1942 


training for periods of several months without 
interruption. Routine measurements of Po on 
the foils and calculations are made by a trained 
laboratory worker but require very little time. 


THE ALLOYING OF PO IN ELECTRODE METAL 


The usual spark plug electrode alloy consists 
of about 95 percent Ni, 4 percent Mn, and 
1 percent other elements. Hence, the published 
values of the temperature of vaporization of Po 
which ranges between 100°C and 600°C, de- 
pending on the type of metal surface plated and 
the method and medium of deposition,’ are 
rather discouraging from the standpoint of 
forming a Po electrode alloy. However, it was 
found that Po is relatively stable on Ni under 
the conditions of deposition employed in this 
work. For example, a Po-plated Ni foil lost 13 
percent of its activity when maintained at 
600°C for two hours. Subsequent heating for 
18 hours at an average temperature of 535°C, 
followed by 53 hours at 606°C caused no further 
loss of Po from the foil. In other cases, only a 
small loss of Po occurred at 700°C. 

The method finally adopted for producing Po 
electrode alloys may be described as follows: 
Po-plated nickel foils are rolled into a cylindrical 
form and attached to a Ni rod. The Ni-Mn melt 
is made in the usual way in an induction furnace. 
Near the end of the melting period the tightly 
rolled Po-plated foils are plunged into the center 
of the melt. The high frequency current is allowed 
to flow for 30-60 sec. after immersion of the 
foils. The melt is then poured. Production melts 
vary from 90 to 700 kg of electrode metal. The 
weight of Ni foils varies from 1 to 5 percent of 
the total weight of the melt. The metal is 
processed into electrode wire by the usual 
rolling and drawing methods. 

A number of experiments have been made on 
the ‘‘alloying efficiency’”’ with laboratory size 
melts with Po plated on Ni, Cu, and Cu-plated 
Ni foils. In no case was more than 25 percent 
of the Po originally introduced into the melt 
lost in the process. In general, with factory 
melts, the alloying efficiency was higher than 
for laboratory melts. Recent results obtained 
with a small factory melt were as follows: 
Five kg of pure Ni foil with a total surface of 
7.2X10* cm? was plated with Po and RaE in 
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the factory unit. The initial Po on the foils was 
measured with the calibrated electroscope and 
found to be 16.6 yg. The average age of the 
active PbCl, used in the process was 150 days. 
Hence, according to the curve of Fig. 1, the 
RaE/Po ratio at the time of plating was 0.0006 
.0085 = 0.07. Since it has been found that nearly 
all the RaE (at equilibrium with the RaD) is 
plated out with the Po, the effective initial Po, 
calculated by exponential extrapolation to the 
time of plating of the curve giving the actual 
total Po as a function of time after plating, was 
17.8 wg. The foils bearing this initial Po were 
introduced into 87.5 kg of normal electrode 
alloy to form a 92.5 kg melt with a Po concen- 
tration of 192 wug-g- alloy or 1700 wug-cm-* 
alloy (density =8.8 g-cm~*). Measurements on 
the ionization of the alloy were made 125 days 
after the Po on the foils was measured. From 
the decay curve for Po, assuming no Po lost in 
the process, the Po concentration of the alloy 
should then have been 0.54 X 1700 = 918 pug -cm-* 
and the total number of alpha-particles emitted 
-cm7*-sec.-! (N,) should have been 180X918 
= 1.65 X10°-cm~*-sec.' (since 1 pyg Po emits 
180 -alpha-particles/sec.). 

The ionization produced by this sample of 
alloy was measured with the laboratory gold-leaf 
apparatus, equipped with microscope and with 
4-cm plate separation so as to include the entire 
hemisphere of ionization. The apparatus was 
calibrated with a series of foils plated from 
standardized radon-tube solutions. In order to 
make the calibration curve (Fig. 5) very accurate 
within the lower range of ionization intensities 
produced by the alloy, several standard foils 
plated and measured about a year previously 
were introduced and their emissions calculated 
from initial values, corrected for Po decay. The 
“zero leakage”’ (with plate potential applied but 
no foil in chamber) was measured several times 
during the experiment and was found to be 
accurately constant. The average value was 
plotted as a negative rate of charge since, due 
to leakage to the shielding case, the leaf moved 
in the direction of discharge. Two thin polished 
plates of the alloy with exposed areas 15.1 and 
7.1 cm’, respectively, were introduced into the 
chamber and charging rates measured. These 
rates were plotted on the graph of Fig. 5 so 
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that the ionization intensities, expressed as 
alpha-particle emission rates from equal areas 
of plated foils, which would give ionizations 
equal to those observed with the alloy samples, 
could be read off as abscissae. These values were 
divided by the respective exposed areas and 
recorded as Np’. The average Ne’ was 23.0 
full-range alpha-particles-cm~?-sec.~'. 

R. D. Evans" has calculated the number of 
ion-pairs g, produced by alpha-particles from a 
1-cm? surface of a solid as: 


gs= eK R.N,, 


where JN, is the total number of alpha-particles 
emitted in all directions-cm~*-sec.—' within the 
solid, R, is the range in the solid, K is the number 
of ion-pairs produced by the alpha-particle in 
air, and ¢ is a constant having a value 0.142 
for Po alpha-particles. 

The corresponding expression for an extremely 
thin layer of Po on a 1-cm? foil surface is 


qr= KN», 


where Ny is the number of alpha-particles 
emitted in the hemisphere above 1 cm? of 
exposed foil surface. For the ionization produced 
by an alloy sample to be equal to that produced 
by a foil of equal surface area, g,=qr and, by 
definition, Ne=N,’. It follows that 


N,=(1/eR,)- Nr’. 


The value of R, has been determined for RaC’ 
alpha-particles in Ni by Traubenberg" as 
0.00184 cm. Multiplying this value by the ratio 
of the ranges in air of Po and RaC’ alpha- 
particles, the range of Po alpha-particles in Ni 
is 0.00104 cm. This procedure appeared justified 
on the basis of results on the variation of 
stopping power along the range of RaC’ alpha- 
particles in Cu.” This value of R, was used for 
the electrode alloy since it is 95 percent Ni. 
Thus, the experimental value of N, for this 
sample of electrode alloy was 


N,=(1/0.142 X 0.00104) X 23.0 
= 1.56 105-cm-*-sec.~! 


10 R. D. Evans, Phys. Rev. 45, 29 (1934). 

1H. Rausch von Traubenberg, Zeits. f. Physik 2, 268 
(1920); Physik. Zeits. 21, 588 (1920); 22, 587 (1921). 

2 E. Marsden and T. S. Taylor, Proc. Roy. Soc. 88, 
443 (1913). 
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and the corresponding Po concentration was 
867 pwug-cm-*. Thus, the fraction of the Po 
introduced which was retained in the alloy was 
867 /918 =0.94. 

Radiographs obtained by the procedure de- 
scribed under Methods of Measurement are given 
in Figs. 6-9. For control work, Po alloy ribbon 
and wire are simply cleaned with fine abrasive 
cloth and placed on the photographic plates. It has 
been found, however, that the cold drawing pro- 
cess occasionally produces slight apparent aggre- 
gation of the Po in the surface of the ribbon or 
wire. This phenomenon probably results from 
imperfect incorporation of the pure Ni foils in the 
Ni-Mn melt. Thus, tiny inclusions of Po-bearing 
Ni with hardness slightly different from that of 
the master Ni-Mn alloy respond to the drawing 
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Fic. 10. Microradiograph of filed and polished sample 
of Po alloy ribbon. 


Magnification—250 x , Exposure—33 min. 


stresses in a manner differing slightly from that 
of the main body of the alloy. The resulting 
“grainy’’ appearance of the radiographs obtained 
with unpolished drawn ribbon and wire is 
occasionally visible to the naked eye. It does not 
appear when the surface is filed or ground to a 
depth greater than 0.006 cm. No variation in 
apparent Po content of ribbon or wire, as judged 
by radiographs has ever been observed over the 
entire length obtained from a given melt.** In 
_ **It is obvious, of course, that the incorporation of Po 
into certain elements of an alloy can be applied to the 
study of the behavior of those elements in the melting, 
rolling, and drawing processes (see work of G. Tammann, 
Zeits. f. Electrochem. 38, 530 (1932), and G. Tammann 


and H. von Lowis, Zeits. f. Anorgan. Allgem. Chem. 205, 
145 (1932). 


VOLUME 13, MARCH, 1942 


making laboratory studies, however, the pre- 
caution of filing and polishing to a depth of 
about 0.01 cm is observed. The samples used in 
obtaining the radiographs of Figs. 6-9 were all 
treated in this manner originally; that is, before 
other treatments were imposed. 

The radiographs of Fig. 6 were obtained with 
ground electrode ribbon of two melts A and B, 
exposed for periods ranging from 12-96 hours. 
The initial amount of Po introduced into melt B 
was 4 times that of melt A. Comparison of the 
48 hours and 96 hours exposures of melt A with 
the 12 hours and 24 hours exposures of melt B, 
respectively, illustrates that this method is 
effective for approximate comparison of alloys 
of differing Po concentrations. This principle is 
also illustrated in Fig. 7, where two alloys of 
the same initial Po concentrations but of 
different ages were used. Calculations based on 
the known rate of decay of Po gave the Po 
concentration at the time of measurement of 
melt A as about 3 times that of melt E. Com- 
parison of the 48-hours exposure of melt A with 
the 144-hour exposure of melt E indicates 
approximately equal densities, consistent with 
the calculation. 

It should be mentioned that reasonable caution 
must be used in comparing samples of differing 
Po concentrations since alpha-emission pho- 
tography is particularly susceptible to the 
phenomenon of photographic reversal. By avoid- 
ing long exposures of intense sources, however, 
photographic reciprocity can be maintained. 

Examination of Fig. 8 shows that no appreci- 
able loss of Po occurs in the annealing processes. 
A sample of ribbon of melt C containing no Po 
was also placed on the plate to show that the 
blackening is due entirely to the alpha-particles 
of Po. The effects of various types of heat 
treatments are shown in Fig. 9. In this case, 
2.5-cm lengths or ribbon were filed and polished 
and then heat-treated as indicated. They were 
then cut into 1.25-cm lengths, one of which in 
each case was filed down to 3 the original 
thickness whereas the other was given no further 
treatment. Thus, the effect of heating on the 
surface layers is given in the top row of radio- 
graphs whereas the effect at the centers of the 
specimens is given in the lower row. It is clear 
from these radiographs that oxidation of the 
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surface of the alloy (produced in a few seconds 
by means of an oxygen-natural gas torch) did 


not reduce the alpha-particle emission (see 
radiographs 2 and 3 of Fig. 9). Prolonged 
heating in a muffle furnace appeared to produce 
a slight loss of Po in the surface layers (radio- 
graphs 5 and 6). A very severe heating at a 
temperature just below the melting point of the 
alloy (radiograph 4) eliminated the Po from 
the surface layers but produced no appreciable 
loss from the interior. It should be noted that 
none of these heat treatments exactly reproduces 
conditions of spark plug electrode service. 
Surface irregularities developed in actual elec- 
trode service prevent uniform contact with the 
photographic plates. Hence, results obtained 
with electrodes cut from used plugs are not as 
definite as would be desired. In general, however, 
it has been found that little, if any, loss of Po 
from the electrodes occurs in service. 

A “microradiograph”’ of an exposed region of 
a “fine-grain alpha-particle plate,’’ furnished 
through the courtesy of Eastman Kodak Com- 
pany is given in Fig. 10. A length of filed and 
polished Po alloy ribbon was laid on the plate 
for 33 minutes. The exposed region of the 
developed plate was photographed at a magnifi- 
cation of 250X. The dispersion of individual 
tracks is very uniform. Many of the tracks are 
not in focus because of their varying depths in 
the finite layer of emulsion. A similar micro- 
radiograph (Fig. 11) taken with unpolished 
Po-alloy ribbon at a magnification of 500, 
shows an extreme case of aggregation of Po in 
the drawn surface of the ribbon. The dark 
shadow across the plate resulted from incorrect 
operation of the camera shutter. 


DISCUSSION 


No attempt has been made in this paper to 
include all the experiments involved in this 
development which has extended over a period 
of some seven years. Only those results which 
appeared of practical significance have been 
described. A number of laboratory experiments, 
designed to explain the Po deposition phenomena 
observed on the basis of fundamental data in 
the literature, were made but are not reported 
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here since this would necessarily involve a 
rather lengthy discussion. It might be said that 
the ‘‘engineering approach”’ was used where the 
complexity of the materials employed excluded 
purely scientific studies. 
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A 
Fic. 11. Microradiograph of unpolished drawn surface 


of Po alloy ribbon. ; 
Magnification—500 x, Exposure—25 min. 








Efforts to develop methods for utilizing a 
greater fraction of the available Po in the spark 
plug are being made but all aspects of this 
program have not been worked out and conse- 
quently it does not warrant consideration at 
this time. 
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Streamering in Negative Corona 
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Corona from sharp points in air can be divided into two ranges, a diffuse range and a stream- 
ering range, this division being applicable for each polarity of discharge although the details 


are different for the two polarities. 


ORONA from sharp points in air can be 
divided into two ranges, a diffuse range and 
a streamering range. 

In the diffuse range of the negative point 
discharge, the current at onset is characterized 
by pulsations! with pulses of the order of magni- 
tude of 10-'° coulomb. With increasing current 
the pulsations become more rapid. When the 
interval between pulses is short compared with 
the time required for the ions to cross the gap, 
the current becomes much steadier as observed 
with an oscillograph. 

The initial onset pulses produce a slight dis- 
turbance at broadcast frequencies on a radio 
with loop antenna about two feet away from the 
gap. As the current is increased to about 10-7 or 
10-® ampere, the radio interference becomes 
much more intense, and then with still further 
increases in current, the interference goes back 
to approximately the same level as that at onset. 

As the voltage and current are further in- 
creased, there occurs an abrupt transition to a 
new regime, in which the current as shown on an 
oscillograph is fluctuating irregularly 10 to 100 
times as much as anywhere in the diffuse range, 
and the radio interference is much more intense. 
This higher current regime will be called the 
streamering range. 

Decreasing the radius of curvature of the point 
increases the value of current at which the 
transition from diffuse to streamering occurs, 
down to a radius of 0.01°™ but further decreases 
to 0.0008 radius of curvature do not increase 
the transition current appreciably. 

The use of a high resistance behind the point 
has no effect on the transition current, but if the 
point consists of a short wire segment whose 
electrical capacity divided into the 10-!°-coulomb 


'G. W. Trichel, Phys. Rev. 54, 1078 (1938); L. B. 
Loeb and A. F. Kip, J. App. Phys. 10, 142 (1939). 
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pulsation charge gives a voltage fluctuation of 
the order of magnitude of the total voltage 
across the gap, such a ballast resistance can 
diminish the violence of the streamering regime, 
considerably. 

The value of the transition current is greatly 
affected by the shape and composition of the 
collecting electrode. The discharge from a point 
towards a toroidal ring is diffuse to a higher 
current than towards a plane. The discharge 
from a row of points towards a symmetrically 
placed pair of rods is diffuse to a higher current 
than towards only one of the rods. In a row of 
points, the total maximum diffuse current is 
greater than from one point but the current per 
point is less. Bringing points closer to each other, 
in a row, reduces the diffuse current limit per 
point. 

The use of semi-conducting material in the 
collecting electrode increases the diffuse current 
limit, providing the material is 0.01°™ thick or 
less, over insulating backing, and has a resistivity 
of about 150 megohm centimeters per centimeter. 

A non-uniform resistivity or a dust speck on a 
metallic collecting electrode greatly reduces the 
diffuse current limit. In measurements of ozone 
production? ijn corona, it was noted that a speck 
of dust on the collecting electrode abruptly in- 
creased the rate of production of ozone ten times 
or more, at constant current. 

In an experiment to examine the effects of 
back ionization, a point was arranged to dis- 
charge towards a toroidal ring, on the opposite 
side of which there was placed another discharge 
point along the axis of and directed back towards 
the main discharge point. The auxiliary point 
was connected to the source of potential on the 
ring through a ballasting resistance which could 


2 W. H. Otto and W. H. Bennett, J. Chem. Phys. 8, 899 
(1940). 
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be varied. Ballasting with about 100 megohms 
caused streamering to set in at a lower current 
than with the auxiliary point disconnected, and 
when streamering began, the current increased 
several-fold for a constant applied potential. 
Ballasting with one megohm or less resulted in 
spark-over to the auxiliary point instead of 
streamering. 

The transition from diffuse to streamering 
discharge occurs at a much lower current from 
the edge of a ribbon 0.002°™ thick by 0.3°" wide 
than from wires 0.002°™ or 0.005° in diameter. 

In each of these experiments, the feature of 
the system which caused a reduction of the 
diffuse current limit was either the production of 
ionizations near the collecting electrode by 
locally intense electric fields, or too slow a de- 
crease of field with distance from the point ‘in 
the vicinity of the electron attachment sheath. 

The sheath of negative ions over the top of 
the point regulates the electric field underneath 
in the ionization sheath and thus controls the 
current. A partial neutralization of the space 
charge of the negative ions correspondingly 
intensifies the field underneath and serves both 
to let the electrons evade attachment until they 
have moved farther away from the point, and 
also it serves to increase the number of ioniza- 
tions exponentially with the increase in field. 

There is a certain amount of photo-ionization 
generally produced in the vicinity of the sheaths 
by the radiation from the ionization sheath but 
this is diffused throughout the negative ion 
sheath. However, the positive ions originating 
near a spot on a collecting electrode would 
arrive along a channel through the negative ion 
sheath, because of the converging form of the 
field. The additional electrons produced by ion- 
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izations underneath the channel would also move 
out along a region where electron attaching 
collisions were postponed so that when they do 
become attached, their greater distance from the 
point makes them less effective for suppressing 
the electric field underneath them. 

In this way, a small pulse of positive ions 
entering the sheath can produce a much larger 
momentary increase of current of negatives. If 
the increase in negatives in turn produces another 
pulse of positives, large fluctuations in cur- 
rent can perpetuate themselves, once they get 
started. 

The fact that the attachment, occurring.at a 
greater distance from the point renders the 
negatives less effective for suppressing the elec- 
tron field underneath, introduces an essential 
instability of the channelled path of positives as 
compared with the uniformly diffused positive 
ion current from photo-ionizations which occurs 
in the diffuse range. The channelled form occur- 
ring in the streamering range is not necessarily a 
self-sustaining spark or brush-discharge channel, 
but may be considered a preliminary form of it. 
The term streamer is here introduced to distin- 
guish between the geometries of the distributions 
of the positives in the two ranges of negative 
corona discharge. 

The electric field near a ribbon varies inversely 
approximately with the square root of the dis- 
tance instead of inversely with the distance as is 
the case with a point or wire. This less rapid 
decrease of field with distance for the ribbon 
makes any interfering action of positives more 
effective the distance effect, and 
experiment shows that the diffuse current limit 
is much lower for a ribbon than for a wire with 
the same curvature. 


because of 
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On Relaxation Effects in Amorphous Media 


ROBERT SIMHA 
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(Received December 19, 1941) 


The Maxwell relaxation theory leading to the concept of a viscous-elastic body is briefly 
stated. Extensions to systems built up from several such viscous-elastic components with 
different relaxation times are discussed. In the case of high molecular weight materials, it is 
more appropriate to assume a continuous distribution of relaxation rates. Equations are 
formulated whose solution gives the distribution function in terms of the rate of deformation 
and of the applied stress or the elongation in terms of the two other quantities. The connection 
with the corresponding distribution of dielectric relaxation times, molecular inhomogeneity 
and chain branching, and effects at mechanical frequencies large compared with the relaxation 


rates, are pointed out. 


HE investigation of substances whose me- 

chanical behavior follows neither Hooke’s 
elastic law nor Newton’s viscosity law has 
received considerable attention during the last 
hundred years. Deviations may consist in a 
non-linear relationship between applied stress 
and resulting deformation in elastic bodies and 
in a non-linear relationship between applied 
stress and resulting rate of deformation in a 
medium with finite viscosity. Empirical equa- 
tions for such cases have been given by Bing- 
ham,! Rabinowitsch,? Reiner,’ Szegvari,‘ Weis- 
senberg®> and others, while the mathematical 
theory of non-linear strain has been developed 
by Murnaghan.® In the following, however, we 
shall not deal with phenomena of this kind, but 
with those arising from particular time effects 
as collected under the names of relaxation or 
elastic after-effect. In an ideally elastic body, 
the deformation is entirely determined by the 
stress, in such a way that a constant load results 
in a time independent unique value of the 
deformation. In a purely viscous material on 
the other hand, as characterized by a constant 
value of the viscosity, a constant stress results 
in a deformation increasing linearly with time. 
The behavior of amorphous bodies follows 
neither of these patterns. An attempt to account 
for this behavior is represented by Maxwell's’ 


in E. C. Bingham, Fluidity and Plasticity (New York, 
22). 

? B. Rabinowitsch, Zeits. f. physik. Chemie A145, 1 (1929). 

> M. Reiner, Kolloid Zeits. 50, 199 (1930). 

* A. Szegvari, Zeits. f. physik. Chemie 108, 175 (1924). 

°K. Weissenberg, Kolloid Zeits. 46, 277 (1928). 

°F. D. Murnaghan, Am. J. Math. 59, 235 (1937). 

7J. C. Maxwell, Phil. Trans. 157, 49 (1867); Phil. Mag. 
35, 134 (1868). See also the excellent account given in 
Burgers’ First Report on Viscosity and Plasticity (Amster- 
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well-known approach to a relaxation theory of 
viscosity. According to these ideas, a system 
held under constant deformation will, in the 
course of time, due to the heat motion of its 
constituent particles, lose potential energy. This 
is equivalent to a decrease of the internal 
stresses, in other words, these stresses relax. In 
a system held under constant external stress, 
on the other hand, the deformation will continu- 
ously increase in order to compensate for the 
relaxation of the internal stresses. This is 
equivalent to viscous flow. Assuming in particu- 
lar that the relaxation follows an exponential 
time law, and that it occurs for any value of the 
stress larger than zero (no residual stress), 
Maxwell postulated : 


dx 1dX 1 
aT Ef (1) 
dt Gdt Gr 


G is the shear modulus, x denotes the strain, 
X the corresponding stress. The first term on 
the right-hand side reproduces the elastic part 
of the volume element, while its viscous behavior 
corresponding to a viscosity Gd is represented 
by the second term. According to (1), the 
influence of the elastic component vanishes if 
the load is applied for a sufficiently long time 
period, i.e., for a period large compared with 
the relaxation time \. This may be actually so 
in the case of fluids which to a certain extent 
possess elastic properties, like certain solutions 
of high polymers at high concentrations. How- 
ever, silk threads investigated by Weber,® glass 
dam, 1935); Second Report on Viscosity and Plasticity 


(New York, 1939). 
8 W. Weber, Pogg. Ann. 54, 1 (1841). 
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fibers examined by Kohlrausch® and Wiechert,'® 
and recently extensively by Taylor and _ his 


collaborators," and finally data of various 
investigators” on different natural and synthetic 
materials show that Maxwell’s relation (1) is 
not capable of reproducing the actual results. 

Two ways of improvement have been visual- 
ized. One may assume with Weber and Maxwell 
that the relaxation time \ is not a constant 
characteristic for the material, but depends upon 
the stress. This is to be expected if values of the 
stress yielding deviations from the linear laws 
are considered. For other cases, a second possi- 
bility has been suggested which appears more 
promising. Thomson" and especially Wiechert" 
were the first to assume that the behavior of 
the amorphous substance under stress can be 
characterized by a whole set of relaxation times, 
giving a whole set of stresses which combine to 
a resulting internal stress. An interpretation of 
the above mentioned experimental results could 
be given on this basis. This theory, as shown by 
Wiechert, includes Boltzmann’s" formulae for 
the phenomenon usually designed as elastic 
after-effect. 

Recently, this concept of a multitude of elastic 
relaxation times has been reconsidered by several 
authors. Thompson,'® Gemant,'? and Burgers’ 
have discussed various mechanical and electrical 
models, consisting of elastic springs and pistons 
exhibiting friction, and of capacities and re- 
sistances respectively. These models are illus- 
trative of the behavior of viscous-elastic matter 
and serve to derive the differential equation 
connecting deformation and stress. A hydro- 
static model has been given by Taylor." Further- 
more, Isakovich, Leontovich and Frenkel and 
Obrastzov'® have discussed the propagation of 


°F. Kohlrausch, Pogg. Ann. 119, 337 (1863); 158, 337 
(1876); 160, 225 (1877). 

"10 E. Wiechert, Ann. d. physik. Chemie 50, 546 (1893). 

" Nelson W. Taylor, J. App. Phys. 12, 753 (1941). 
Nelson W. Taylor and Robert F. Doran, J. Am. Ceram. 
Soc. 24, 103 (1941), where further references are given. 

“See for instance S. S. Kistler, J. App. Phys. 11, 769 
(1940). 

J. J. Thomson, Applications of Dynamics to Physics 
and Chemistry (London, 1888). 

4E. Wiechert, Dissertation, Kénigsberg i. Pr. 
Ann. d. physik. Chemie 50, 335 (1893). 

4% L. Boltzmann, Pogg. Ann. Erg. Bd. 8, 624 (1876). 

‘6 J. H. C. Thompson, Phil. Trans. A231, 339 (1933). 

17 A. Gemant, Trans. Faraday Soc. 21, 1583 (1935). 

18]. Frenkel and J. Obrastzov, J. Phys. Acad. Sci. 
USSR 2, 131 (1940), where further references are given. 


1889; 
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waves in such media. However, none of these 
modern investigators have explicitly formulated 
equations holding for more than two or three re- 
laxation mechanisms at the most. 

Kuhn'® expressed the Young modulus E in 
the form of a sum 


E=)>); Ee (2) 


to be extended over all relaxation times \,; 
present. E; the moduli as observed at 
intervals of time which are small compared 
with all occurring relaxation times. It follows 
from this equation that the value of E measured 
after a time ¢ comparable with some of the \; 
present will be smaller than the initial value. 
On this basis Kuhn concluded that in the case 
of rubber at room temperature for instance, the 
small value of & (order of 10° dyn/cm? as com- 
pared with 10" dyn/cm? for crystalline bodies) is 
due to the presence of a relatively large number 
of small A; whose contribution to E cannot be de- 
tected in an ordinary extension experiment be- 
cause its duration is too great in respect to these 
relaxation times at room temperature. The differ- 
ence in this connection between, e.g., glass and 
rubber lies in the fact that the relaxation times in 
glass do not differ from each other by so many 
order of magnitudes as in the other case. In other 
words, if we plot the distribution of relaxation 
times as a continuous curve, we should obtain 
a wider spread in the case of rubber and similar 
high polymers. 

As mentioned briefly before, Maxwell’ has 
put forward a molecular picture of the relaxation 
process. The applied stress causes a shift of 
particles into positions of higher potential energy 
corresponding to the stress energy stored up 
in the body. Due to the thermal motions there 
is a finite probability for each molecule to leave 
its position and to jump into a more stable one. 
Such a repeated decrease of potential energy on 
a molecular scale can ultimately lead to a 
decrease of the stress energy of the body. The 
internal stresses will relax if the deformation is 
kept constant. The extent of relaxation will, at 
a given temperature and external force, depend 
upon the relative stability of the various mo- 
lecular configurations. The relaxation time will 


are 


19 W. Kuhn, Zeits. f. physik. Chemie B42, 1 (1939). 
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be connected according to this concept with the 
mean lifetime of the particle or group of particles 
in the equilibrium position shifted by the stress. 
Its magnitude can be estimated with the aid of 
a relation by Frenkel :?° 


A= oe” ae (3) 


where U is the activation energy necessary for 
a jump from one position into another one, and 
do denotes the period of vibration of the molecule 
in the equilibrium position. We assume for the 
sake of simplicity that the temperature de- 
pendence of Xo is negligible. Setting, for instance 
Ayo =10-" sec., U=18 kg cal per mole, we find 
\~1 sec. at T=300°K.*! 

A series of relaxation times in high polymers 
instead of a single one can be due to two causes: 
first to the inhomogeneity of the molecules in 
respect to size which is almost always present; 
second, and more important perhaps, to various 
mechanisms of relaxation in an_ otherwise 
homogeneous assembly of molecules as character- 
ized by varying values of A» and U. Both will 
be determined by the type and intensity of 
molecular interactions and also by the kinetic 
effects connected with the statistically deter- 
mined shape of a long chain molecule possessing 
internal flexibility.22 This results in a multitude 
of frequencies 1/X» and energies U. 

Kuhn's considerations, however, do not lead 
to quantitative relations between the observable 
changes of quantities like the deformation and 
the relaxation times. This question was con- 
sidered recently by Bennewitz and Rétger™* and 
Holzmiiller and Jenckel.** The first two authors 
formulated and solved the equations of motion 
in one dimension of a body consisting of N 


20 J. Frenkel, Zeits. f. Physik 35, 652 (1926). E. Jenckel, 
Zeits. f. physik. Chemie A184, 309 (1939) finds that the 
temperature dependence of the two relaxation times of 
glassy selenium does not follow the simple exponential 
law (3). 

2tN. A. de Bruyne, Proc. Phys. Soc. 53, 251 (1941), 
derives an expression for U in terms of the ultimate co- 
hesive strength and the external stress, thus leading to a 
dependence of \ upon X as previously discussed. 

2,W. Kuhn, Kolloid Zeits. 68, 2 (1934); E. Guth and 
H. Mark, Monatsh. Chem. 65, 93 (1934); J. Frenkel, 
Acta Physicochim. USSR 9, 235 (1938). 

23K. Bennewitz and H. Rétger, Physik. Zeits. 40, 416 
(1939). 

*4W. Holzmiiller and E. Jenckel, Zeits. f. physik. Chemie 
A186, 359 (1940). 
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Maxwell-type species, obeying an equation: 


dK; dx; 
—=a;——y7;K;; j=1,2---N. (la) 
dt dt 
The a; are elastic moduli multiplied by a 
form factor o with the dimensions of an area if 
K; denotes the force. yj=1/\; is the rate of 
relaxation of species 7. In this form, the Eq. (1a) 
holds for extension, while the original Maxwell 
Eq. (1) is formulated for shear. However, if we 
assume that Poisson’s ratio yp is not influenced 
by the relaxation effect, the procedure is justified. 
Equation (la) together with the equations of 
motion: 
d*x 


x ON 
m—+)> K;=K,, (4) 
dt? i=1 


where K, denotes the external force, determine 
the displacement x as function of time. In this 
way, the above authors found for x, if the inertia 
is neglected and if a constant load is applied: 


K. N-1 
x= “|i y 2 s(leu) +rt] (5) 


aj ye 


where s;, g;, and r depend upon the relative 
frequencies of occurrence of each species. 

The first term corresponds to a purely elastic 
reaction with an elastic modulus E: 


N 
z a;=cE. 
j=l 
It is predominant at a time ¢ negligibly small as 
compared with all \;. A measurement at a time 
t larger than all A; allows a determination of 
the last term. It represents a flow process with 
a viscosity 
1 N E; 
ie mora é4 — 
2(1+ 4) i=1 y; 
if we consider the relation between shear 
modulus G; and E;: 
E, 
G;=———_. 
2(1+y) 


The second term gives the transition between 
the two extreme cases of a purely elastic and 
purely viscous behavior, with all y; equal to 
zero and infinity respectively. The sum, there- 
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fore, measures the elastic after-effect. Figure 1 
shows a typical deformation-time curve as 
obtained from Eq. (5). The second half repro- 
duces the elastic recovery, that is the gradual 
decrease of the strain if the applied load is 
suddenly removed. 

The equations have been explicitly solved for 
the case of two and three species respectively. 
In the case of a greater number of different 
relaxation rates, an evaluation of the coefficients 
s; and g; becomes too cumbersome. In view of 
our previous discussion, it appears desirable, 
however, to consider such a case in a study of 
the plasto-elastic properties of a substance built 
up from a polydisperse mixture of giant mole- 
cules. For this purpose, it will be appropriate to 
describe the distribution of relaxation times or 
rates by means of a continuous function ¢(y) 
such that g(y)dy gives the fraction of species 
with relaxation rates between y and y+dy and 
having an elastic modulus Eg(y)dy. Therefore 


| g(y)dy=1. 
0 


Such a procedure was adopted by Fuoss and 
Kirkwood™ in their studies of dielectric dis- 
persion in polyvinyl chloride. In this case, 
the relaxation time is according to Debye’s 
theory*® entirely determined by the size and 
shape of the polar molecule and the temperature 
and internal friction of the solvent in which the 
solute tends to become oriented by the electric 
field. The above investigators were able to show 
that a linear chain molecule exhibiting Brownian 
movement of its parts,” necessarily possesses a 
multitude of electric relaxation times. They 
furthermore formulated an integral equation 
whose solution gave the distribution function if 
the dielectric dispersion curve was known and 
could be expressed analytically. 

While the mechanisms operative are not 
identical in the mechanical and electric case, 


‘because in the latter we are dealing with rota- 


tional movements only, and in the former there 
are also dilatational movements,”’ one can expect 


% Raymond M. Fuoss and John G. Kirkwood, J. Am. 
Chem. Soc. 63, 385 (1941); J. Chem. Phys. 9, 329 (1941). 

26 P. Debye, Polar Molecules (New York, 1929). 

27G. B. Jeffery, Proc. Roy. Soc. A102, 163 (1923); 
R. Simha, J. Phys. Chem. 44, 25 (1940); J. App. Phys. 13, 
147 (1942). 
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a connection to exist. For in both processes the 
rate of relaxation is characteristic for the reaction 
of the assembly of molecules acted upon by an 
external field. Hartshorn, Megson, and Rushton?* 
used an average value for the dielectric relaxation 
time of several phenol derivatives formaldehyde 
plastics, as determined by means of the simple 
Debye theory and calculated therefrom reason- 
able values for the elastic modulus from Eq. (1). 


MATHEMATICAL DEVELOPMENT 


The following equations are derived for the 
calculation of g(y) from a knowledge of the 
deformation-time curve. As in the development 
of Eqs. (1a) and (5), it is assumed that one and 
the same Poisson constant is characteristic for 
each relaxation species. Furthermore, the inertia 
term in the equations of motion is neglected. 
This is justified if the rate of change of the 
external force is small compared with the 
relaxation rates. An integral equation for ¢(y) 
is found and a solution in terms of the acting 
force and of the rate of elongation is given. This 


Fic. 1. 


can be done by following partly a procedure 
adopted by Wiechert.'* This author, however, 
does not formulate the problem in our way, but 
rather assumes the deformation to be given and 
computes the force in some special cases. 

Let X; be the stress tensor and x; the deforma- 
tions, such that in Voigt’s?® notation in the 

28 L. Hartshorn, N. J. L. Megson, and E. Rushton, Proc. 
Phys. Soc. 52, 817 (1940). . 

2°W. Voigt, Lehrbuch der Krystallphysik (Leipzig, 
1928). X:, i=1, 2, 3 denotes the tension components, Xi, 
1=4, 5, 6 the shears and correspondingly for the x;. The 
elastic potential is in this notation 


, “J 
w(x1, X2° + -X—) =} Zz CijXiXj. 
i,7 
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limit {-0: 


6 
X:=¥D. c;jx;. 
j=1 


The c,;; are the elastic constants of the body . 


reducing to two independent ones in an isotropic 
material. If X;“ denotes the i-component of the 
internal stress experienced by species with a 
relaxation rate y“, we have 


N 


X= XM. (6) 


n=1 


Maxwell’s Eq. (1) for each species reads as 
follows: 


6 dx; 
=> Cc ij —— YX > n=1, 2---N. 
j=! dt 


dX ;™ 
dt 


The solution of this differential equation for 
q 
X;, if we consider the time interval 7;=t=T, 
is: 
= me - poe (n) o_o 
X M(T) =X (Tye y\"?(T—T}) 
6 T oe dx, 
$40 fo et (7-0¢,,—dt. (7) 


I=17T) dt 


Combination of (6) and (7) gives: 


N 
X ;( T) = =. X (Tye (2-1 
n=1 
- . r (n)omr dx; 
4X J ciyMer™t-9—at, (8) 
j=1 n=1 T) dt 


In this way we can find the stress tensor from a 
knowledge of the initial conditions, the tensor 
of the deformation rate and the magnitude of 
all y™. With Wiechert, we introduce a new 
function: 


N 
¥ij(u) => ce (9) 
n=1 


and find: 


N 
X(T) =D XM (Tye F-T 
n=l 
6 T 
+2. 
7=1 T1 


dx; 
¥i;(T—t)—dt. (8a) 
dt 


Equation (8a) together with the equations of 
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motion 
X(T) =X (T), 


where X;‘? stands for the i-component of the 
system of external stresses, determine the 
behavior of our viscous elastic body built up 
from N partial systems each obeying Maxwell's 
relation with a particular value of the relaxation 
time. 

In order to arrive at simpler results, let us 
consider an isotropic body, in particular a thread 
of length / to which a uniform tension X“ is 
applied, resulting in a_ relative elongation 
Al/l=x. Let E™ be the Young modulus of 
species n, then the modulus found at the begin- 
ning of the extension experiment 7,;=0, when 
all relaxation effects vanish, equals 


N 
E=> Em. 


n=1 
Equation (9) reduces to 


Vv a 


y(u)=>> E™Me ieee | go(yje“"dy. (9a) 


n=1 0 
The contribution of species » to the internal 
stress at (= 7,=0 will evidently be 


XM (0) =X (O)E™/E, 
In this way (8a) and (10) give with 
k(t) =dx(t) /dt, (10) 


T X“(0) 
xo(r)= f ¥(T —t)k(t)dt+- ¥W(T). (8b) 
0 


(8b) represents an integral equation either in 
¥(u) or in dx/dt. Let us first assume that we 
know the elongation-time curve. Then we have 
an inhomogeneous integral equation with a kernel 
k(z) proportional to the rate of elongation. The 
solution can be written in terms of Fourier 
transforms.*® The Fourier transform Ly(u) of 
¥(T) is defined for 0<7T< = as 


1 x 
Ly(u) =——— f V(y)e™dy 
(2x)} 0 


30E. C. Titchmarsh, Introduction to the Theory of 
Fourier Integrals (Oxford, 1937). 
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and therefore 


1 a 


f Ly(uje~** “du for T>0. 
(2x)! J_, 


¥(T)= 


We further write 7—‘t=z, T7=y, and obtain 


X‘°(0) d 
¥(y) =X“ o)- f ¥(z)k(y—z)dz. (8c) 


Let us consider first the equation in the limits 


0=y: Then forming the transform of p(y): 
X‘°(0) 
Ly(u) 
E 
—f |x wf R(y—2z (ds |evdy 
° 
1 £ Lx 
= DLyt(u) — . f Va)ds f et Oudt 
(29)? 0 0 
a = Lxi(u) — (24) Ly Li (u). 
hus 
Lxe(u) 
Ly(u) = 
X (0) 
+(27r)'L;(u) 
E 


and we find the solution 


” 1 “ Lx «©(u) 
¥(T) =- { e~'Tdu 
(27)! X‘(0) 


+ (2x)!L,(u) 
E 
for 7 >O0. 


In order to find the solution for finite limits 
we consider that the above integral equation 
retains the same form if we replace ¥(y), X“(y) 
and k(y) by the same quantities each multiplied 
by a factor e~*. If a is chosen as a sufficiently 
large positive quantity, 
before and extend the 
infinity. Now the 


V(y)e~™ is: 


then we can proceed as 
limit of integration to 
transform of the 


1 L 
: W(t)e-* ‘ei 'dt 
3)? 0 


— fv eiutiv) dt = Ly (w), 


function 


Ly(uje—vu = 


~ (2m)! 
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if w=u+iv. Remembering that L,(w) is the 
transform of y(t)e~*' we have: 


1 +e 
¥(y)= * J Ly(w)e— edu 
(29)? J_, 


1 a 
_ f Ly(w)e—4t du 
(29)? J_, 


1 ia+z 
= - f Ly(w)e "dw 
(2r)} ia—x 


and similarly for the other functions. In this 


manner we obtain the solution of our integral 
equation 


‘ 1 {3 ota Lx‘ (w) 
(y)= - ” lla 
(29)} J ia—x (0) 


- —+ (2 wr) Li(u 
- 


It can be shown* that this can also be written as 
E 


¥(y) = X(y) | 
X (0) 


—— i) X‘(s)m(y—sz)dz 
X6°(O) Jo -(11a) 


| 
L.(w) | 


1 iat+z 
m(y) =—— f ager e~ 'vedw. 
(29)? Jia-x 1—(2r)* Li (w) 


Having found y(y) 


with 


, we get the distribution of 
relaxation rates by solving the integral equation: 


x 


¥(y) = E{ eM o(y)dy. 
0 


This is Laplace’s integral equation*® with the 
solution 
¢g(y) = 


—f- V(ty)et*dy fory>0. (12) 


Equations (11) and 


(12) represent the final 
solution. 
Suppose now that we know in (8b) the 


distribution function and therefore ¥(y). By a 
reasoning similar to that applied above, we find 
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the rate of deformation as function of time :3° 


dx 1 oto L(w) - 
—(T)=— f —_——e~iTudy, 
d i 


(13) 
t 2 Ly (w) 


a—® 


if L,(w) denotes the Fourier transform of the 


function 
ae. ee 
Ree ee 


Equations (11), (12), and (13) thus contain 
the solutions of the respective problems: Given 
the deformation-time curve of an_ isotropic 
thread under tension, find the distribution of 
relaxation times or rates and, given the distribu- 
tion of relaxation times, calculate the elongation- 
time curve. 

Wiechert™ proposed for the distribution func- 
tion a logarithmic formula 


b Y max “| 
¢(y) =—— exp ; —57{ In —— } ¢. 
Jr ort 


There is no a priori reason for such a relation. 
For the orientational relaxation of polyvinyl 
chloride Fuoss and Kirkwood™ found an entirely 
different representation. At present we must 
consider the problem from the opposite point of 
view, that is, carry out extension experiments 
similar to those of Taylor and his collaborators" 
and Jenckel.*°*! Then formulating an integrable 
analytical expression for the elongation-time 
curve, we find the distribution of relaxation 
times or more correctly those parts of the 
distribution which contribute in a measurable 
manner to the total stress when the first reading 
is taken. The above experiments on glass and 
glassy selenium, respectively, can be interpreted 
on the basis of two relaxation mechanisms only. 
Taylor’s empirical expression is identical with 
Eq. (5) for two ;.%-*4 The results of these ex- 
tension experiments should be discussed from 
several points of view. For polar materials, it 
would be interesting to compare the distribution 
of mechanical and dipole relaxation rates. 
Furthermore, comparison of different distribu- 
tion curves may allow correlations with the 


a E. Jenckel, Zeits. f. Elektrochem. 43, 796 (1937). 
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conditions under which the samples were pre- 
pared, e.g. polymerization temperature and 
branching, polydispersity, previous treatment 
of textile fibers, and extent of crystallization. 

Kuhn!® and Frenkel and Obrastzov?® have 
pointed out the possibility of producing shear 
waves in media with finite relaxation times, if 
the wave frequency is large compared with the 
relaxation rate. In this case, the partial mecha- 
nism connected with this particular relaxation 
time or groups of relaxation times would behave 
like an elastic solid (infinite relaxation time). 
A viscosity of about 107-10° poises (glass 
blowing viscosity) corresponds, according to 
Maxwell’s relation between shear modulus, re- 
laxation time and viscosity, to a relaxation rate 
y ~10* sec.—', while a viscosity of 105-10° poises 
as appearing, for instance, in injection molding 
leads to y ~ 10° sec.—'. The limiting frequency of 
ultrasonics on the other hand, is of the order of 
10° sec.—!. In such media, therefore, transverse 
waves should experience practically no damping. 
Raman and Venkateswaran® showed, by investi- 
gating the nature of the scattered light in benzene 
and glycerine, that a solid-like behavior is ex- 
hibited for disturbances with sufficiently large 
frequencies. In a medium having a whole distri- 
bution of values, it may be possible in this way to 
draw some conclusions as to this distribution by 
measuring the dispersion of the absorption of 
shear waves. An extension of Isakovich’s, Leon- 
tovich’s and Frenkel’s and Obrastzov’s”® highly 
interesting investigations on the propagation of 
waves to a medium possessing a relaxation spec- 
trum would fall in this line. 

It may be noted again that the above con- 
siderations hold in the range of validity of 
Hooke’s law only. In the more general case the 
elastic moduli £; will depend upon the rate of 
flow. Entirely different phenomena will occur if, 
for instance, some of the species are brought to 
the limit of rupture. The author had the privilege 
to present this work before a meeting of members 
of the Research Laboratories of the Dow Chem- 
ical Company, Midland, Michigan. 


32 C, V. Raman and C. S. Venkateswaran, Nature 143, 
198 (1939). 


207 








Innovations in Instruments 








Rheostat Cages 


A variety of rheostat cages has been developed by the 
Ohmite Manufacturing Company for use with Ohmite 
rheostats. A cage enclosure is a convenient form of table 
top or surface mounting and is advisable where there is 
possibility of mechanical injury to the rheostat or human 
contact with electrically “live’’ parts. Surface mounting 
cages are widely used in connection with fractional h.p. 
motor controls; heating device controls; dental and medical 
equipment; and for use in laboratories. 





Ventilated cages with perforated sides are most generally 
applicable. Other types available are: dustproof cages, 
cages to house rheostats in tandem, explosion proof cages, 
ventilated cages with one-half closed as a splash guard, 
and a laboratory type semi-enclosed table-top cage. Various 
types of terminals can be supplied. For details write to the 
Ohmite Manufacturing Company, 4835 Flournoy Street, 
Dept. 7C, Chicago, Illinois. 


Index for Chemical Analysis by X-Ray Diffraction 


A 4000-card file index of x-ray diffraction data for use 
in the Hanawalt method of chemical analysis by x-ray 
diffraction has recently been published by the American 
Society for Testing Materials. This compilation is spon- 
sored by a joint committee of the A.S.T.M. and National 
Research Council, under the chairmanship of Professor 
Wheeler P. Davey of Pennsylvania State College. The 
data include not only Hanawalt’s original published data, 
with his later corrections, but also additional data that have 
been contributed by the Aluminum Company of America 
and the New Jersey Zinc Company, together with data 
taken from the technical literature in the English language. 

The cards give all pertinent data found in the sources 
with provision for insertion of accessory data such as 
crystal structure, density, etc. The index identifies the 
three strongest lines in the x-ray diffraction pattern of 
some 1300 crystalline compounds, the chemical names and 
symbols of which are as given by the various sources. 

The Hanawalt method has been described in the tech- 
nical literature, and an A.S.T.M. Committee is perfecting 
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a tentative recommended practice for the identification of 
crystalline materials by the x-ray diffraction method, with 
early publication anticipated. 

Copies of this 4000-card index packed in finished con- 
tainer boxes can be obtained from A.S.T.M. Headquarters, 
260 South Broad Street, Philadelphia, Pennsylvania, at $50 
per set. 


New Rotary Type Control Switch 


A new double-break, direct-acting rotary switch for the 
control of circuit breakers has been developed by the 
Roller-Smith Company, Bethlehem, Pennsylvania. The 
switches are insulated for 600 volts and have a continuous 
current carrying capacity of 10 amperes. 

This switch is built up of units, each unit consisting of 
two single break contacts connected in series by means of 
a jumper. The switch has a spring return arrangement so 
that when the handle is released the switch returns to the 
“off”’ position, opening the contacts. The handle is of the 
pistol grip type, facilitating operation of the switch. 

A red and green target is furnished as a part of each 
switch to indicate the last position of the switch. For ex- 
ample, to operate the circuit breaker control switch in 
conjunction with indicating lamps, the switch is turned to 
the ‘‘Close” position, 45° clockwise. This operation closes 
the breaker and energizes the red lamp through a breaker 
auxiliary switch. When the handle is released, the spring 
return feature causes the switch to return to the “Off”’ 
position, opening the contacts. The red target shows that 
the last position was to close the breaker. To trip the 
breaker, the switch is turned to the “Trip’’ position, 45° 
counterclockwise. This brings up the green indicator and 
lights the green lamp through another auxiliary switch. If 
the breaker trips by some means other than by the switch, 
the green lamp and the red target are seen at the same 
time. This indicates that the breaker tripped on some fault. 


Fibronized Koroseal Tubing 


Irvington Varnish and Insulator Company announces 
their new Fibronized Koroseal Tubing designed for use as 
electrical insulation by manufacturers in the electroni¢, 
instruments, aircraft, electrical appliances, and power in- 
dustries. This new insulation was developed from Koroseal, 
a product of the B. F. Goodrich Company. 

The advantages of Fibronized Koroseal include inside 
and outside smoothness, exceptional elasticity and close 
manufacturing tolerances; excellent resistance to acids, 
alkalis, solvents, and heat. It is fireproof and possesses an 
insulation resistance of infinity after 16 hours at 90 percent 
R.H. and 105°F. As a result of the A.S.T.M. test of being 
subjected to 225°F for approximatély 1000 hours, it was 
found to have retained its flexibility. It has a tensile 
strength of 2845 pounds per square inch, a dry dielectric 
strength (0.022” wall thickness) of 1050 vpm, a wet di- 
electric strength (0.022” wall thickness) of 817 vpm after 
24 hours’ immersion. 

Fibronized Koroseal Tubing meets or excels all A.S.T.M. 
specifications. It is available in A.S.T.M. sizes and in a 
variety of colors including a transparent finish. 
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